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On the Uniqueness and Non-existence of Stokes Flows 


ROBERT FINN & WALTER NOLL 


Communicated by C. TRUESDELL 


1. Introduction 


It has been known since the time of Stokes that under suitable assumptions 
on the behavior of the velocity field at infinity, there exists no steady two- 
dimensional flow of a viscous fluid past an obstacle, in which the velocity com- 
ponents are infinitesimals of the first order. On the other hand, flows of this 
type in three dimensions are explicitly known (cf. [2]). An explanation sometimes 
advanced for the discrepancy in results is that in the known flow past a sphere [2], 
the assumption on the ratio of inertial to viscous forces which is used to derive 
the equations of motion becomes violated in a neighborhood of infinity. In the 
view of the authors, this phenomenon brings into serious question both the 
physical and the mathematical validity of the known results, and makes it 
imperative to investigate the sense in which the boundary problem is correctly 
set. More precisely, it is natural to ask whether there exist two-dimensional 
flows in which the velocity tends to its limit more slowly than has in the past 
been assumed, and whether a three-dimensional flow is unique in a (physically 
reasonable) class of flows for which the usual uniqueness proofs break down*. 
In this paper we prove that the answer to the first question is no, to the second 
yes. In fact, we show that there are no two-dimensional flows for which the 
velocity is bounded, and we prove the uniqueness of a three-dimensional flow 
under the single assumption that the velocity tends to a limit at infinity **. As 
a corollary, we obtain representations, valid in a neighborhood of infinity, for 
the velocity field of the most general such flow. 


Once the uniqueness of a flow past an obstacle with prescribed velocity at 
infinity is established, the continuous dependence of the velocity field on the 


* In [2], SOMMERFELD dismisses this question with the remark, ,,DaB die so 
gewonnenen Gln. (13) und (13a) die einzig méglichen Lésungen unseres Problems 
sind, haben wir zwar nicht bewiesen. Wir folgern es aber aus dem Axiom, da jedes 
richtig gestellte Problem der mathematischen Physik nur eine Lésung haben kann.“ 
This reasoning is evidently circular. In order to establish that the problem is cor- 
rectly set it is necessary to prove the uniqueness of a solution. 

*xx After preparation of this manuscript our attention was directed to a paper of 
CHARNES & KRakowskI, Carnegie Inst. Tech. Technical Report No. 37 (1953), in 
which a rigorous proof is given for the non-existence of a plane Stokes flow. The 
proof we present here seems simpler, and applies without essential change to the proof 
of uniqueness of a three-dimensional flow and to the derivation of asymptotic 
representations. 
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limiting speed is a direct consequence of the linearity of the equations of motion. 
The existence of such a flow, despite extensive studies in this direction, is yet 
to be proved. 

2. Notation and definitions 


Throughout this paper we use G1BBs’ notation of vector analysis. The under- 
lying space is of dimension 7 = 2 or n =3 (plane or ordinary three-space). Points 
are denoted by the letters P, Q, etc., vectors by bold face letters. By P+r 
we mean the end point of the vector r if it starts at P. The symbol e is reserved 
for a variable unit vector or, equivalently, for a variable point on the unit sphere 
(circle) w. The surface element (element of arc) of the unit sphere (circle) will 
be denoted by dw. In the plane (n= 2) wxv and curl u have to be regarded as 
scalars. 

By a slow flow, or Stokes flow, in a region, we shall mean a vector field q(P) 
which satisfies the system 


(2.1) Aq=gradf, divg=0 


for some scalar field #(P). In the text we shall refer to a Stokes flow simply 
as a flow. We assume that within the flow region q is three times and # twice 
continuously differentiable. # is determined only up to an arbitrary additive 
constant. 

By an obstacle 8 we shall understand a finite number of piecewise smooth 
nonintersecting simple closed surfaces (curves). A solution of (2.1) defined 
throughout the exterior of an obstacle 8 will be called a flow past 8 if q is con- 
tinuous up to 8 and q=0 on ¥ and if # and the first partial derivatives of q are 
bounded up to &. 

The region exterior to a sphere (circle) of radius 7, and center Q shall be 
denoted by € and called a neighborhood of infinity. Let v=v(P) be a vector 
or scalar field defined in ©. Then we write v=o/(r*) if limrv—*v(Q+r7e) =0 
uniformly in e, and we write v=O(r*) if |y~*v(Q+,7e)| <const.< o for large 7. 
These properties are evidently independent of the choice of Q. 


A flow defined in a neighborhood of infinity will be called uniform at infinity 
if q=q)+o(1), and q is said to be its velocity at infinity. 


3. The basic lemmas 


Lemma 1. Leth be a field which is harmonic in a neighborhood of infinity ©. 
Then h has a unique decomposition h=h,+hg, 

where h, 1s harmonic in the entire space (plane) except at co, and 

where hy is harmonic in © and, for n=}, hy=O(r>) while for n=2 it is of 
the form 


(3.1) hy(Q + re) =alogr+h}(Q4+re), h¥ =O(r2)*, 


h is harmonic at oo if and only if h=O(r—"*?), 


* hy, h, and a are uniquely determined by h and independent of the choice of Q, 
but hF, of course, is not. 
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Lemma 2. Let h be as in lemma 1. Then the following three statements are fully 
equivalent: 


(a) ese eee eco PS 2 
‘ ==TA OG ae if N= 3 
(b) (DON ge 
(c) h has an expansion of the form 
(3.2) h=h(Q+ re) => r* s,_, fe) *, 
k=v+1 


uniformly convergent in every compact subregion of ©. 


The serves obtained from (3.2) by taking the gradient term by term is also uniformly 
convergent 1m every compact subregion of © and hence it represents an expansion 
for grad h. We have gradh=O(r~’—*) and, conversely, if grad h=O(r-’—*), then 
h=O(r-*-*) + const. 

These two lemmas are corollaries of well known theorems on isolated singu- 
larities of harmonic functions (cf. [7] §§ 93—95). 


Lemma 3. Let h be as in lemma 1. Then the following two statements are fully 
equivalent: 


(a) b= Oy"), 
(b) fh(Q+re\dw=0 for any point Q and any r>1,**. 


Proof. Assume first that (a) holds. Then by lemma 1 / is harmonic at ov. 
Hence the function h(P)=A(Q+se)=s~"t?h(Q+s te) is harmonic in the 
neighborhood of Q and at Q. Therefore by the mean value theorem for harmonic 
functions we have (w = 22 or 421) 


wh(Q) = fh(Q+ re) dw =r" fhn(Q+re) do. 


On the other hand, (a) implies 4(Q) = lim7”~?h(Q+yre)=0, hence that (b) 
holds. ane 


Assume now that (b) is true. We use the decomposition h =h, +h, of lemma 1. 
Since hf, is harmonic everywhere the mean value theorem implies wh,(Q) = 
fh.(Q+re)dw. Using the same method as before and lemma 1, we get 


fhg(Q + re) do = 


@ 


ees ton 2 


wbrt for 4 =%3), 
where b= limrh,(G+yre) for »=3. Therefore (b) implies 


22a logy for #=2 
by for = 3: 


—i4(Q) -| 


* For n=2 the s,(e) are trigonometric functions while for ~=3 they are 
spherical harmonics of order k. 
xk This lemma, as well as its proof, is valid also for n>3. 


7* 
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These relations are possible for arbitrary Q and arbitrary large 7 only if h,(Q) =0 
and a=b—0. Hence h=h, and, for n=3, hg=0(r1). Lemma1 for »=2 and 
lemma 2 for »=3 show that (a) holds, g. e. d. 

Lemma 4. Let q be a three-dimensional flow defined in a neighborhood of 
infinity & and uniform at oo. Then 


(3.3) curlgaeOG =); 


Proof. We put w=curlg. By (2.1) we have 4w=curl dq=0, 1.e. w is 
harmonic in ©. We denote the sphere of radius 7 around Q by G, and the volume 
between ©, and ©, by &,,,. Then Green’s identity yields 


_ few ow 
0= [ dwav {FZ as 4 [ #as, 
Sr 


B, To Sro 


where 0/dn denotes the normal derivative in the direction away from &,,.. The 
first integral on the right is independent of 7 and will be denoted by a. For the 
second we have we = ee w(Q+re) and dS=rda. Thus 


(3.4) a [w(@ +re)do = —ar. 


Integration with respect to 7 yields 


(3.5) fw(Q+re)dw=ar+b. 
We multiply this by 7? and integrate again with respect to 7, obtaining 
(3.6) fwava=Fat tote, 
2 
Bror 


where € is independent of 7. Another Green’s identity yields 


(3.7) feulgadV=fqxndS+4+fqxndS, 

ror ro Sr 
where nm is the unit normal vector directed away from &,,. Again, the first 
integral on the right is independent of 7 and will be denoted by d, the second 
can be written in the form 


faxndS=P7 fq(Q+re) xed. 


e 


Hence from (3.6) and (3.7) we get 


(3.8) fale+re) xedo = ‘a4 b+ (e—ad)r*, 


By hypothesis, q is uniform at ov, i.e. lim g(Q-+-re) = uniformly in e. Hence, 
letting 7—> co in (3.8) yields b=0 and qyx fedw=0=4a. From (3.5) and 


lemma 3 for 7=3 follows then w=O(r-), q. é. a. 
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Lemma 5. Let q be a plane flow defined and bounded in a neighborhood of 
infinity ©. Then 


(3.9) CMG 0+) 
The proof of this lemma is analogous to that of lemma 4. Surface integrals 


are to be replaced by line integrals and volume integrals by area integrals. 
w=curlq is now a scalar field. In place of (3.4) we get 


a [#0 +re)dw=-—ar, 
so that 2 
(3.10) fw(Q+re) dw =—alogr+b. 


Instead of (3.8), we now have 


fal@+er) xedo = ar(5 logy z)+ 52 L (c —d) 77, 


It suffices here to assume that q is bounded in order to conclude that a=0 
and 6=0. From (3.10) and lemma 3 for »=2 follows w=O(r4), q. e. d. 


4. Uniqueness of three-dimensional flows 


Theorem I. For any obstacle 8 there is at most one three-dimensional flow q 
past 8 which 1s uniform at infinity and has a prescribed velocity at infinity qo. 

Proof. Since the equations (2.1) are linear, it is sufficient to prove that 
Qo =0 implies q=0. 

By lemma 4 we have 
(4.1) OG =O) 7) 


and hence, by lemma 2, 
— curl w = — curlcurl g = Aq = gradp = O(r'4). 


From lemma 2 follows that =O(r-?) + const., but since we are free to adjust p 
by an additive constant, we may assume 


(4.2) p=O(r?). 

Consider now the Green’s identity 

(4.3) f (curlg)?dV=fq-AqdV + fcurlq:-(qxn)dS, 
B B S 


* Professor TRUESDELL has pointed out to us that our proof of non-existence 
of plane flows can be viewed as an application of a theorem of HAMEL and KampE 
DE FERIET, which asserts that in a plane motion adhering to a fixed boundary the 
vorticity is orthogonal to every harmonic function. But in a slow plane motion the 
vorticity is harmonic. Therefore the vorticity is zero, from which one concludes 
easily that the velocity is zero. To make this argument rigorous in our case it is 
necessary to find a suitable estimate on the behavior of the vorticity at infinity. 
This estimate is provided by lemma 5. 
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valid for any smooth vector field q with div q=0 defined in a finite region B 
with piesewise smooth boundary S. We apply this identity to the flow q and 
the region B, between the obstacle 8 and a sphere G, so large that % lies interior 
to G,. By (2.1) we have q-4q=q- grad p=div(pq), so that (4.3) can be 
transformed into 


(4.4) PL maps io —curlg:-(q xn)]dS. 


Here we have made use of the fact that q=0 on B%. It follows from (4.1), (4.2), 
and g=o(1) that the surface integral in (4.4) tends to zero as roo. Hence 


(4.5) ep (curl q)?dV = 0, 


where ¥ is the region exterior to 8. Since curl q is continuous, it follows that 
curl g=0 throughout &. Together with div q=0, this implies 4qg=0. But 
q=0 on % and q=o(1); hence, by the maximum and minimum principle for 
harmonic functions, applied to each component of q, it follows that q=0, g.e.d. 


5. Non-existence of plane flows 


at 


Theorem II. For any obstacle 8, any bounded plane flow q past % 1s the 
state of rest q =0. 

Proof. An argument analogous to that used in the proof of theorem I shows 
that from lemma 5 follows 


(5.1) w= cul G¢-= O04), PSO) 

Also, we have 

(5.2) J (curlg)?dA =f [Pq-n—curlq(q xn)]ds, 
I, G, 


where Y, is the region between % and a large circle ©,. But here we can conclude 
from (5.1) and the boundedness of q only that the integral 
(5.3) I= f (curlg)2dA = fwd A 

U U 


over the exterior YU of 8 is convergent. However, we know that w is harmonic 
in % and that w=O(r1). Hence, by lemma 2, (3.2), w is of the form 


(5.4) w=w(Q+re)=r's(e)+O0(r%), 
and therefore 
(5.5) w® = 7-255 + O(r-3), 


Let Y,,, be an area between two circles C,, and ©, both exterior to B, so that 
Y,,, is contained in YW. We then have from (5.5) 


fwd = log - [ do + for dA. 


Wrox w Tae 


This must tend to a finite limit as 7-> co, because the integral (5.3) converges. 


The second term on the right certainly converges and hence, since log ”~ is 
Y 


not bounded as 7 , we can conclude that f{s2dwm=0 and hence 0 


w 
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Therefore, by (5.4), w=O(r-’), from which, as before, =O(r-2) follows. Going 
back to (5.2) we see now that the integral on the right in fact tends to zero as 
y—> co and hence, as in the proof of theorem I, we conclude that curl qg=0, 
div qg=0 and 4q=0 throughout Y%. Thus q is harmonic in a neighborhood 
of infinity. Since q is also bounded, it follows from lemma 1 for »=2 that q 
is harmonic at oo and hence uniform at infinity. Therefore, for some choice 
of a Cartesian coordinate system, we have g,—9), g,—>0, where qg, and q, are 
the components of q. As before, the maximum principle for harmonic functions 


09x cE; 
oy Ox 
But, since g,=0 on ¥, we must have g,=0 everywhere, q. ¢. d. 


implies g,=0. Therefore curl q= =0, divq= = 0, and hence g, =const. 


6. Representation formulas 


The following theorem establishes a one-to-one correspondence between three- 
dimensional flows with given velocity at infinity and vector fields which are 
harmonic at infinity. 


Theorem III. Let q be a flow defined in a neighborhood of infinity © with 
velocity at infinity Gy. Then there is a harmonic vector freld u, defined in © and 
harmonic at infinity such that 


(6.1) Gq=Qtu-+t grady, 
where 
(6.2) y=y(Q+re) =—iffstdivu(Q+se) ds. 


Conversely, if u ts an arbitrary harmonic vector field defined in © and harmonic 
at infinity, then (6.1) and (6.2) give a flow with velocity at infinity qo. 
Proof. Consider the differential equation 


(6.3) A (ry) — Ay = 64 + 47 4 = p(Q +12). 
A special solution is 
(6.4) 4=1(0 tre) =47 J sip +se)ds. 


The integral converges because of (4.2) ; moreover, 
(6.5) Wa Oat) Aye Oll7e2), Segtadi(y2y\i== O77): 
Applying Laplace’s operator to (6.3) and observing that 4f=0 we get 


a) — 
1447+ 4r ae Ay =0, 


hence Ay=c(e)r* +. This is compatible with (6.5) only if c(e)=0. Thus 4y=0 
and, by (6.3), A (727) =p. Hence A [q — q — grad (r? 7) ] =grad p — grad A(r?y) =0 
i.e. U=Q—Q,)—grad(r27) is harmonic. From (6.5), ¢Y—Q)=o(1), and from 
lemma 1 it follows that w is harmonic at infinity. Also, we have 


div u = — div grad (7? ¥) = — A(r?v) = — #, 
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hence y=7?y in (6.2). Thus the first part of the theorem is proved. The second 
part can be verified immediately. 
By use of the series expansion (3.2) for wu, 7.¢. 


(6.6) u= Yrs), 


(ak 


it follows easily as a corollary of theorem III that any flow which is uniform 
at infinity can be represented by an expansion 


— 1 
(6.7) q=H%+>,7% Is.-s(€) — Bak — ty Set (2) , 

h=1 
where the spherical harmonic vectors s;(e) are defined in terms of the s,(e) by 
(6.8) Si(e=7') grad div [z7"h* $ze9le) a are 


Conversely, for any set of spherical harmonic vectors s,(e) such that (6.6) con- 
verges uniformly in some neighborhood of infinity, (6.7) and (6.8) define a flow 
with velocity at infinity q,. The classical Stokes flow past a sphere of radius 
7, [2] corresponds to 


Sy = — 270%» $,=$7a[G— 3€(Q°2)], S$, = 'S3= 85 +"-= 0: 


Theorem IV. Let q denote a flow defined in a neighborhood of infinity © with 
velocity at infinity Qo. Then there are harmonic vector fields 84 and H, defined 
in © and harmonic at infinity, such that 


div H =0 
(6.9) ee 
div. A (7? $2) == 0, oS Ota 
and such that 
(6.10) G= Qt curl (r? Q) + H. 


Conversely, if H and 82 are harmonic vector fields defined in ©, harmonic at infinity 
and satisfying (6.9), then (6.10) represents a flow with velocity at infinity qo. 
Proof. Define 


(6.11) Q = —Zr i fsscurlg(Q+se)ds. 
Then, as in the proof of theorem III, 4Q=o0, A (r? $2) = —curl q, 


A(q — curl? Q) = Aq + curlcurlg = 0, 
hence 


Gq=Qtcul?fQA+ H 


for a harmonic vector H. Since curl q=O0(r), we have Q =O (7) and curl 7? Q = 
O(r+). Hence H =o(1). By lemma 2, Q and H are harmonic at infinity. 
If, conversely, H and Q are harmonic in & and at infinity and satisfy (6.9), 
then if we set 
G=Q+tcul?fvO+ HA, 
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then div g=0 and —curl dg =4AA?r?Q=0, hence there is a scalar p such that 
Aq=grad p. Also, one can see that q=q,+0(r), g. ¢. d. 


7. Representation of a flow past an obstacle 


If the flow is defined throughout the exterior B of an obstacle 8 and q=0 
on & the results of the preceding section can be made more specific. For we 
may set 


(7.1) @(0) = -4 fr curl (P) dV 
B 


where 7 is the distance between Q and P. This integral in general does not 
converge absolutely, but it exists in the sense of a principal value as the limit 
of integrals over concentric spheres with center at Q or at any arbitrary point. 
This is the case because by (6.1) and (6.6) we have 


fr curlG =-7 -curlu—=7 *s, x e-- O(7 3), . S$, = const., 


and fedw=0. Since n- curlg=0 on &, one can show that A@ =curl q and 


div p=0, hence —curl curlee = Ag =curl q, which implies the existence of a 
scalar h such that q=q,)—curlp-+gradh. The function / may be multi-valued 
if GB is not simply connected, but gradh’ is single-valued in B. Evidently, 
div gradh=Ah=0. Also, simple estimates on (7.1), using curl q=O(r%), 
show that curl =O(r1). Hence gradh=O(r4), and, for any given deter- 
mination of / in a (simply connected) neighborhood € of infinity, 4 =O (log 7). 
By lemma 2, / tends to a limit at infinity in ©. Thus follows 


Theorem V. Every flow past 8 with velocity at infinity q, admits the repre- 
sentation, valid in the exterior & of %, 


(7.2) G=Q— cule + gradh, 


where ep is defined by (7.1), and where his harmonic in B. If B ts simply connected, 
h is single valued and may be chosen to be harmonic at infinity. Otherwise, any 
determination of h tends to a finite linut at infinity. 

We may combine this result with the results of §6 to obtain still other 
representations. 

Theorem VI. If % is star-shaped relative to a point Q, every flow past 8 which 
is uniform at infinity admuts the representation 


(7.3) q=q(0+e7 =qtcul/’?Q + curl H+ gradh, 


where H and hare harmonic in & and at infinity and Q is defined by (6.11). If B 
is not star-shaped, (7.3) are valid in any star-shaped neighborhood of infinity in 
the flow. 

Proof. Tf ¢~ is defined by (7.1), Q by (6.11), then A(p+7?7Q)=0. Hence 
= —(7?Q+H) for a harmonic vector H. Using (7.2), 


G=Q% + curl’? Q + curl H+ gradh. 
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From (6.9) and (7.1) we find 7? —=0O(1), p=O(log 7). Hence H=O (log 7) and, 
using lemma 2, we may assume H is harmonic at infinity. This yields (7.3),. 
For the classical Stokes flow past a unit sphere with center Q we obtain 
for q(Q+er7) 
ae 3 2 % : Qo 
Cire Deri curl (r curl ze — g curl (curl y i" 
3 


By Mee 2 Ch ede eq 
=I, curl (r curl * ‘i grad ( 7) )- 
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Stress-Deformation Relations for Anisotropic Solids 


G. F. SmitH & R.S. RIVLIN 


1. Introduction 


We consider a body of material to be subjected to a deformation in which 
a point of the material initially at X; in a rectangular Cartesian coordinate system 
x; moves to x; in the same coordinate system. In the theory of finite elasticity 
it is assumed that there exists a strain-energy function W which is a polynomial 
function of the displacement gradients 0x,;/0X; only. It can then be shown that 
W is a polynomial function of the quantities g;; and /g defined by 


Sea, Cea: . 
87 aX, aX? g =|8i,|- (T=1)) 


If the materialis isotropic, then the dependence of W on the displacement gradients 
can be further restricted to polynomial dependence on J,, J, and g, where 


I, = 83, Te, = 84; &js- (1.2) 


If the material is anisotropic, then W can be expressed as a polynomial in a 
number of polynomials in g;; each of which is invariant under the transformations 
characterizing the class of anisotropy. These invariants have been obtained for 
each of the classes of crystal symmetry by SmiTH & RIVLIN (1957). 


In the present paper we shall assume that no strain-energy function exists 
for the material but that the components of stress ¢;; in the coordinate system ¥; 
are polynomial functions of the displacement gradients 0x,/0X, and shall in- 
vestigate the way in which the particular symmetry of the material restricts 
the manner in which the stress components can depend on the displacement 
gradients. The calculations are carried out only for the monoclinic and rhombic 
systems. However, the corresponding calculations for the remaining crystal classes 
can be carried out by a similar method. The appropriate procedure for each 
class can be readily deduced from that employed by SmitH & RIVLIN (1957) in 
the case when a strain-energy function exists. 


2. Invariance under Rigid Body Rotations 
We assume that the stress components ¢;; at a point of the body, in the co 
ordinate system x,;, are polynomial functions of the displacement gradient 
0x,/OX,, thus 


1 


te; = t;;(0%,/0X,). (2.1 
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We now consider that the body undergoes, in addition to the deformation con- 
sidered, an arbitrary rigid body rotation in which points at x; move to %;, where 


and A;, are constants satisfying the relations 
AjAy, =A; Ane = Oi © \Ayj| =T. (2.3) 


The stress é; , in the coordinate system 4%; due to the resultant deformation is, 
of course, given by 


t; , = t,;(O%p/OX,). (2.4) 
Let x* be a rectangular Cartesian coordinate system which rotates with the 
rigid body rotation described by (2.2). The stress in this coordinate system is 


t;,(@x,[0X,) and is related to ¢,, by 
t,, = Aj, Ajrter- (2.5) 
From (2.1), (2.2), (2.4) and (2.5), we see that 
1,,(0%p/0X,) = Aggy Ajitys(O%)/OX,)- (2.6) 


The restrictions imposed on ¢;; by the requirement that it satisfy (2.6) for all 
A,, satisfying (2.3) have already been examined by GREEN & RIVLIN (1957) in 
a slightly different context. Following an argument similar to that previously 
used, it can be seen that ¢;; must be expressible in the form 


Ox, 0%; 


ome i Wig ke pee oe: 
tis f 6;; a ax, ax, FASS) (2.7) 
where f and /,, are polynomials in g,, and 1 defined by 
TONER OOH: _ | &%p bs 
oS OX, OX, and Y= lax, | (2.8) 


and /,, is symmetric for interchange of y and s. The relation (2.7) may be re- 
written in the form 


ax, O%; 1 
a a (ax aa Srs 6i;] fess (2.9) 
where /* is defined by 
[eat fore By lak (2.10) 


and is thus a polynomial in g,,. 


3. The Effect of Symmetry 


In this section we shall consider the further restrictions imposed on the 
relations (2.9) if the material possesses certain types of symmetry in its unde- 
formed state. 

We shall consider that the symmetry of the material is defined by three 
preferred directions in the undeformed material, given by the unit vectors h,,h, 
and h;, and the group of transformations which specifies the equivalent positions 
of the triad of vectors. The stress-deformation relation (2.9) must be form- 
invariant under each of the transformations of this group. 
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We shall consider only symmetries for which the associated group of trans- 
formations is a sub-group of the full orthogonal group, so that if S (=|| s.,|I) 
is a transformation of this group, 


Sig Sip = S55 Spi = 0; z- (3.1) 
Let x; be the coordinate system into which the system x, is transformed by 
this transformation. Then, 


45, = S55 x; and GENT OGD (3.2) 


where x; and X;, are the coordinates in the system %,; of the points at x, and X, 
respectively in the system x;. We have 


Sj TOF. = Si1Sjm Elm: (3.3) 
f 
Also, the components of stress in the coordinate system %; are given by 
ti; S77 Sim Lim: (3.4) 


The form-invariance of (2.9) under the transformation S requires that 


Ox, 


: : Bre ij) fre Boe?) (3.5) 


ti; =f (Bbq, 7) 0:54 Gee 


where 7 =| 0%,/@X,|. From equations (3.1) to (3.5) we obtain 


a Ox, ° Ox; 1 Fie he 
Ue =f* (Bp qT T) 0; Sh lax aX, 2 Emn 6,,) Sim Ssn Te hones T) ‘ (3.6) 


Comparing equations (3.6) and (2.9), we have 


Ox; O%; 1 2 
ip (2) 0;; lee aX, 3 Smn 6,;) Sip Serato ea) 


Ox; Ox; 
= f* (pq) 0 a7 Le ae = brs S:j| fre Coe) (3.7) 


We have omitted the arguments t and 7, since t=7 for any orthogonal trans- 
formation and we may understand polynomial dependence of /* and /,, on t 
on both sides of equation (3.7). From (3.7), we have 


f* (Bp q) = F* (Bp) (3.8) 
and 
(4) i Ox; 1 _ y) j Ox; 4 
ae aX, = 3 Smn 6,;) Sinden tre (ee) — Ge OX, — 3 &rs 84] frs(Bp4)- (3.9) 


Equation (3.8) implies that /*(g,,) isa polynomial scalar invariant of g,, under 
the transformation S. 

Let 0,;; and bi be the components in the systems x; and %; respectively of 
an arbitrary second order tensor. We define the second order tensor with 
components B,,; and B;; in the systems x; and %, respectively by 


(3.10) 


OX; Ox; 1 
B= bis (See ox. — 3 bro Oi 
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We note that B,, is not, in general, symmetric. Then, with the notation 


DEC SPS Bel 8 (3.11) 
equation (3.9) becomes = | 
F =B,,f,s(Bpq) = Brsfrs (8pa) - (3.12) 


This expresses the fact that F is a polynomial scalar invariant under the trans- 
formation S. It can therefore be expressed as a polynomial in a set of polynomial 
scalar invariants of B,, and g;;. Since F is linear in the elements of B;,, the 
only invariants involving B,; which need be included in the set are those in- 
volving B,, linearly, and these will occur only in the first degree in the expression 
for F. The form for /,,(g,,) can then be obtained from the expression for /’ by 


the relation * E eF 
Licleee) ae OB. = (3.13) 


4. The Monoclinic System 


If the material considered has monoclinic symmetry, the unit vectors h, 
and h,, defining two of the preferred directions in the undeformed material are 
not at right angles, while the third unit vector h, is perpendicular to the plane 
defined by h, and h,. We choose as our reference system the rectangular Car- 
tesian coordinate system x;, the x, axis of which is parallel to h,. The axes x, 
and x; may be in arbitrary perpendicular directions in the h, h, plane. The 
monoclinic system contains three classes which, together with the symmetry 
transformations defining them, are: domatic (I, R,), sphenoidal (I, D,), prismatic 
UI, C, R,, D,), where 


12(4,4,1),. C=(—4, 49a) OS 


D, = “y— 1, = 4); (4.1) 

and (a,b,c) denotes the diagonal matrix with diagonal elements a,b and c. 
We see from §3 that the stress ¢;; is given by equation (2.9) where /* is a 
polynomial in t and a polynomial scalar invariant of g,, under the transformation 
group defining the symmetry properties of the material. Thus, for the mono- 
clinic-domatic material, it must be a polynomial scalar invariant under the trans- 


formation R,. It is, of course, automatically invariant under the transformation 
I. With S=R,, we see from (3.3) that 


and 


P* (B11, 822) S38» as» Sar> Sie 7) = F* (811) Goo» Sas» Sas» (S515 == has Te ee) 


It follows from a theorem in classical invariant theory (see, for example, WEYL 
1946, p. 36 et seg.) that the necessary and sufficient conditions for (4.2) to be 


satisfied is that /* shall be a polynomial in 
811» 822) 833» 23, 831» S12, 8s kia) T- (4.3) 


* Note added in proof (September 1957). The result (3.13) can be obtained in the 
following manner. From (3.12) and (3.10) 


CH Ok, A OF oF 0B Oni, Oz Gl OF 
pia Si Eg eS SN One pe eS rs a i a 
(ax ey ee al Irs (8pq) 6b;;  0B,, 0b;; (ox exGe 6 brs 6) BD ras 


Without loss of generality, we may therefore take f,s im (2.9) to be given by (3.13). 


, 
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We also see from §3 that F, defined by (3.12), must be a polynomial in rt 
and a polynomial scalar invariant of B,, and &pq under the transformation R,. 
This requires that 
#611, 82a» bes. Bes, G31» £12, P11,P as, Bes, Bos, Bar, Bro, Baa, By3, Byy,T) 

= F(8115 82028939 023, = 8a1n —€19> Pir, Bass Bys, Bos, 
eT Ripe ety hme el ose Whee POPE nea) (4.4) 
As in the case of equation (4.2), together with the consideration that F is linear 


in the elements of B;;, it follows that F is a polynomial in the quantities (4.3) 
together with 


By, Boo, Bss, Bog, Bao, €31 Bai, 831 Bre; 812 Bai, S12 Pro: 
831 Bis, 831 Bar, 812 Pig, 812 Bar- (4.5) 


The quantities (4.5) can, of course, enter only linearly into the expression for F. 
We see from (3.13) that 


A141» 812% + 831 Pra, S12%13 +b £31 P13 || 
ll f-sI| = || S1e%12 + 831 Bio, Ago, Xo3 , (4.6) 
£12%13 + S31 Pia, Ao, X33 | 


where the «’s and f’s are polynomials in the quantities (4.3). We thus have the 
result that, for a material with monoclinic-domatic symmetry, the stress /;; is 
given by equation (2.9) in which /,, is given by (4.6) and /*, the «’s and f’s are 
polynomials in the quantities (4.3). It is readily seen that this result is also 
valid for the remaining classes of the monoclinic system since /* and F are auto- 
matically invariant under the transformation C and (4.2) and (4.4) are the condi- 
tions for their invariance under the transformation D,. 


5. The Rhombic System 


If the material considered has rhombic symmetry, the unit vectors h;, de- 
fining the preferred directions in the material, are mutually perpendicular. We 
take as our reference system a rectangular Cartesian coordinate system x;, the 
axes of which are parallel to the vectors h;. The rhombic system contains three 
classes which, together with the symmetry transformations characterizing them, 
are: rhombic-pyramidal (J, R,, R,,D,), rhombic-disphenoidal (I, D,, Dz, Ds), 
trhombic-dipyramidal (J, C, R,, R,, R,, D,, D,, D3), where I, C, R, and D, are 
defined by (4.1) and R,, R,, D,, and D, are defined by 


R,= (1, —1, 1), R, = (1,1, — 1), 
D,=(—1,4,—1) and D,=(—1,—1,1). (5-4) 


It is readily seen that for each of these classes, the restrictions imposed on 
f* by the symmetry of the material are 


{* (Gog, 831» £12) =1*(Se3, — 81> — 819) 
= f*(— go5, 831, — 812) = 1*(— Sea, — 81> 819): (5. 


b& 


NP G. F. Smiru & R. S. Rrviin: Anisotropic Solids 


Similarly, the restrictions imposed on /’, defined by (3.12), by the symmetry 
of the material are 
F (G3, S31» 812, Bes, Bsr, Biz, Bse, Pris, Bo1) 
= F (go, — 8s1» — f12» Bas, — Bair — Bia, Bsa, — Bys3, — Ba3) 
== F(— 85, 831» — 12, — Bos, Bai, — Bia» — Baz, Biz, — Bo1) 
= F(— go5, — Sar» G12» — Bos, — Bai, Bis, — Bae, — Fis, Boy). (5-3) 
In equations (5.2) the polynomial dependence of /* on gy1, £a2, 833 and t is under- 


stood and in equations (5.3) the polynomial dependence of F on gy, So0, 833, T; 
B,,, Bzz and By, is also understood. 


From (5.2), it follows (cf. SmitH & RIvLin, 1957, § 7) that /* is expressible 
as a polynomial in the quantities 


811» §22> &33> S33, £31 Sie, 03 831812 and Tt. (5.4) 


Also, by a method similar to that used in deriving the result (5.4) and bearing 
in mind that F is linear in B,,, it is seen that the restrictions (5.3) imply that 
F is expressible as a polynomial in the quantities (5.4) together with 


By, Bog, Bgg, 803 Bag, 831 Bay) 812 Pie, 823 Bs2, 831 Piz, 812 Par, 
831 812 Bos, 819 803 Bai, 823 831 Pia, 831 819 Pso, 812803 Pigs, 893831 Bei. (5-5) 


Then, employing the relation (3.13), we see that /,, may be expressed in the form 


11> 812%12+ S23 831P12, &31%31 + 81289331 
el &12%12 + Sos 831 Pia, Hoe, 03 %e3 + £31 812he3||, (5.6) 
£31%31 + 8129 823Ps1» 223 %a3 + 891 812 Pas, X33 


where the «’s and f’s are polynomials in the quantities (5.4). We thus have the 
result that for materials with rhombic symmetry, the stress is given by (2.9) 
in which /* is a polynomial in the quantities (5.4) and f,, is given by (5.6). 
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Divergence Formulas Involving Vorticity 


Louis N. HOWARD 


Communicated by C. C. LIN 


I. Introduction 


There are many “conservation theorems” in hydrodynamics involving the 
vorticity vector. Perhaps the most familiar is the theorem of PoINCARE that the 
“center of gravity” of vorticity in a viscous incompressible 2-dimensional flow 
at rest at infinity remains fixed. Such theorems are of interest because they 
provide information of an overall character without requiring knowledge of specific 
solutions to the flow equations. A thorough discussion of these theorems, with 
derivations and references, is given in C. TRUESDELL’s monograph [7], particularly 
in Chapter VI. Almost all of these results rest on a formula of the type: 


f(x, q, &) = div V(a, q, 2). (1) 


Here q is the velocity field and 2 = V xq is the vorticity. Once such a formula is 
obtained, one gets a conservation theorem by applying the divergence theorem 
and formulating sufficient conditions for the boundary integral to vanish. In 
this paper I shall call a formula of the type (1) a “vorticity divergence formula”’ 
if it is zdentically true for all (sufficiently differentiable) fields q, or, part of the 
time, for all such fields as are also solenoidal. 


The purpose of this paper is to investigate and classify all such vorticity 
divergence formulas, and also the analogous formulas in two dimensions. In 
this way one obtains concise and unified proofs of the formulas on which the 
conservation theorems rest; but the main interest here is in finding all possible 
formulas of type (1), and thus, consequently, the form of all conservation theorems 
which can be obtained by these methods. It turns out that practically all the 
possibilities lead to known conservation theorems, although there are a few 
which seem to have been overlooked. 

While almost all known conservation theorems involving vorticity can be 
traced back to a vorticity divergence formula of type (1), I do not mean to imply 
that these are the only ones of interest. For example there is the formula of 
ByJorcuoM (see [2]) which may be expressed as 


g-Q = div (hQ) 


where h is the “Monge potential’”’ occurring in the representation of q as q= 
ViA+ Veg. This is not of the form (1). 
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Some divergence formulas are commonly written in a tensor form which 
appears to be different from (1), for example 


Q = div(Q~a). 


But any such formula can be reduced to the scalar form (1) by contracting with 
an arbitrary constant tensor, ¢.g. k»-&—=div(Qxk), and conversely from the 
scalar form the tensor form can be recovered by using the arbitrariness of the 
constant tensor. Thus one gets all such tensor formulas as by-products of the 
study of (1). 


II. 3-dimensional Vorticity Formulas for General Flows 
We attempt to find those functions / (a, q, 8&2) which are divergences of vector 


functions of the same arguments identically in q (and 2=Vxq). Now 


DV inate av; 
=—=- Ve a QQ (2) 


f(x, g, 2) = div V(x, q, 2) ery. ‘, 20, Le 


Since this is to be an identity and there are no second derivatives of q on the 
left hand side, one gets as a first necessary condition 


av; 
62; pee (3) 
Now Q; ; is arbitrary except that 2,;;=—0. Thus a must have the form 
lies gy 6;;. Contraction yields gp = 4A OV and $0 : 
02; if oF Bag: 
OVe. tM 
Eee 82; 6, 3. (4) 


av, 
This shows that 2, 


terms are the same; thus 


=0 if a=+-7, while if «=7 one sees that the three diagonal 
0V,(Q2,) __ @V, (22) __ @V5(Q3) __ : 

, | TS NMMEALT Maur Ca A. Since these are 
functions of different Q,’s, they must all be constants in Q, 7.e. functions only 
of # and q. Thus finally 


Vi (a, q, 8) = A (x, q) Q; + B;(ax, q). (5) 


This relation has been proved necessary for (3) to hold and is easily seen to be 
also sufficient. 


(5) must now be put back into (2), yielding 


oA OB; 0A OB; 
Q) = <Q. t aoe 4 

To make the right-hand side a function of a, q and Q it is necessary that 
the quantity in brackets be antisymmetric in 7 and 7, and since this must be 
true identically in x, q, and Q it follows that 


aA OA er, 

2g ok sa , (7a) 
OB Ore 
0g; 1 Og; 7) 
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O ; : 
oe =0, thus A=A(«), while (7b) implies, as is 


v 
familiar in kinematics, that the vector B has the form: 


(7a) clearly implies that 


B= F(x)xq+G(a). 


(“If the deformation tensor is zero, the motion must be a translation and in- 
finitesimal rotation.’”’ A very simple direct proof is the following: Write (7b) 
as B, ,+ B, ,=0 and differentiate to get B; ;, = — B, ;,. One also has, of course, 
Pi Bia, 

. By jn She By Zainal ol B; ni SiS By 53 rj By i; 7 Bi gj = See, By jr: 
..B; ;,=0, showing that B; is a linear function of the g;. Putting a general 
linear function into (7b) one gets the above result almost at once. Cf. [4]). 


Finally, putting these back in (1) one obtains 
f(a, g, 2) = div[ A(x) Q4+ F(a) xq + G (x) | 


=VA-2+ (VXF)-q—(VXq)-F+V-G 
or 


f(a, g, 2) =(VA —F)-Q4 (VXF)-q+V-G. (8) 


This may be simplified by observing that first of all the term V-G is of 
no interest since it does not involve q or ®, and second that VA may be in- 
corporated into F if one is interested only in the form of /. Thus the most general 
form of f is 


f(x, g, 2) = — F(x)-Q + (VXF)-q. (9) 


I 


If a “pure” vorticity divergence formula is defined as one in which the 
second term of /, that which does not involve Q, is absent, it is clear that a pure 
vorticity formula results only if F is irrotational, and in this case one may as 
well write F=—WA. The most general pure vorticity divergence formula can 
thus be expressed in either of the forms: 


VA-Q=div[4 Q] 
VA-Q=div[qxVA]). (10) 


For some purposes one of these equivalent forms may be more useful than the 
other. 


The first form of (10) is closely connected with the vorticity theorems which 
TRUESDELL [3] has obtained from the ‘“‘symmetrical moments of vorticity”’ (cf. 
also BERKER [4] Thm. II to which (10) is essentially equivalent). Assume that 
A(a) is the general homogeneous polynomial of degree m+ 1: 


A(H) = Ae ising iy Mig ++ Minas 


where A; ;,,, is Symmetric in all its indices. Then (10) becomes 


m+ 


(W4A) Ag inines i, ++ %ig 23,4, =A ? 25%, 82; 


In+1 ty we Ie 41 OX: xi, in+1 
7 


8* 
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Since the A;,;,,, are arbitrary except for symmetry, one sees that the sym- 
metric part of 


{(n +1) Mies Hy (oe = a4 [eee es H;,. 25] 


in 
or 


Ais aialye eg te Kp wn Nigga hoy, Vimy chat ei eens rte ana ee = [Az ca Bera eeule 

These are exactly the formulas leading to the vorticity theorems connected 
with the symmetrical moments of vorticity given by TRUESDELL (‘‘First General 
Conservation Theorem”). In the formal sense one may regard “‘any’’ function 
A(x) (expanded in a power series) as a linear combination of homogeneous poly- 
nomials and so regard any pure vorticity divergence formula as a “‘linear combi- 
nation” of TRUESDELL’s formulas. In this sense, the vorticity theorems connected 
with the symmetrical moments of vorticity are the only pure vorticity theorems 
(in three dimensions) valid for all flows. Of course, the “impure” vorticity 
divergence formulas given by equation (9) may also be expressed in terms of 
a polynomial basis, if desired. 


III. Vorticity Divergence Formulas for Solenoidal Flows 


If the identity (2) is required to hold only for all solenordal fields q (rather 
than all q), one can expect to find more vorticity theorems. I remark first that 
the condition V-q=0 imposes no (algebraic) restriction on 22=V xq, and that 
consequently the first part of the argument above, leading to equation (5), goes 
through as before. A difference appears, however, at the next step, since the 
q;,; are no longer arbitrary but are restricted by the condition g; ;—0. Thus one 
can no longer conclude that the symmetrical part of the quantity in brackets 
in (6) is zero but now only that it is a scalar multiple of 6;;. Separating the part 
involving & from that independent of Q yields 


0A 0A 
and 
aR, , ane 


0A 2) : : : 

where yp = = Bar Q, and y= > a , as is easily seen by contracting (11) and (12). 
k 

As before, (11) implies oa =o hence A = A(x). For example, by differentia- 


ting (11) with respect to Q, follows 


oA aA). <hdip. epee 
0g; ee 0G;  GQN Brgy 
0A 
So 
0”; 


(12) is more interesting, for one now really gets something different. For the 
moment, denote partial derivatives with respect to g; by subscripts preceded by 
Mays 


7 


a comma: B,;. Then, differentiating (12), one gets 


By int By th = Qin 6;;- (13) 
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Now solve these for B; ;, in the manner familiar in connection with the Chri- 
stoffel symbols. By interchanging 7 and k in (43) follows 


Q By ji t+ By iz = Pi On- (14) 
Subtracting (14) from (13) yields 
Bi, 75 = Orn 0,5 5 Ope (15) 
Interchange 7 and k in (15) and add the result to (13): 
2 Bi pk = Pin 9:3 + Pj Orn — Pri On: (16) 
One can now show that the ,; are independent of q. Differentiate (16): 
2B, pnt = Poni 9:5 + P51 Ok — Prit Ojn- (17) 


Antisymmetrize this in k, J, obtaining 
O= 51 Os Ps ik OT "0,41 Ojn as ee Oj 


Contract this by setting 7 =k and summing: 


0 = 3,51 — Piji—- Pig + Pris HF 
or 
P51 = — Pris Oj1- (18) 
Contracting (18) yields ,;;= —39,;;, hence y,;;=0, and then from (18) 9, ;;=0. 
Thus the »,; are independent of qg. Set y,;=2K;=2K;(a), and return to equa- 
tion (16), which now becomes 

B 
Integrating this leads to 


By, = K,%, 0 


Ky 03; ME ee. 0; % ee K; Orn 


17k — 


K; 9; KG 4 L 


yom 49? 


where the L;, are independent of q. Inserting this in (12) shows that L,;+L;;= 


2a 6;;, where « =a (a) is a scalar independent of q. Thus L;; must have the form 
Lj; = « 0;; + By (&) &x5;- 
Therefore 
B; ;(@, Q) = Ky (@) 9, 0;; + Kj 9; — Ki q + &(x) 0;; + Be (®) ei; 
= (Ki, % + &) 0;; paalS PC Se esl ER iq: (19) 
Integrating (19), one obtains 
B, (x, Q) = Ky % 4 — & Ki GG + %(@) 9; + Ba (&) &n 5 9; (20) 


(dropping the “‘constant’’ of integration which is independent of q). 
Finally, putting (19) and A = A(x) into (6) yields 


} (a, q, &) = Fe fe ae t+ [Exim G1 Kin + By] 2, 


= oAie) a5 Exlm HK, py Q,+ 


OX, 


' OK, <l OK; Oe (led OBp A 
b Qi Ox, Th > ax; qj an; Git ax; Rij Vj 
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or 


f=(VA+qxK+B]-2+¢4-VK-q—2V-Kq-¢+Va-q—G-VXB. (21) 


This is the most general vorticity divergence formula valid for arbitrary 
solenoidal flows. 
If one looks again for “pure’’ vorticity formulas, one must have 
OK Uh bak) eee OK pe eae 5 
ee sie se} qi qj Oxp 0; qi qj =— ( ) 
and 


(Va —VxXB)-q=0. (23) 


From (23) follows VX B=Va, hence Vx (Vx B) =0 and 4a=0. From (22) 


one gets 
ok; 


aK; OK, 
Ox; OX, 


0; 
mee KG 
and, contracting, obtains ae 0; hence 
k 


CE Oe 


OX; Ox; Os 


Therefore 
K(x) = —[kxx+1]], 


where k and lI are constant vectors. Putting these results into (21) and incor- 
porating VA into B, we find for the most general pure vorticity divergence 
formula valid for solenoidal flows: 


f(a, q, &) = [(kxa +) xq + B]-R (24) 
where VX (VXB)=0. This may be rewritten as 
f(x, q, 2) =k-xx (qxQ) + 1-qxQ+ B-Q (25) 
where k and U are arbitrary constant vectors and B = B(#) satisfies the equation 
Vx (Vx 8) =0. 


It is perhaps more interesting, however, to consider (25) modulo the formulas 
which are true for arbitrary flows. It is clear that then there emerge only the 
three independent formulas 


k- (a@x (q¢xQ)) = div[k@kxa — (kxax)-qq] (26) 
(which leads to Poincarsk’s integral), 
I-qxQ = div [lt 7 —l.qq| (27) 


(which leads to Lamp’s integral), 


and £§-Q2=div[«q+ qx], (28) 


where VX ®=V«a, modulo such formulas, with ® irrotational, as are known 
to hold for arbitrary flows. Since B is determined up to a gradient by «, one 
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may expect to denominate such formulas (28), modulo those that hold for general 
flows, by the scalar function «, which is an arbitrary harmonic function. One 
way of doing this explicitly is as follows: Given a harmonic function «(a), de- 
fine B by 


B (x) Loiviteaxeva® 


This solves VX ® = Va; when supplemented by an arbitrary gradient, it gives 
all solutions. Two different «’s give linearly independent f’s modulo gradients 
unless %,=aa,+b, so that except for this rather trivial transformation one 
gets a different formula (28) for each harmonic function « (modulo the formulas 
that hold for general flows). In particular, looking for a basis of polynomial 
vorticity formulas, one sees that besides the special formulas (26) and (27) and 
the generally true symmetrical moment formulas of TRUESDELL, in the solenoidal 
case there exist exactly as many more linearly independent homogeneous for- 
mulas of degree ” as there are homogeneous harmonic polynomials of degree n, 
1.€. 2n+1. (n=1, since for »=0 follows «=const. and B=0. The additive 
arbitrary constant b above destroys the polynomial of degree 0, but of course 
has no effect for ” >0.) 

For example, in the case 7 =1, besides the six independent formulas given 
by TRUESDELL in the dyadic form 


rQ+Qr = div[rrQ] 


there are three more independent formulas for solenoidal q obtained by choosing 
a=a-r. By use of the arbitrariness of a these may be combined into the vector 
formula 


rxQ = div [2rq—qr-+ (q-r) 1]. (29) 


IV. Two Dimensional Formulas 


I look first for formulas valid for arbitrary fields gq. Writing out equation 
G2 ON 


OX, OX» 
immediately that V must actually be independent of ¢, and so obtains 


(1) for the two dimensional case with € = replacing &2 one sees almost 


OV, “eV, 
f(a, q; ¢) = div V(x, q) a Ox; ie aq; Gj, is (30) 
re ee! eee 
For this to depend only on x, q, and ¢ it is necessary that oan =f eqn 0. This 


equation implies as before that V, must be linear in q, and from the antisym- 
metry follows 
Vi (a, @) = 9 (@) €:5 9) + Yi (@). (31) 


Dropping y; as leading to nothing of interest, one then obtains the most general 
possible formula, with 


f= 22) 2,54 + p(@)E. (32) 


* This expression for the vector potential is perhaps the simplest. For further 
discussion of it see Amer. Math. Monthly 58, 487—489 (1951). 


— 
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=, so that there is 


ae a 
For ‘‘pure”’ vorticity formulas it is necessary that ae 
t 


essentially only the familiar one: 
C = div (¢2, — 4). (33) 


Turning now to formulas valid for arbitrary solenoidal fields q, we find more 
interesting results. As before we start with (30), but now can conclude only 


a, OV, ag Leh 
These equations are equivalent to the following two: 
eV Ep 
04s Oy 
ONT 
0q, O92" (35) 


Since these are the Cauchy-Riemann equations, V, and V, must be conjugate 
harmonic functions in g, and q, (for each point (x, x,)). Thus follows 


which is the most general vorticity divergence formula in two dimensions for 
solenoidal fields. 


To discuss the case of pure vorticity formulas it is convenient to use complex 
notation. Let W¥+71V, =v, 9, +7 4.=4, +1 *,=z. Then the Cauchy-Riemann 
equations (35) take the form 

ov 
a (37) 
(Here 9,9, 2,2 are regarded as independent variables.) For a pure vorticity 


theorem it is now necessary also that a =0, which in the complex notation 
becomes 


Re \ae} = (): (38) 


Now dv/éz (assuming continuity of the mixed derivatives) considered as a 
function of the complex variable g is analytic, and since by (38) its real part 
vanishes, it must be a (pure imaginary) constant in q, 7.e. a function only of z 
and Z, or a. Integrating yields 


v=h(z,2) + g(q, 2). (39) 


The h term may be dropped since it gives a term in V which is solenoidal by 
itself and leads to nothing of interest in f(x, q,C), so finally one gets for the 
most general pure vorticity divergence formula for two dimensional solenoidal 
flows 

v = v(q, 2) 
and 


-OU), 4. fu+tv v—d 
Re {i ate div ( eae ). (40) 
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The simplest case of (40) is that with v=w/(z) g, which gives f(x, q,¢)= 
Re {7 w (Z)} € = (a) € where 9 is an arbitrary harmonic function. This is a result 
pointed out by HaAmeEL. It is interesting that there are in fact many more such 
formulas analogous to HAMEL’s. The other extreme case is of course f= @(q)¢ 
where here w can be an arbitrary harmonic function of gq. A polynomial basis 
for these formulas is given by the real and imaginary parts of 

fom = (ty — 4 5)" (Gy + 4 qa)"C = div ata)" a tae)” 


m+1 


i,—1) (41) 


forn+m=1. For n=m=0 one gets simply the one independent case /=¢, 
which is in fact a divergence whether or not div g=0. 


V. Remarks on Time Derivatives 


The formula ¢=div(q¢,, —q,) shows that if g—0 sufficiently rapidly at oo 
then the total vorticity [¢dA must be zero. ‘‘Sufficiently rapidly” in this case, 
however, means faster than 1/7, and this is perhaps rarely the case. Neverthe- 
less, as PoINCARE showed, if for an incompressible viscous fluid one assumes 
only that €-+0 sufficiently rapidly at oo, one can establish the constancy in 
time of the total vorticity. This can be shown as follows. For such a fluid the 


Navier-Stokes equation leads immediately to vee vAC or 


se 
S+V-(G0) =r V-VE 
or ac 
ieee Eee! (42) 


Let 2 be a fixed area with boundary S. Then 


a A tea A : 
x ft44= [ae4= fee ql) -ndl. 
QI 5) S 


On letting %& expand to the whole plane and using the assumption that ¢ 
and its derivatives tend rapidly to zero at oo one arrives at the constancy of 
the total vorticity. 

Now, keeping this argument in mind, try to generalize it to expressions of 


the form 
af f@caa = [ tw) 4a. 
Df ») 


By use of (42) this becomes 
ac 4 NS. 
[ras frv-tve—aciaa 
w w 


= [V-oVe—gs) —v6Vf\da + f p(Ne+a-Vitldd. 
a oT 
Thus to express this as a line integral, identically in q, one must express 
q-V/¢ asa divergence. But it has been seen that this can be done only ifq-V/ 
is the real part of a function analytic in g,+79, and x,—1x,. Call this function 
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u+iv, so that qf, +%,f,,=4. Forming the harmonic conjugate in q yields 


v=—hr +dte +C (a), and then the Cauchy-Riemann equations in x, and x, 
give j 
VN Pees ne G3 Tse: = Ot binoh ; LAP ey = Cy 
Giles alts Gator Sh ipete = Vatican: a Cs. 
from- which: C =.const.j) 9. = fy, ta Oh 
f=a(xtt x4) +b-x+k. (43) 


This is the only possible form for f, and in fact the three independent parts of 
(43) lead exactly to the three familiar formulas of PoINcaRE: 


é is 

5 [64a 25 

foal eh GAO 

ot 

= fw-wtdA =4y [CaA. (44) 


For these one must assume q—>0 at co and €-0 sufficiently rapidly. It is not 
necessary here to go into more detail; the important point at present is that 
(43) shows these formulas of PorncarE to be the only such formulas so obtainable. 

In the three dimensional case, TRUESDELL has shown that the time deriva- 
tives of the symmetrical moments of vorticity can be expressed as surface inte- 
grals in such a way that certain order conditions at infinity weaker than those 
required to give the vanishing of a given symmetrical moment will nevertheless 
imply its constancy. Besides these formulas, however, there are also, in thecase 
of solenoidal q, two formulas of MorEAu (see [/] for references) which may be 
written : 


= [ex CNS: ih adV + boundary integral 
B 


B 
a — (x -x)QdV =2 [wxaav + boundary integral. (45) 
B B 


where @ is the acceleration or the acceleration plus an arbitrary gradient and B 
is a fixed volume. Under suitable conditions (possibly difficult to satisty) one 
may expect the surface integrals to vanish as 8 expands to the whole space. 
I now show that the formulas of TRUESDELL and, in the solenoidal case, those 
of Moreau are the only possible ones of this type. As before, we start with 


ay | f@) -Qav = fiw. dV. 
B B 


But from 


a) 
a=“T+9.Vq=“14+Qxq +V(>-7) 
follows 
OQ 


ap 7 VX@+ VX (GxQ). (46) 
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(At this point one may modify a if desired by adding a gradient.) Then 
a [ fav = [7-0 a+ 9x2] dV 
B B 


a | -V-[fx(a+qxQ)] dV + [ Vxf-la + qxQav. 
B B 


Formulas leading to TRUESDELL’S are obtained by choosing f to be the gradient 
of a homogeneous polynomical, and it has been seen in the general case that 
Vxf+qxQ is not (identically) a divergence unless VXf is zero. However, in 
the solenoidal case (25) shows that Vxf-qx& will also be a divergence only 
in case 


Vxf=kxx+l. 


(The term involving a must be left as it is since @ is unknown until one makes 
some dynamical statement about the flow.) A suitable vector potential for this 
is f= —$(x-x)k-+4UX«, and using the arbitrariness of k and I one gets exactly 
the formulas (45). Thus TRUESDELL’s and MorEAv’s formulas are the only ones 
of this type. 


VI. Index of Formulas 
1. Three Dimensions 


a) General Flows: (9) 

b) General Flows, pure case: (10) 

c) Solenoidal Flows: (21) 

d) Solenoidal Flows, pure case: (24), (25) 

e) Same, modulo b): (26), (27), (28) and following discussion. 


2. Two Dimensions 


a) General Flows: (32) 

b) General Flows, pure case: (33) 

c) Solenoidal Flows: (36) 

d) Solenoidal Flows, pure case: (40) 
e) Same, polynomial case: (41). 
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Nachweis periodischer Losungen bei gewissen 
nichtlinearen Schwingungsdifferentialgleichungen’ 


HANS EHRMANN 


Vorgelegt von L. COLLATZ 


§ 1. Einleitung 


Unter den bei eingliedrigen Schwingern haufig vorkommenden und in bezug 
auf die praktische Bestimmung von periodischen Lésungen viel untersuchten 
Schwingungsdifferentialgleichungen gibt es zahlreiche, bei denen die Existenz 
von periodischen Lésungen noch nicht nachgewiesen wurde. So ist z. B. bei der 
Differentialgleichung ¥-+« «+f x8 =c cosw t, B>0, m. W. bisher nur fiir kleine 
B die Existenz einer periodischen Lésung mit der Periode 2/m gesichert. Ferner 
ist bisher, soweit mir bekannt ist, nur in sehr speziellen Fallen bei nichtlinearen 
Schwingungsdifferentialgleichungen die Existenz von mehreren periodischen 
Lésungen derselben Periode nachgewiesen worden. Ein bekanntes Beispiel hier- 
fiir ist die sog. Duffingsche Differentialgleichung ¥-+ksin x=csin?#, die von 
zahlreichen Mathematikern untersucht wurde (DuFFING, HAMEL, HAMMERSTEIN, 
IGLISCH u. a.). SchlieBlich gibt es, soweit mir bekannt ist, wenn man von speziellen 
Schwingungsgleichungen absieht, noch keine allgemeinen Existenzsatze fiir sog. 
subharmonische Lésungen, d.h. solche periodische Lésungen, deren fleinste 
Periode ein Vielfaches der Periode der erzwingenden Kraft ist, obwohl die Exi- 
stenz solcher Lésungen auch fiir die Anwendungen von Interesse ist und sie 
daher in bezug auf die praktische Bestimmung viel untersucht wurden. 

In dieser Arbeit soll nun gezeigt werden, daB es eine groBe Klasse von nicht- 
linearen Schwingungsdifferentialgleichungen zweiter Ordnung gibt, die sowohl 
mehrere periodische Lésungen der Periode der erzwingenden Kraft als auch sub- 
harmonische Lésungen besitzen. Es zeigt sich dabei, daB zu dieser Klasse im 
wesentlichen solche Differentialgleichungen gehéren, bei denen auch noch kein 
Existenzsatz fiir mindestens eine periodische Lésung vorlag. Dabei handelt 
es sich hier zunachst nur um die Anwendung eines allgemeinen Existenzsatzes 
fiir die Lésungen von Randwertaufgaben, der in einer fritheren Arbeit aufgestellt 
wurde, auf den Nachweis periodischer Lésungen. 

Zu diesem Zweck werden in § 2 zunachst einige elementare Untersuchungen 
durchgefithrt, um die Bedingungen aufzustellen, unter denen aus der Existenz von 


S Einige der hier vorgetragenen Ergebnisse sind bereits in meiner Dissertation, 
Existenzsatze fiir periodische Loésungen bei eingliedrigen nichtlinearen Schwingern, 
enthalten; Hamburg, Juli 1954, Hauptreferent Professor Dr. L. Corratrz, dem ich 
auch hier fiir einige wertvolle Ratschlage danke. 
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Lésungen gewisser Randwertaufgaben, insbesondere der 1. Randwertaufgabe 
X(t) =%, %(4)=%,, bzw. der 2. Randwertaufgabe x(t.) =x), %(t4) =, die 
Periodizitat der Lésung folgt. Es handelt sich dabei um einfache Verallgemeine- 
rungen von Uberlegungen, wie sie bereits F. LEttENMEYER! bei der Duffing- 
schen Differentialgleichung angestellt hat. Die Ergebnisse sind in Satz 8 zu- 
sammengefa8t, der es erméglicht, unmittelbar gewisse Existenzsatze fiir Rand- 
wertaufgaben auf den Nachweis periodischer Lésungen anzuwenden. 

In § 3 wird der Existenzsatz fiir Randwertaufgaben noch einmal ausfiihrlich 
formuliert, wodurch diese Arbeit auch unabhangig von der vorhergehenden lesbar 
ist. Darauf werden einige Existenzsatze fiir periodische Lésungen aufgestellt, 
wobei zunachst von der allgemeinen Differentialgleichung 2. Ordnung % = f(t, x, %) 
ausgegangen und dann erst die Spezialisierung auf nichtlineare Schwingungs- 
gleichungen durchgefiihrt wird. Dabei werden auch rheonichtlineare Schwinger 
der allgemeinen Form 


H+ a(t, x, 2% +e x) =Pl0 


mit erfaBt. Es zeigt sich, daB unter gewissen allgemeinen Bedingungen fiir 
die Funktionen d, g und # unendlich viele verschiedene periodische Lésungen 
existieren. 

In § 4 gelingt es dann, fiir dieselbe Klasse von Schwingungsdifferentialglei- 
chungen auch die Existenz von subharmonischen Lésungen nachzuweisen. Dabei 
werden wiederum gewisse Uberlegungen der fritheren Arbeit verwendet, die auch 
hier in Form von Satzen zunachst noch einmal formuliert werden. Die Unter- 
suchungen iiber subharmonische Lésungen in § 4 stellen das Hauptergebnis dieser 
Arbeit dar. 


§ 2. Bedingungen fiir die Periodizitat einer Losung 


Wir treffen zunachst folgende allgemeine V oraussetzung: 
A: Fiir alle in dieser Arbeit auftretenden Differentialgleichungen 


(1) & = f(t, x,%) 
gelte: 


1. Die Funktion f(t, x,%) set fiir alle (endlichen) Punkte des t, x, %-Raumes 
eindeutig definiert. 
2. Sie sei so beschaffen, daB die Gleichung 


(2) % (t) —% (t)) cadet, x(t), %(t)) dé 


fiir alle Anfangswerte ty, x (to) = %,% (to) =%o genau eine fiir alle t Sty stetig diffe- 
venzierbare Losung x(t) hat, die die Anfangsbedingungen ty, X%), Xo erfullt und sich 
stetig mit diesen Anfangswerten dndert. 

3. Es existiert eine positive Konstante T derart, dag die Gleichung 


(3) (t+ T, x, %) =f (t, x, %) 


fiir alle t, x und x tdentisch erfiillt rst. 


1 LETTENMEYER, F.: Uber die von einem Punkt ausgehenden Integralkurven einer 
Differentialgleichung 2. Ordnung. Deutsche Math. 7, 56—74 (1942). 
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Die Bedingungen 1. und 2. sind z. B. erfiillt, wenn die Funktion /(é, x, v) 
stetig ist beziiglich ihrer 3 Argumente und entweder stetige partielle Ableitungen 
nach x und v hat oder in jedem endlichen Gebiet des ¢, x, v-Raumes einer Lip- 
schitzbedingung der Form 


(4) \¢(é, «, ») —f(,%, | <M,|*% —%| + M|v —9| 


gentgt. 

Jedoch ist die Stetigkeit von f(¢, x, v) beztiglich ¢ nicht notwendig, damit 
die Voraussetzung A gilt, wenn (4) in jedem endlichen Bereich des ¢, x, v-Raumes 
erfiillt ist und noch gewisse einfache Zusatzbedingungen gelten®. 

Ist x(t) eine Lésung von (2), so erfiillt x(¢) tiberall dort, wo der Integrand 
stetig ist, die Differentialgleichung (1). Die Differentialgleichung (1) ist dann 
bis auf eine Punktmenge vom MaBe Null erfiillt oder, wenn wir uns im Hinblick 
auf die Anwendungen auf stiickweise stetige Funktionen f beschranken, gilt (1) 
dann in jedem endlichen Bereich bis auf endlich viele Unstetigkeiten. 


Aus den Voraussetzungen folgt sofort, daB mit x=g(t) auch x=p(t+kT) 
mit ganzem & Lésung von (1) ist: 


(5) GLAT) =f(t LAT, p(t LAT), PLAT) =F (4 PEHAT), PE LAT). 


Hieraus ergibt sich weiter der elementare 


Satz 1. Notwendig und hinreichend dafiir, daB eine Lésung von (1) bzw. (2) 
periodisch mit der Periode kT, k =>1, ganz, ist, ist das Erfiilltsein der Randbe- 
dingungen 


(6) (i qObT Sxlile wa tak ees 


Beweis: Die Notwendigkeit ist unmittelbar klar. DaB die Bedingungen hin- 
reichend sind, folgt aus der Eindeutigkeit der Lésung der Anfangswertaufgabe 
und aus Gleichung (5): die Lésungen x= g(t) und x=@(t)=q(t+T) fallen 
zusammen. 


Zur Anwendung vieler Existenzsatze fiir die Lésung von Randwertaufgaben 
auf den Nachweis periodischer Lésungen von (1) bzw. (2) sind jedoch haufig die 
Randbedingungen (6) nicht unmittelbar geeignet, sondern miissen durch solche 
ersetzt werden, von denen sich jede nur auf eine Randstelle ¢) bzw. #, bezieht. 
Zu diesem Zweck sind noch gewisse zusdtzliche Symmetrieeigenschaften der 
Funktion /(¢, x, %) erforderlich, die wir jetzt herleiten wollen: 


Satz 2. Die Differentialgleichung (1) sei so beschaffen, daB mit x =q(é) auch 
x= g(—?) Lésung ist. 

Dann ist cine Lésung mit den Randbedingungen % (0) =%(2kT) =0, (R= 1, 
ganz), periodisch mit der Periode kT und symmetrisch zu den Geraden t=0 
und ¢=$kT, und eine Lésung mit den Randbedingungen % ($ T) = « (kT) =0, 


k = 2, ganz, ist periodisch mit der Periode (k —1)T und symmetrisch zu den 
Geraden ¢=3T und ¢=2AT. 


2 Vegi. z. B. Kamxg, E.: Differentialgleichungen, Losungsmethoden und Lésungen, 
Bd. 1, Gewéhnl. Dgln., 3. Aufl., S. 33f. Leipzig 1944. 
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Beweis: a) Nach (5) ist mit p(t) auch p(t{—kT), k=1, ganz, Losung von 
(1) und nach Voraussetzung ist daher mit (¢) auch p(t) =q(kT —1#) Lésung. 
Man erhalt daher wieder eine Lésung von (1), wenn man q/(é) an der Geraden 
t= 3hkT spiegelt. Eine Losung mit der Randbedingung % (47) =0 fallt wegen 
der Eindeutigkeit der Lésung der Anfangswertaufgabe mit ihrem Spiegelbild 
beziiglich der Geraden t= $kT zusammen: p(t) =@(kT —2). Daher gilt fiir sie 
x(0)=x(kT), x(0) =—2%(RT). Wegen %(0)=0 sind also die beiden Randbe- 
dingungen (6) erfiillt, d. h. die Lésung ist periodisch mit der Periode RT. Wegen 
x%(0)=0 fallt sie auBerdem mit ihrem Spiegelbild beziiglich der Geraden t=0 
zusammen und ist daher auch symmetrisch zu t= 0. 


b) Entsprechend fallt eine Lésung mit den Randbedingungen « (} T) = % (ERT) 
=0O mit ihren Spiegelbildern beziiglich der Geraden t=47 und t=$hkT zu- 
sammen, ist also symmetrisch zu diesen Geraden. Daher gilt x($7) = x(kT— 4 ; 
E (iT 4 (k—1)T) und #G7)=—*%($7+(k—1)T) und wegen %($T)= 
also auch %*($7)=%(37T+(k—1)T). Nach (6) hat die Lésung daher die 
Periode (k —1) T. 


Satz 3. Die Differentialgleichung (1) sei so beschaffen, daB mit x= q(t) auch 
x— —o(—t) Losungiist. 

Dann ist eine Lésung mit den Randbedingungen x(0) = x($kT) =0, R=1, 
ganz, periodisch mit der Periode kT und punktsymmetrisch zu den Punkten 
t=0, x=0 und t= ZkT, x=0, und eine Lésung mit den Randbedingungen 
4 (37)=x($kT)=0, R=2, ganz, ist periodisch mit der Periode (k —1)T und 
punktsymmetrisch zu den Punkten ¢=37, x=0 und t=$kT, x=0. 


Beweis: a) Aus der Voraussetzung und Gleichung (5) folgt, daB mit « = »(?) 
auch w(t) = —@(kT —2) Lésung ist. Wegen »($kT)=0 und der Eindeutigkeit 
der Anfangswertaufgabe fallen beide Lésungen zu einer Lésung zusammen, die 
spiegelbildlich zum Punkt ($AT, 0) verlauft: w(4)=—qg(kT—?#). Daher gilt 
~(0)= (RT) und v(0) =—g(kT) =0= (kT). Die Lésung ist also periodisch 
mit der Periode kT. Wegen (0) =0 fallt sie auch im Punkte (0, 0) mit ihrem 
Spiegelbild beziiglich dieses Punktes zusammen. 


b) Wie in a) gilt p@) =— (kT —2), p(t) =p(kT —2), wobei der Punkt die 
Ableitung nach dem Argument bedeutet. Daher gilt (3 T) =9 (2 T+ (k —1) T) 
und wegen 9 (2 i) =O p (2 = (2 T + (k —1) zy, also Periodizitat mit der 
Periode (k—1)7. Nach Voraussetzung und Gleichung (5) folgt ftir sie auch 
y(T —t) =— y(t), also Symmetrie zum Punkt ($7, 0). 


Satz 4. Die Differentialgleichung (1) sei so beschaffen, daB mit «= q(t) 
4. auch — y(—2), 2. auch y(zT —1#) Lésung ist. 

Dann ist eine Lésung mit den Randbedingungen x(0)=%(;k7T)=O0 mit 
k=2n+1, n=0, ganz, periodisch mit der Periode kT, punktsymmetrisch zu 
den Punkten (0,0) und ($kT,0) und symmetrisch zu den Geraden (= 4kT und 
t—2kT. Und eine Lésung mit den Randbedingungen x (37) =% (7+ ¢kT) 
=(Qmit k=2n+1,n=0, ganz, ist ebenfalls periodisch mit der Periode kT, punkt- 
symmetrisch zu ($7, 0) und [Ss i 0); sowie symmetrisch zu t=$7+ kT 
und ¢=347-+ 2kT. 
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Beweis: a) Nach Voraussetzung und (5) ist mit @ (2) auch 9(¢7+nT—?t) = 
p (kT —1) (n=0,1,...) Lésung. Beide Losungen liegen spiegelbildlich zur Ge- 
raden t=1kT und fallen wegen % (LRT) =0 und der Eindeutigkeit der Anfangs- 
wertaufgabe zusammen: 9 (t) = p($kT— ?). Daher ist mit y(0) =0 auch p(s kT) =0. 
Aus Satz 3 folgt die Periodizitat mit der Periode kT und die Punktsymmetrie 
zu (0,0) und (4&7, 0). Hieraus und aus der Symmetrie beziiglich t=$hkT folgt 
p (ERT +t) =—p (ERT —) = — et) =—v(kT +1) = (RT —9), also Symmetrie 
beziiglich der Geraden t= $kT. 


2m+1 7 
2 

Mit k=2n,—1 erhalten wir aus (¢) die Lésung p(T + 3T —12). Beide legen 
symmetrisch zueinander beziiglich der Geraden t=} 7-+ 4hT und fallen wegen 
m(&XT+4ERT)=0 zusammen: p($T+4kT—-)N)=o(ET+GghT+. Wegen 
gy (4 T) =0 ist daher auch y(} 7 +37) =0. Nach Satz 3 folgt die Periodizitat 
mit der Periode kT und die Punktsymmetrie zu den Punkten (27,0) und 
(*e4 iG 0). Aus der Punktsymmetrie zu (4 i a 0) und der Symmetrie zur 
Geraden T=47T+H1hkT folgt die Symmetrie zur Geraden ¢=$7 + 2hT. 


b) Wie in a) ist mit g(¢) auch o( —t) Lésung fiir jedes ganze ,. 


Satz 5. Notwendig und hinreichend dafiir, daB mit einer beliebigen Lésung 
x=q(t) auch x=q(—#) Lésung von (1) ist, ist das Bestehen der Gleichung 


(7) f(—t, x, —v) =f, x, ») 
fiir alle Werte ¢, x, und v. 
Bewets: a) Mit p(t) sei auch y(t) =q(—2) Lésung von (1). Dann gelten die 
Gleichungen 
GO =1% 90,¢0) und ¢(—)=f(t, 9(-4, —¢(—-9) 


identisch in ¢, wobei der Punkt jeweils die Ableitung nach dem Argument be- 
deutet. Ersetzt man in der letzten Gleichung ¢ durch —z?, so folgt aus beiden 
die Identitat in ¢ 


H(—4, 9), —9@) =f 90,9 (0). 
Da ¢ eine beliebige Lésung ist, fiir die die Anfangswerte tf), @ (to), P (to) noch 


beliebig vorgeschrieben werden kénnen, muB8 die letzte Gleichung auch identisch 
in jedem der 3 Argumente ¢, p und®@ sein. Die Bedingung (7) ist also notwendig. 


b) Es gelte (7) und es sei y(#) eine Lésung von (1). Dann folgt aus (1) und 
(7) mit »(/) = ~(—2) 


v4) =G9(—1) =1(—-4, o(—4,¢(—9) =f (6 e(— 9, —¢(— 9) =F (69, 0). 


Also ist p(—?) Lésung von (1), und (7) ist somit auch hinreichend. 


Satz 6. Notwendig und hinreichend dafiir, daB mit einer beliebigen Lésung 
x=p(t) auch x= —(—d2) Lésung von (1) ist, ist das Bestehen der Gleichung 


(8) f(t, =). 0) = Ae v) 


fiir alle Werte ¢, x und v, 
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Beweis: a) Mit p(t) sei auch y(t) = —gy(—?) Lésung von (1). Dann gelten 
die Gleichungen 


$) =f eO,¢) und —G(—1) =f (4 ——(—4,¢(—9) 


identisch in ¢, wobei der Punkt jeweils die Ableitung nach dem Argument be- 
deutet. Ersetzt man in der rechten Gleichung ¢ durch —#, so folgt aus beiden 


die Identitat in ¢ ee oe 9. @(0) Pan. i (t (0), 000). 


Da ¢ eine beliebige Loésung ist, fiir die die Anfangswerte t), y(t ), q(t) noch 
beliebig vorgeschrieben werden kénnen, muB die letzte Gleichung auch identisch 
in jedem der 3 Argumente ¢, y und @ sein. Die Bedingung (8) ist also notwendig. 

b) Es gelte (8), und es sei g() eine Lésung von (1). Dann folgt aus (1) und 
(8) mit p() =—9¢(—d: 


oi) = —¢(—) = — 1-4, 9(— 9. 9(—9) 
=—f(-4,—9, 9) =1(6 vO, v@)- 


Also ist auch w(¢) = —q(—#) Loésung. Die Bedingung (8) ist somit auch hin- 
reichend, 


Satz 7. Notwendig und hinreichend dafiir, daB mit x =  (¢) 1. auch x= — y(—2) 


und 2. auch x =g(} 7 —12) Loésung von (1) ist, ist das Bestehen der Gleichungen 
(8) und 


(9) f(¥T —t, x, —v) =f it, x, v) 
fiir alle ¢, x und v. 


Bewets: DaB (8) notwendig und hinreichend fiir den ersten Teil der Behauptung 
ist, folgt aus Satz 6. Es ist also nur noch zu zeigen, dal fiir (9) das Entsprechende 
gilt. Zu dem Zweck fiihren wir in der Differentialgleichung (1) und in den L6- 
sungen die Transformation t=1t+47Z der unabhangigen Variablen ¢ durch. 
Dann wird o(f)=9(t +32) =y(z) und oT —4) =—g(ET—1)=y(— 20), und 
(9) geht mit f(t, x, v) =/(¢7T +t, x, v) =g(t, x, v) tiber in 
(Qa) g(—t, *, —v) = g(t, x, v). 

Nach Satz 5 ist aber diese Bedingung (9a) notwendig und hinreichend dafiir, 
daB mit einer Lésung p(t) der Gleichung y=g(t, y, y), p= ie , auch wy(—7) 
Lésung ist. Durch Riicktransformation folgt damit die Behauptung. 

Wir fassen jetzt die einzelnen Satze 1—7 zu folgendem Ergebnis zusammen: 

Satz 8. Fiir die Differentialgleichung 
(1) % = f(t, x, %) 


set die Anfangswertaufgabe ty, X (to) = Xo, X (to) =H fiir alle ty, Xo, %o eindeutig 
durch eine fiir alle t>t, stetige und stetig differenzierbare Funktion x(t) ldsbar. 
Ferner gelte 


(3) (t+ T, x, %) =f, x, %) 


fiir alle t, x und % mit einer Konstanten T >0. 
Arch, Rational Mech. Anal., Vol. 1 
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Dann ist hinreichend dafiir, daB eine Lésung x(t) die Periode kT mit natiir- 
licher Zahl k hat 


Tabelle 1 
er sn tn aan ean 
Im Falle folgender zusdtzlicher V oraussetzungen Das Erfiillisein der Randbedingungen 
fiir f(t, x, Vv) 
Sr (esr 
1 Mate, =) ={16 70) # (0) =4#(zkT) =0 
in , . = . -(/R+1 
2 ta, 0) =i 7,0) #(ET) =#(“2*T)=0 
3 f(t, =45 0) = — 16 4,2) (0) =%(2kT) =0 
| 1 k+1 
4 f(—t, —4#,v) = —f(t #, v) x (ET) =2(=* T)=0 
| [ooh ee iene ae x (0) =%#(tkT) =O, k ungerade 
5 | #(RT —t, x, —v) =f (t, ¥, v) x ( ) “(4 ) , to} 
| f(-4—")=—fe%), | a7) —2GT+2RT)=0, h ungerade 
ON ES taxon ee a) ala Thea OG Teele 


Es sei noch bemerkt, daB diese einfachen Untersuchungen nichts dariiber 
aussagen, ob die sich ergebenden Lésungen mit den speziellen Randbedingungen, 
falls sie existieren, jeweils die einzigen periodischen Lésungen mit der Periode 
kT sind und ob kT die kleinste Periode ist, was im allgemeinen nicht der Fall 
sein wird. 

§ 3. Existenzsatze flr periodische Lésungen 


Wir wollen die Untersuchungen des letzten Paragraphen nun auf folgenden 
Existenzsatz fiir die Losungen von Randwertaufgaben anwenden?: 


Existenzsatz: Die Differentialgleichung (1) erfiille die allgemeinen Voraus- 
setzungen unter A. Ferner gelte fiir die Funktion f(t, x, v): 


B. Es existieren stetige Funktionen h(x) und k(t), 7=1, 2,3, derart, daB 
fun alle t, %, v 
(10) i, x, 2) —h(x)| Sh [ x] + AoW [ol + AsO 
gilt. 

C. In jedem Intervall t) <t <t#, existieren nach Vorgabe einer Konstanten 
C >0 positive Konstanten e, c und X derart, daB in fy <t< ¢, fiir alle lesen 


(14) i ee 


Fe = 


gilt, sofern | «| =X ist. 


3H. EnrMANN: Uber die Existenz der Lésungen von Randwertaufgaben bei ge- 
wohnlichen nichtlinearen Differentialgleichungen zweiter Ordnung. Math. Annalen 
134, 167 (1957). Es handelt sich um den 2. Existenzsatz. Die Voraussetzungen A sind 
dort noch etwas allgemeiner. Auerdem verwenden wir hier wegen der Anwendung 


auf Schwingungsdifferentialgleichungen die Bezeichnungen ¢ und x statt * und y 
fiir die unabhangige und abhangige Variable. 
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Dann existieren unendlich viele verschiedene in <f,,¢,) stetige und stetig 
differenzierbare Lésungen von (2), die iiberall dort, wo der Integrand stetig ist, 
der Differentialgleichung (1) geniigen und die Randbedingungen 


(r,) , A, X (to) + ag % (fo) =a, 
(r2) by x(t) + 6, % (4) = 6 


erfiillen, wobei @,, a,, a, b,, b, und b (irgendwelche) Konstanten sind, die ledig- 
lich der Bedingung unterworfen sind, daB sowohl a, und a, als auch 6, und 0, 
nicht zugleich verschwinden. 


Man erkennt sofort, daB alle in der Tabelle 1 rechts auftretenden Rand- 
bedingungen Spezialfalle der allgemeinen Randbedingungen (r,) und (r,) sind. 
Daher k6nnen wir unseren Existenzsatz unmittelbar auf den Nachweis periodi- 
scher Lésungen anwenden, sofern die Funktion / den Symmetrieeigenschaften 
einer der 6 Zeilen der Tabelle 1 geniigt. Da die Funktion /(¢, x, v) beziiglich 
ihres 1. Arguments periodisch mit der Periode TJ ist, geniigt es, die Voraus- 
setzungen fiir ein ¢-Intervall der Lange T z. B. fiir 0<¢<7T zu fordern. Da 
wir noch keine Aussage dariiber machen kénnen, ob kT im Falle k >1 die kleinste 
Periode ist, beschranken wir uns zundchst auf den Fall k=1 und fassen die 
bisherigen Ergebnisse in folgendem Satz zusammen: 


Satz 9. Die Funktion f(t, x, v) erfiille die Voraussetzungen A, B und C in dem 
Bereich OXt<T, —w<x<+to, —o<v<+too. Ferner geniige sie den Be- 
dingungen einer Zeile dey Tabelle 1. 

Dann besitzt die Differentialgleichung (1) unendlich viele verschiedene periodi- 
sche Lésungen x(t) mit der Periode T: x(t+T)=x(t). 

Spezielles Beispiel: Gi gy ORD) ah ys 5. 

Bei der Anwendung auf eingliedrige Schwinger hat die zugehérige Schwingungs- 
differentialgleichung im allgemeinen die speziellere Form 


(12) H4d(t,x,%)%+elt,x) =p), 


wobei d(t, x, x) x das ,,Dampfungsglied", g(¢, x) die ,,Riickstellkraft’ und 4 (é) 
die ,,erzwingende Kraft‘‘ bedeutet, wenn man Dimensionsfragen unberiicksichtigt 
1aBt. 

Beschranken wir unsere Betrachtungen auf diese noch sehr allgemeinen 
Schwingungsdifferentialgleichungen, so kénnen wir aus Satz 9 leicht spezielle 
Bedingungen fiir die Funktionen d, g und # herleiten, um die Existenz von 
unendlich vielen periodischen Lésungen zu sichern: 


Satz 10. a) Die Funktionen d(t, x, x), g(t, x) und p(t) secen so beschaffen, 
daB die Lésungen der Differentialgleichung (12) eindeutig und stetig von den An- 
fangswerten ty, x (to) = %o, % (to) =%y abhangen*. 

b) Die Funktionen d und p seien beschrankt, d. h. es existieren zwet Konstanten D 
und P derart, daB fiir alle t, x, x 
(13) |4(é,%,%)|<D und |p| <P 


gilt. 


4 Siehe Anm. 2 


Q* 
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c) Ferner seien die Funktionen d, g und p in t periodisch mit einer Periode 
L0 
(14) d¢4+T, #4) =dt 44), 204 Dee zg), Dea TY =P). 


d) Es existieren positive Konstanten X, a und o derart, dap fir alle t 


(15) eet) 
gilt, sofern |x|=X ist. Auferdem existiert eine stetige F unktion h(x) derart, dap 
fiir alle t und x 
(16) Ig(é, x) —h(o)|<K,|x| + Ke 
gilt, wobet K, und K, Konstanten sind. 

e) Gelten dann noch entweder die Gleichungen 


(17) d(-t,.%,-v) =-dx0), e(—t =e), b(-9=20 
oder die Gleichungen 

(18) d(—t,—x,v)=—d(t,%,v), g(—t, x) =— eg %), p(-—)=— Pil), 
so hat die Differentialgleichung (12) unendlich viele verschiedene periodische Lé- 


sungen der Periode T. 


Bewets: a) und c) ergeben die allgemeine Voraussetzung A (§ 2). Die Glei- 
chungen (13) und (16) ergeben die Voraussetzung B [Gleichung (10) | des Existenz- 
satzes fiir Randwertaufgaben, und aus (13) und (15) folgt die Voraussetzung C 
dieses Satzes, denn es gilt nach (13) und (45) 


{(t,x,v) _ —d(t, x, v) v—g(t, x) + pid) g Plein lepett fir |x[>X. 
= x Sit: 


x x 


Wahlen wir ¢=o/2, so geht der Quotient auf der rechten Seite wegen a >0, 
o>0 fiir alle |v] <|x|*** gegen — oo mit | x|—> oo. Daher gibt es nach Vorgabe 
einer Konstanten C >0 positive Konstanten ¢(= 0/2), c(=1) und X*, mit denen 
die Gleichung (411) fiir alle |v] <c|x|*** erfiillt ist, sofern | «|= X* gilt. Damit 
ist also auch die Voraussetzung C des Existenzsatzes fiir Randwertaufgaben 
erfiillt. 

Nach Satz 9 brauchen also nur noch die Bedingungen einer Zeile der Tabelle 4 
erfiillt zu sein, damit die Differentialgleichung (12) unendlich viele periodische 
Lésungen der Periode T hat. Dies sichert die Gleichung (17) bzw. die Gleichung 
(18). Ist (17) erfullt, so gilt 


also die Voraussetzung der 1. Zeile der Tabelle 1, und ist (18) erfiillt, so folgt 
H(—t, — x, v) = —d(—#, — x, v)v —g(—4, — x) $A t) 
=d(t, x, v)u+ g(t, x) — p(t) = — Ft, x, v), 


also die Voraussetzung der 3. Zeile der Tabelle 1, w. z. b. w. 
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Spezielle Beispiele: 
1.%+%sint+ xcost+yx2=c-coswi, y>O0, rational, oa” 
T= 2B: ie 


oe . . ’ . n 
2.% +x costsinx + Ginx =c-:sinwt, ae ne ae gan7y. sO sr. 


3.% +%(c, +c. x*— PcosQMt)=0, c,>0, us. 


Zum SchluB diese Paragraphen wollen wir unsere Betrachtungen noch weiter 
auf den am haufigsten untersuchten Fall der Schwingungsdifferentialgleichungen 
spezialisieren, bei dem die Dampfung und die Riickstellkraft nicht mehr von 
der Zeit abhangen: 


(19) % + d(x, %)% + a(x) =p), PE+T) =P). 


Zu diesem Zweck schreiben wir die Tabelle 1 um in folgende Tabelle, die ein 
unmittelbares Nachpriifen der Symmetriebedingungen fiir die Differential- 
gleichung (19) leichter erméglicht: 


Tabelle 2 


Zusatzliche Voraussetzungen fiir die Funktionen 


g (x) P (2) 


Randbedingungen fiir die Periodizitat 
mit der Periode kT 


d (x, v) 


keine 


= —d(x, = p(t) 
ey 4(ET) i(“F*7)\— 
2 eae &(-"2) is) p(—t) . Moca ok ie 
ad (es) =—g(x == 
i (ET) = (*£*7)- 
ee), F Aaa x (0) =2#(£RT) =0, k ungerade 
3 a Be ) : i, oder 
gic oe v) le rane) Ae a «($T)=%(3T+4kT) =0, kungerade 


Aus den vorhergehenden beiden Satzen folgt dann unmittelbar der speziellere 

Satz 11. Die Funktionen d, g und p seien stetig und so beschaffen, daf die 
Lésungen der Differentialgleichung (19) eindeutig und stetig von den Anfangswerten 
ty, X (to) = Xo, % (to) =%q abhangen. 

Fiir die Funktion p(t) gelte mit 2 positiven Konstanten T und P fiir alle t 

pt+T) =p) und |p) <P. 
Fiir die Funktion d(x, %) gelte 
|d(x,%)| <D_ fiir alle x und x. 

Ferner existieren positive Konstanten X, a und o derart, dap 


a = a)|x\° gilt, sofern |x| =X ist. 
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Geniigen auperdem die Funktionen d(x,v), g(x) und p(t) den zusatzlichen Vor- 
aussetzungen einer Zeile der Tabelle 2, so besitzt die Differentialgleichung (19) un- 
endlich viele periodische Lésungen der Periode T. 


Spezielle Beispiele: 
2n+1 


1k ae eo '0 COS. OE ha gy Oe rae alle anderen Konstanten a,, ¢c, @ 
v=0 


beliebig. 


2. # + asin (xx)% + Ginx =fBsinwé mit beliebigen Konstanten «, 8, a. 


§ 4. Subharmonische Lésungen 
Definition: Die Differentialgleichung (1) ¥ =/(t, x,%) habe die Periode T>0, 
d. h. es gelte 
f¢+T, x,%) =f(t, x,%) fiir alle 7, x und~. 


Es sei T >0 der kleinste Wert mit dieser Eigenschaft. 


Dann versteht man unter einer periodischen Lésung der Ordnung k2=1, ganz, 
eine Lésung x(t) von (1) mit der kleimsten Periode kT: 


x(t +kT)=x(t) fur alle ?. 


Ist R=>2, so nennt man x(t) ,,subharmonische Lésung™. 


Bei der speziellen Schwingungsdifferentialgleichung (19) ist also eine sub- 
harmonische Loésung eine solche periodische Lésung, deren Periode ein Viel- 
faches der Periode der erzwingenden Kraft ist. Bei der Untersuchung der be- 
treffenden Differentialgleichungen auf subharmonische Lésungen werden wir 
wiederum einige Ergebnisse aus der Arbeit iiber Existenzsatze fiir Randwert- 
aufgaben verwenden. Wir fassen hier zunachst die betreffenden Ergebnisse in 
die folgenden beiden Satze® zusammen : 


Satz 12. Die Funktion /(t, x, v) geniige den Voraussetzungen A, B und C. 

Dann existiert nach Vorgabe einer Konstanten R>0O und einer natiirlichen 
Zahl N eine Konstante R derart, daB fiir eine Lésungsfunktion x (¢) der Gleichung 
(2) [baw. (1), wo f(é, x(t), #(d)) stetig ist,] mit den Anfangsbedingungen x (fj) = %», 
% (to) =%po, fiir die die Bedingung 


(20) xo + 4 = Ro2R? 
erfiillt ist, folgendes gilt: 
Ate 
(21) (x(¢))? + (4())?> RR? fir “Asta. 


2. Die zugehorige stetige Kurve ¥ = v(x), in Parameterdarstellung x(t), % (t) 
mit ¢ als Parameter, umlauft im Intervall ty <t<t, in der x, x-Ebene den Kreis 
x + %? = R? mindestens N mal, wobei ein Umlauf t, <!< T, +1 dadurch definiert 
ist, daB die Punkte P,(x(t,), #(z,)) und Pyia(%(Ta41), %(tr41)) der Lésungs- 


» Siehe Anm. 3. Es handelt sich dort um Satz 5 von Kapitel IV und um Kapitel V. 


Nichtlineare Differentialgleichungen 135 


kurve x(t) auf der Halbgeraden gy, die von (0, 0) tiber (x , %9) nach Unendlich 
verlauft, liegen und das von der Lésungskurve «x (é), % (é) fiir t,, <¢t< t,., und dem 
Geradenstiick P, P,,, gebildete geschlossene Kurvenstiick den Kreis x? + 42 = R? 
umschlieBt. 


3. Fiir jeden der N, => N Umlaute gilt 


t,—t 
amie saraey et aa 


50), Weed. 0.2, NG de 


4. In jedem Teilintervall von <¢),t,5 der Form a<t<b, b—a=o, <a, 
b<t, liegt mindestens ein Punkt 1,,, wie er in 2. definiert ist. 


Wir betrachten jetzt einen Punkt P=(% 9, %)) auf der ,,Anfangsgeraden“‘ 
(r,) ax +a,%=a 


in der x,%-Ebene. Dieser Geraden ordnen wir einen MaBstab zu, indem wir 
von einem gewahlten festen Punkt P, auf (r,) eine Richtung als positive Richtung 
auszeichnen und einen Mafistab mit fester Einheitslange auf (r,) abtragen. Dann 
bedeute P, > P,, daB die MaBstabszahl von P, gréBer oder gleich der von P, ist. 
Geht P= (4%, %9) gegen Unendlich, so geht auch R§=x§+%§ gegen Unendlich. 
Ferner gilt ® 


Satz 13. Es existiert ein Punkt P’ > P, derart, daB fiir alle P > P’ die Schnitt- 
punkte S = (t,, «(¢,), ¥(¢,)) der Lésungskurve x(é), % (¢) von (2) mit den Anfangs- 
bedingungen x (t)) = %9, ¥ (tp) =%, mit der ,,Endgeraden“ 


(ro) b,x +b,% =b 


sich im Intervall 4) <¢<z#, stetig mit dem Punkt P=(x9, x9) auf (r,) andern. 


Setzen wir die letzten beiden Satze 12 und 13 der zitierten Arbeit als richtig 
voraus, so kénnen wir folgenden Satz iiber die Existenz subharmonischer Lésungen 
beweisen : 


Satz 14. Die Funktion f(t, x, v) erfiille die Voraussetzungen A, B und C in 
acmeporech On tat, Hwa yet wo, =~ocvesco. Es se T 0 de 
kleinste Periode der Funktion f beztiglich t. Die Gleichung (1) besitze keine pertodt- 
sche Lésung mit einer Periode T*+=mT, m=1, ganz. Ferner geniige f den Be- 
dingungen einer Zeile der Tabelle 1. 

Dann existieren nach Vorgabe einer ganzen Zahl k=1 unendlich viele ver- 
schiedene periodische Lésungen der Ordnung k von (1) bzw. (2), deren kleinste Periode 
also kT ist: x(¢+kT) = x(i). 

Bemerkungen: 1. In Satz 14 ist fiir den Spezialfall k =1 der Satz 9 enthalten, 
der sich oben unmittelbar aus dem Existenzsatz fiir die Lésung von Randwert- 
aufgaben ergab. 

2. In den meisten Fallen erkennt man unmittelbar (ohne Rechnung) an der 
Form der Differentialgleichung (1), daB héchstens periodische Lésungen mit einer 
Periode mT, m>=1, ganz, méglich sind. Dies ist z. B. bei den Schwingungs- 


6 Siehe Anm. 3. Der Beweis findet sich im Kapitel V der betreffenden Arbeit. 
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differentialgleichungen der Form (19) der Fall, wenn T >0 die kleinste Periode 
von p(t) ist. Die obige Voraussetzung lieBe sich noch etwas verallgemeinern. 


Beweis: Da die Voraussetzungen des Satzes 9 erfiillt sind, existieren unend- 
lich viele periodische Lésungen mit der Periode T, die die beiden Randbedin- 
gungen einer der Zeilen 1 oder 3 der Tabelle 1 fiir k=1 erfiillen; denn sind die 
Voraussetzungen einer der Zeilen 4, 5 oder 6 erfiillt, so auch die der Zeile 3, 
und die Voraussetzung der Zeile 2 stimmt mit der von 1 tiberein. Wir kénnen 
uns daher auf die Zeilen 1 oder 3 beschranken. Damit sind entweder im Falle 1 
die Randbedingungen % (0) =%(42T) =0 oder im Falle 3 die Randbedingungen 
x(0) = «(4kT) =0 hinreichend fiir die Periodizitat der Losung x(¢) mit der Peri- 
ode kT. 

Es sei nun die natiirliche Zahl k vorgegeben. Dann wahlen wir eine Kon- 
stante C>0 und bestimmen fiir das Intervall O0<¢<kT die nach Voraus- 
setzung C existierenden positiven GréBen e, c und X, mit denen die Ungleichung 
(11) erfiillt ist. Sodann setzen wir R=X und bestimmen nach Satz 12 die 
GréBe R derart, daB die Lésung x(t) in der x, #-Ebene im Intervall 0<f< kT 
den Kreis «2+ %2—R? bestandig umlauft und in jedem Teilintervall ax<t<b 
von <0, kT» der Lange b—a=4T mindestens einen Schnittpunkt 1, mit gp hat, 
wobei T,, und g, die Bedeutung des Satzes 12 haben. Bei unseren Randbedingungen 
fallen Anfangsgerade (r,), Endgerade (r,) und gy mit einer Koordinatenachse bzw. 
Halbachse des x, ¥-Systems zusammen, namlich entweder mit der x- oder mit 
der *-Achse je nach dem Fall 1 oder 3 der Tabelle 1. 


Es existiert dann auf g, ein Punkt P= P’, wobei P’ die Bedeutung von 
Satz 13 hat, derart, daB fiir alle P=P erstens (20) gilt und zweitens sich die 
Schnittpunkte S; der Loésungskurve x(é) mit der Endgeraden (r,) *=0 baw. 
x=0 mit P stetig andern im Intervall O<¢< AT. Fiir diese Schnittpunkte gilt 
wegen (21) und R=X |%,|>X bzw. |x,|>X. Wegen |%,|>0 bzw. (11) fiir 
|x| >X kann die Lésungskurve in beiden Fallen die Endgerade (r,), 4 =0 bzw. 
x=0, auf jeder Seite des Kreises x?-++ %? = R? nur in einer Richtung schneiden. 
Die Schnittpunkte mit (r,) folgen also abwechselnd auf den zwei Halbachsen 
von (rz). Nach Satz 12 nimmt die Zeit zwischen je zwei dieser Schnittpunkte 
mit wachsendem Ry also auch mit wachsendem P ab. Betrachten wir den ersten 
auf¢—=3T folgenden Schnittpunkt mit (r,), 4, > 47, so kann daher bei wachsendem 
P nicht stets ¢,>3T7 bleiben, sondern es muB fiir ein P>P ¢,=3T werden’. 
Dann sind beide Randbedingungen der Zeile 1 bzw. der Zeile 3 der, Tabelle 4 
erfiillt. Wir haben eine Lésung mit der Periode T. Durch weitere VergréBerung 
von P erhalt man eine weitere Lésung, usw. 

Von einer solchen Lésung x(t) =x(t+ 7), fiir die P= P>P ist, gehen wir 
Jetzt aus. Natiirlich hat sie auch die Periode RT. AuBerdem gilt fiir sie wegen 
der Periodizitaét im Falle der Zeile 1 (Tabelle 1) (nT) =0, n=0, 1, 2, ... bzw. 
im Falle der Zeile 3 x(n} T)=0, n=0,1,2,... Wegen der Symmetrie zu den 
Geraden t=n3zT nach Satz 2 baw. zu den Punkten t=n4T7, x=0 nach Satz 3 
(in der ¢, x-Ebene) liegen in jedem der Intervalle n3T <t<(n+1)4T, n=0,1, 
2,... genau gleich viele Schnittpunkte S, der Lésung x(¢) mit der Endgeraden 


? Nach diesem Prinzip wurde der zu Anfang von § 3 angegebene Existenzsatz fiir 
Randwertaufgaben bewiesen. Vgl. die in Anm. 3 zitierte Arbeit. 
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(r.) (in der x, #-Ebene), d.h. Stellen % (¢)=0 bzw. x(t)=0. Ihre Zahl in einem 
der Intervalle n3T<t<(n+1)$T sei s. Dann liegen im Intervall 0<t<kiT 
genau ks Schnittpunkte mit (r,). Der Schnittpunkt an der Stelle kT ist der 
(ks-+-1)-te. Wir betrachten sodann den (ks+2)-ten Schnittpunkt S,..,, fiir 
den also gilt: ¢,,:.>R4T und vergréBern nun P iiber P hinaus. Dann dndert 
sich nach Satz 13 S,,.. und damit auch ¢,,,, stetig mit P. Nun kann nach 
Satz 12 bei fortwahrender VergréBerung von P und damit auch von R, nicht 


standig t, 4. >hk4T bleiben, sondern es muB einmal fiir ein P = P~P Eee Sra 


sein nach denselben Uberlegungen wie oben. Fiir dieses P sind wiederum beide 
Randbedingungen % (0) = x (kT) =0 bzw. x (0) = x (4 RT) = erfiillt. Nach Satz8 
hat daher die Lésung «(t)=~(t, P) die Periode kT: x(t-+kT)=<x(t). Sie be- 
sitzt offenbar im Intervall O<¢<4$kT genau ks+1 und in 0<t<kT wegen 
ihrer Symmetrieeigenschaft beziiglich des Punktes (ERT, 0) (Satz 3 und 6) bzw. 
der Geraden t=}hT (Satz 2 und 5) genau 2ks+2 Schnittpunkte mit (r,). 

Wir haben nun noch zu zeigen, daB kT die kleinste Periode unserer so kon- 
struierten Lésung ist. Zu diesem Zweck nehmen wir an, es gabe eine kleinere 
Periode kT. Dabei sei RT die kleinste Periode unserer Lésung. Dann muB k 
ein Teiler von k sein, was wir jetzt zunachst zeigen wollen: 

Nach Voraussetzung ist R>1 eine ganze Zahl. Mit RT und RT ist auch jede 
GréBe T* = T(c,k —cgk) mit ganzen Zahlen c, und cy Periode von x(t). Die 
Zahlen c,k—cyk sind aber gerade die Vielfachen nd, »=0, +1, +2,... des 
groBten gemeinsamen Teilers d=(k,k) von k und k8. Damit ist auch dT Periode 
von x(¢) mit d<k und d Teiler von k und k. Aus der Minimaleigenschaft der 
Periode kT folgt daher k=d und damit, daf & ein Teiler von k ist: k=qk mit 
ganzem g. Da nach Annahme k <k ist, folgt g= 2. 


Wir zeigen jetzt, daB q keine gerade Zahl sein kann: 

Annahme q ware gerade: k=2mk. Dann ware auch 4kT =mkT Periode 
unserer konstruierten Lésung x(¢). Nun ist aber ¢=3kT nach unserer Kon- 
struktion der (ks-+ 2)-te Schnittpunkt der Lésung mit der Geraden (r,). Da 
(r,) und (rg) zusammenfallen, ist 4=0 der 1. dieser Schnittpunkte. Diese legen, 
wie wir oben gesehen haben, abwechselnd auf der einen und der anderen durch 
den Kreis x2-+-%2— R? getrennten Halbgeraden von (r,). Da k=2mk als gerade 
angenommen wurde, liegt der Schnittpunkt ¢,,..—3kT auf der anderen Halb- 
geraden von (r,). Die Lésung x(f) kann daher nicht die Periode }kT =mkT 
und also auch nicht die Periode kT haben. 

Es sei jetzt g ungerade, wegen k<k also g=3. In jedem Intervall der Form 
nkT <t<(n+1) kT liegen gleich viele Schnittpunkte S; von x(é) mit (r,). Ihre 
Anzahl in einem solchen Intervall der Lange kT sei 7. Dann legen im Intervall 
O<t<gkT=hkT genau qr solcher Schnittpunkte. Nach den obigen Betrach- 
tungen sind es andererseits genau 2ks+2. Es miiBte daher gy =2ks +2 sein, 
wobei die GréBen g,7, k und s ganze positive Zahlen sind und g >3 Teiler von k 
ist. Dies ist aber nicht méglich. Damit ist kT in der Tat die kleinste Periode 
unserer Lésung, w. z. b. w. 


8 Vel. z.B. Hasse, H.: Vorlesungen tiber Zahlentheorie, S. 24, Satz VI. Berlin 
1950. 
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Genau wie in §3 bei Satz 9, 10 und 11 1a4Bt sich auch der Satz 14 leicht auf 
Schwingungsdifferentialgleichungen spezialisieren. Es ergeben sich dann die 
beiden Satze: 

Satz 15. Fiiy die Differentialgleichung (12) gelten alle Voraussetzungen des 
Satzes 10. Ferner existiere keine periodische Lésung mit einer Periode T* == mT, m 
gan. 


Dann besitzt. (12) unendlich viele subharmonische Losungen jeder Ordnung. 


Satz 16. Es gelten fiir die Differentialgleichung (19) alle Voraussetzungen des 
Satzes 11, und es sec T>0 die kleinste Periode der Funktion p(t). 
Dann besitzt (19) unendlich viele subharmonische Lésungen jeder Ordnung. 


Beide Satze folgen unmittelbar aus den vorherigen. 
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1. Introduction 


The parametric theory of ordinary differential equations has achieved a high 
level of development. In particular the asymptotic theory has been studied 
extensively by many authors, including G. D. BrrkHorr, LANGER and others. On 
the other hand, the parametric theory of partial differential equations has barely 
begun to be investigated. 


The existing literature on this subject falls into several categories. In one 
of these categories an appropriate limit partial differential equation is associated 
with an original partial differential equation. The original equation contains a 
parameter which is considered large. The limit equation is usually but not 
necessarily the formal limit of the original equation as the parameter tends to 
infinity. A statement is then made about the difference of appropriate solutions 
of the original and limit equations as the parameter tends to infinity. In 1922 
W. STERNBERG [6] considered the equations 


Cu C pe ea 
(1) x Oy + [C(x,y) + #]u=0 
and 
ow 
(2) 6x Oy + k?w=—0. 


STERNBERG shows that 


where u is a solution of (1) and w is an appropriate solution of (2). In a second 
paper STERNBERG [6] shows that u—w=0(,,] for the pair of equations 
Au + [A(x,y;k) + k?]u=0 
Aw+k?w=0. 


Here A is bounded uniformly with respect to k. In 1944 W. Wasow [8] considered 
the pair of equations 


Au+ko" =k} (x,y) 
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In 1950 N. Levinson [4] considered the pair of equations 


é ou 
5 Au + A(x,y) 5 + B(x,9) gy + € (x,y) 4 = D(x,9) 


5 


0 
A(x, y) ge + B(x, 9) Gy + C(% 9) @ = D(x,9). 


In all of these results the domain of definition D of the solutions is bounded. In 
STERNBERC’S results there is no discussion of boundary values of the solutions. 
In the results of Wasow and LEvinson it is shown that the solution w of the 
limit equation takes on the values of the solution w of the original equation on 
appropriate segments of the boundary 8 of D. This is as much as can be expected 
since the original and limit equations are of different orders. 

A second category of the literature deals with formal series expansions of 
solutions. In 1933 G. D. BirKHorF [1] suggested a procedure for obtaining formal 
series expansions of solutions of the reduced wave equation 


(3) Aut ku=o. 


co 
The procedure is to insert the series aghey GAP into (1) and equate coef- 
0 


ficients of the various powers of k to zero. This gives rise to a recursive system 
of first order partial differential equations for the determination of the coef- 
ficients Y% and o,,, of the terms of the series. BIRKHOFF conjectured that there 
exists an appropriate solution of (3) to which the formal series is asymptotic. 


In 1944 W. TryjiTzinsky [7] generalized BIRKHOFF’s procedure by con- 
sidering a more complicated formal series and a linear partial differential equation 
of the m” order. The domain © of definition of the solution of the equation 
is bounded. Tryjirzinsky shows that under appropriate conditions, the formal 
series exists in a subdomain Q, of D. In the case that the equation is of elliptic 
or hyperbolic type, he shows also that there exists an appropriate solution of 
the differential equation to which the formal series is asymptotic in a subdomain 
De Ol Dye 


In 1954 M. KiinE [3] considered the hyperbolic wave equation Au = ae 
from which Au-+k?u=0 is reduced. He shows that if a solution of Au — ae 


has certain jump discontinuities, then the corresponding solution of Mu + k2u =0 
possesses an asymptotic expansion of the type conjectured by BIRKHOFF, 

In 1956 J. KELLER, R. Lewis & B. SECKLER [2] applied the procedure of 
Birkuorr to a large number of specific diffraction problems. They compare 
their series with those obtained from the exact solution of these problems. In 


every such comparison there is complete agreement between the known correct 
result and the formal series. 


In this paper we investigate the validity of the method suggested by BIRKHOFF. 
We consider a solution w of the reduced wave equation defined in the exterior 
® of a bounded or unbounded domain. The values of w on the boundary %$ 
of D are prescribed. w satisfies the radiation condition given below as equation (5) 
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In Section 2 we shall illustrate the method of BirKHorr and show how to 
obtain a formal series expansion for u. The discussion in this section follows 
closely some of the material in [2]. In Section 3 we define the notion of regular 
domain and then establish two lemmas upon which the proofs of the theorems 
of Section 4 will depend. In Section 4 we discuss the asymptotic properties of 
the formal expansion. We show that the expansion is asymptotic to a solution 
of the reduced wave equation. Some theorems are given which illustrate the 
relationship between this solution and the exact solution of the boundary value 
problem. In the case that the domain is regular and is the complement of a 
certain class of bounded domains, we use the estimates of [5] to show that the 
formal expansion is asymptotic to the exact solution of the boundary value 
problem. 


2. The Method of Birkhoff 


We restrict ourselves to the case of two dimensions in order to achieve clarity 
in the manipulations to follow. 


Let (x,y) be a solution of the reduced wave equation 
(3) Aut k?u=0 


in the exterior, D, of a bounded or unbounded domain. Let %, the boundary 
of D, be given parametrically by x= x,(t), y=ypo(t). ¢ is arc length along %. 
As one moves along % in the direction of increasing ¢, the domain D should lie 
on the right. On % let « be prescribed by 

M 


1 an (¥, V) 
(4) Uy = ee oe, 9) Ds (i kh)? ise 


n=0 


Finally let w satisfy the radiation condition 


(5) lim {| —ieultas= | 


where Gp is a circle of radius R with some fixed point in space as center and 
where ds is a line element. 

To find the asymptotic expansion of u, the procedure of BIRKHOFF is as 
follows: Assume that w has an asymptotic expansion u™ of the form 


On(%, ¥) 
lo) ik Y (x, y) n\*> J 
(6) Ur = € 2, Ga 


We insert “© into (3) and equate the coefficients of various powers of k to zero. 
This leads to the following system of equations for Y and o,, n=O, 1, 2,... 
(7) Uae) 

(8) VG OV E64" Y= G7, Ca 


Thus Y and g, are solutions of first order partial differential equations. Their 
boundary values are obtained by equating the value of the formal expansions 


* When it is necessary to stress the fact that a function is evaluated on the 
boundary $ of D, we shall use the subscript 0. 
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on % to the boundary values, (4), of w. We obtain 
(9) y=—@ 

w= 0,48.,M 
oe =0, n=M+1, M-+2,.... 


Thus the formal expansion u® is equal to w on &. 

To solve the eiconal equation (7), we let %=#, Y%=gq and use the method 
of characteristics. We denote the arc length along az characteristics by s. 
This method yields the following system of ordinary differential equations 


tee dp 
ei Oe 
hoe ay Te 
ey Pe: he erent 
14) 
Oe e+ =1 
d®@ iG, dV o 


From (14 f) and (11g) we find 


d® d dy, 1/, ad® \2 
a) Pea baba fa os [4 ( 


dt dt dt 

(12) dD dyy dx dP \2 
= Oy bs 

DE dt dt ~ di \: | a : 


From (11a)—(11e) we find 
a) % = PoS + Xp 
b) ¥=4oS + Yo 


for the characteristics or rays, and 


(14) Y=st@. 


(13) 


The recursion equations (8) may be solved as follows. From (7) we see that 
the vector V¥, normal to the surfaces ¥Y% =constant, is a unit vector. From 
(11a) and (11b) we see that the rays are also normal to this surface. Thus 


Vo, -V¥ is the ordinary derivative os of o,, with respect to the length of a ray 


along each ray. Thus (8) is just an ordinary differential equation along each 
ray. The solution of (8) along a ray is 


ai f pewas V?o,_,(t) dt 


0 t 


(16) 9,(6) = 0,(0) [20 © f G-1( P29, (0 at. 


G(s) is the curvature of the surface Y%—constant at the point s along the ray 
emanating from (x (é), yo(t)), and o, (0) =a, (x ee Vo (t)). 
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We shall introduce new coordinates (s’, 8) in place of (s, ¢). s’ is the distance 
along a ray measured from the envelope or caustic of the rays. f is the angle 
which each ray makes with the x-axis. These coordinates are called caustic 
coordinates. Let sy(¢) be the distance along a ray from the caustic to 8. s,(t) 
has the same sense as s. Then (s’, B) are related to (s, t) by 


17 ss=s+s (i), P= tan2 fo. 
The regularity of the caustic coordinates depends on 8 and ©. We now drop 
the prime. In the new coordinates (16) and (14) become, respectively, 


So \2 1 il 
(18) on (6B) = on (508) () — 54 [ 470,16) de 
(19) Y—=s—s,— @. 

Let the caustic be given by 
(20) ~=§(6) y=n(). 
Let a dot denote differentiation with respect to B. Let 
Qe &2 + 7? 

(21) b = Ecosf +77 sinB 


C = £cosB +ijsin Bp. 


With these abbreviations, we may express the o, in caustic coordinates as 
follows *: 


(22) 01, (8, B) = > fin (8) SP 


5 1 F 1 ce ; O 
lin= ae i ah 1,” ee 1,” gi 20 eel) tae SE 
; 3 ‘ Sais 1 : 
tol — 3) h-2 nA ali (i Slee fee 1+ 0, 
3n 
ps SOS ey) n= 0, 
j=1 


foo = %(S0-B) Gi 


3. Prerequisites 


In this section we discuss formal details concerned with the existence of the 
formal series and the regularity of the caustic coordinate system. We shall 
introduce the notion of regular domain. Finally some lemmas will be proved. 


A. Definition and Assumptions. #, and q will be real for all values of ¢ 


; ad@® \2 
provided the radicals in (12) are real. This requires (ar) <1. If the boundary 


values of wu arise from an incident field, this condition is automatically satisfied. 


* For details concerning the derivation of these expressions, see [2]. 
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For in this case, the phase of the incident field satisfies the eiconal equation (7). 
From (13) and the ambiguity of sign in (12), we see that each value of ¢ gives 
rise to two rays. Thus as s increases from sy, two rays will emanate from ¥ 
and enter D or its complement. In principle, we could use either sign and the 
corresponding ray to solve for the coefficients ¥ and the o,, »=0,1, 2,..., of 
the formal series. Indeed we could use one sign for certain values of ¢, and the 
other for the remaining values. 

Now the characteristic differential equations and the choice of sign in (12) 
assign a positive sense to the arc length s on each ray. If as s increases from Sp, 
a ray enters 9, then along this ray the formal expansion behaves like e'*!*!, 
If the ray enters the complement of D as s increases from sg, then along the 
part of this ray which enters D, the formal expansion behaves like e~**'*!. Since 
a solution of Au--k?u—0 which satisfies the radiation condition behaves like 
e’*”* the formal expansion could not be valid along the second type of ray. If 
for each value of ¢t one of the rays emanating from % were of the former type 
while the second of the latter type, this observation would serve to determine 
a single ray for each ¢ along which to compute Y and the o,. The following 
lemma and remark assert that this is indeed the case. 

a(x, 9) 
d(s, t) 
as s increases from Sy, and conversely. 


Lemma. 7 - ->0 for a value of t, then the corresponding ray enters D 


Remark. UE a = 1 (F i the upper and lower signs corresponding 
to the upper and lower signs in (12). 

The proofs of these facts follow from simple computations using (11), (12), 
and (13) and therefore will not be given. 

TP) 
0(s, t) 
that the radicals in (12) vanish. In this case the ambiguity in sign disappears 
and the procedure for constructing the series is well defined. 


= 0 for s=sy and certain values of ¢, the preceding remark asserts 


We introduce the notion of a regular domain in the following: 


Definition. A domain D with boundary & is regular relative to D, if the 
portions of all rays for s>Ss, form a field in D, if ss> e>0, and if 


max {]€)I1n)||E A111 BEB) HB) <4. 


éand A are positive numbers. 


Hereafter we shall consider only the portions of the rays which correspond 
to positive values of s. 


Examples of Regular Domains 


i) The domain to the left of the parabola x—=y? is regular relative to the 
phase @= x (plane wave). The caustic in this case degenerates to the focus of 
the parabola. . 


ii) The domain exterior to a bounded smooth convex curve is regular relative 
to ®=constant. 


wae I. 
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ii) No domain which is the exterior of a bounded domain is regular relative 
to ®=x. If the boundary is smooth, then by deleting appropriate portions of 
Da regular subdomain of ®D is obtained. 


From (22) and (23) we see that o,,(s, 8) is a function of 


2n 2n—2 
Oa Gla. Oa, 0 ay 


(24) do. ap? Ogee Be BB? gpeame 


LSS 
Thus the amplitudes a,, ~=0,1,..., N must be sufficiently often differentiable 
for the formal expansion to exist to N terms. Hereafter we assume that these 
quantities are as many times differentiable as is sufficient to justify the arguments 
which follow. 


B. Two Lemmas. The relationship between caustic and cartesian coordinates is 


(25) x= E(B)-+scosB, y=n(8)+ssinB. 
The Jacobian of this transformation is 
(26) J(s,B) = sin B — 7 cosB — s. 


The Laplacian in caustic coordinates is 


eke a\ 2 1 c GN\ ROR AU OR, a ihe C2 Ore 
(27) V =(1 E =) as | e s? s) és s? @s 0B ire ap a s3 ap - 


a,b, and c are defined in (21). We now state and prove 


Lemma 3.1. Let D be regular relative to D. Let the quantities occurring in 
(24) and the first and second derivatives of these quantities be bounded. Then 


(28) |4o,| and {| 4o,| dS 


are bounded. 


Proof. Applying (27) to (22), we obtain 
3n 


: 1 1 
(29) V¥0,= 9 Dy 


= 


\e 1 9 
HS +4) tint hint e(S +3) bn t 
CE tA) fin t 201 +3) bin, (EFAE FI Gf jn) 
Ss 


s2 


Since D is regular relative to ®, s=> e>0 at every point in . From the recur- 
rence formulas (23) and the hypotheses, all quantities in the right member (29) 
are bounded. Thus |V”2a,,| is bounded in D. From (29) we have a constant C, (m) 
such that 


C,(n 
(30) |V?o,,| <a" 
Then 7 
2m ©co 7% OO 
n dpd 
G1) ff |Va|axays ff OP lrepldsapscm [fA, 
5 ore Ons 


from (26). Here C(m) is some constant. This proves the lemma. 
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We now state and prove 
Lemma 3.2. Under the hypothesis of Lemma3.1, u% satisfies the radiation 


condition. 
Proof. 


N 
N _ jik(P—S) 2 > ee _ 
(32) UE Me 4 1 (7 \p 


Denoting differentiation with respect to 7, B, and s by subscripts, we have from (32) 


rs # 
ul’ —ikuN — uf) B,+ uy s,—tkhu 


3n 3 

ik (D— = 

iks N Ds {e" ( So) Tinks! 
é fe 


es Ga B, + 
(33) ee Sea Are 
mee pa, Gk)" (s,—1) 
hen, Ze MOA p< 
a 5 pls: Gh? cm 
From (25) we have 
(34) ye EU t pet st 25(y sin B + € cos p) 
and 
(35) tan? = ae 
From (34) we have 
(36) lim — = 45 


Differentiating (34) and (35) with respect to 7, we obtain 
y = B,(éé +H + s(7sinB +7 cos B + Ecos B = &sin B)) — 
+ s,(s + ysin B + & cos B), 
0 = B,((m + s cos B) (§ + s cos B) — (E — ssin B) (n + ssin B)) + 
+ s,(sin B(E + s cos B) — cos B(y + ssin)). 
Upon solving (37) for B, and s,, we may deduce, using (36), that 


(37) 


(38) B,=0(—] and” $)-<S4 "Vas — S=sTost 

Then from (33) we see that 

(39) a cee al (1/2) ob 
Ss ye 


which proves the lemma. 
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4. The Asymptotic Properties of the Formal Series 


In the first part of this section, we consider domains D which are regular 
relative to ® (see Definition 3.1) and whose boundaries 8 are bounded. These 
conditions will not be mentioned specifically in the hypotheses of the theorems 
which are stated. In the second part of this section we shall indicate the extent 
to which these results are valid when 8 is unbounded but © is still regular. 
Finally in the third part of this section we shall consider the case when ® is 
not regular relative to @. 


A. Regular Domains with Bounded Boundaries. We shall need the following 


Lemma 4.1. Let w satisfy the radiation condition and be a solution of 
Au+ku=fin®. Let|f| and ff |f| dS be bounded uniformly with respect to k. 
D 


Then for P in the closure of D, we have 


(40) |w(P)+ 2 f |p HY | P — Ol) —w HY | P — Q))|as| so). 
B 


H\ is the Hankel function of the first kind and of the zero“ order. 


Proof. Since w satisfies the radiation condition and JJ |f| dS exists we have 
from Corollary 1 of Chapter 1 of [5] 


(41) 4iw(P) = [52 ele — Ql) we HY) (k| P — Q))] as 


ae ee El i Q|) as 


We now show that the domain integral is bounded by O(k~*). The Hankel 
function decays like (k| P—Q|)~? as k| P—Q|—o, and has an integrable 
logarithmic singularity at P=Q. These facts and the hypotheses show that 
the domain integrand is absolutely integrable. Now break up D into a bounded part 
containing P and its complement. In the unbounded part we have | P—Q|> 
constant. Thus the integral over this part is bounded by O(k~!). STERNBERG [6] 


has shown that 


(42) {S| Ao \(k| P — Q|)|dS<O(k-), 


where Y is a bounded domain and P is any point in the closure of Y{. The constant 
occurring in the order relation (42) depends on the area of 2. Since | /| is bounded 
uniformly with respect to k, (42) completes the proof of the lemma. 

Our results are now given in the form of theorems. 

Theorem 4.1. wu, the formal series to N terms, approximates a solution of the 
reduced wave equation to the order O(k-\~*). This solution is the exact solution 
of the boundary value problem plus a remainder. The boundary values of this 
remainder are at least as small as this degree of approximation. 

Proof. Let 
(43) kNoN =u— wu, 


10* 
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From the recurrence equations (7) and (8), we obtain 
(44) Aw’ + k2wN = (— iN e'** Ao,. 
Lemma 3.1 implies that /=(—1)%e'* "Aa, satisfies the hypothesis of Lemma 4.1. 


Lemma 3.2 implies that w™ satisfies the radiation condition. Thus if 


N 4) 
28” HE) (k| P — Q|) — w' HP (k| P— Q))] as, 


(45) 4eaN RY = f[ 
% 


Lemma 4.1 applied to w’ implies 


(46) | wo” — RN R®| <O(k-5). 
By (43) this is 
(47) |u — uw — RX|<O(k-*-*). 


Since R% is a solution of the reduced wave equation, (47) proves the first two 
assertions of Theorem 4.1. The last assertion of the theorem follows by con- 
sidering (47) for a point on 8 and the following inequality which is implied by 
(9) and (10). 

0 N2=M 
fh = O(N Sy SN. 


(48) Jmol S) 
kh 


N+1 


Theorem 4.2. u% approximates to the order O(k~X~5) a solution of the 
reduced wave equation composed of a single and double layer distribution. The 
densities of these distributions are respectively the values of the normal derivative 
of uN on & and the value of u% on &. 

We know that the boundary values of wv’ are the same as the boundary 
values of u for N= M. For N<M we have from (10) that | #9’ —u)| <O(k-~?). 
The following relationship persists between the values of the normal derivative 

‘0 


N wal ee. Ou 
of uw on &, ae and the values of the normal derivative of u on %, ai 


Theorem 4.3. aa satisfies the integral equation of the first kind for 
to the order O(k-N~5) for N= M or to the order O(k-N—3) for N<M. 

We now furnish the proofs of these two statements. 

Proof of Theorem 4.2. Since Au+k?u=0, (44) becomes 


OU, 
“On 


(49) AuN + k? uN = rae AR DAs 


From Lemmas 3.1 and 3.2, we see that we may apply Corollary 1 of Chapter 1 
of [5] to vw’. This gives 
(50) 47u% (P =f{ eu HO (k| P wn _@ FQ) 
Ne Ayes 0 (f | QO \)i => & Bn 0 (&| P— Ql)|as — 
B 
— gar {fet 40, @|P— elas. 
‘D 


In the proof of Lemma 4.1, we pointed out that the domain integral in (40) 
was bounded by O(k~5). Thus (50) implies the theorem. 
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Proof of Theorem 4.3. In (50), let P approach 8. In the limit we have for 
Ped 


auN é 
(54) siu¥(P)— [ F HM(a|P— Ql) ast fw 2 HM @| P— Ol) ds 
B B 


e*” Ay HO) (k| P—Q|) aS 


Using (9) and (10), we may rewrite this as 


5 ou 1 i 
(52) ete ce Aka = cet ae 


a ik = ne ik Ay 
oe Dae (ik)” ae ay (ck)” alk "ds. 


If N= M, the last two terms in the right member are zero. Thus (52) and the 
estimates for the domain integral derived in Lemma 4.1 imply 


(53) sin(P) — fe AD hee Hp ds soe 
B 

This proves the theorem in this case. For N<M we must show that the last 

two terms of the right member of (52) are bounded by O(k~*~4). The first 

of these terms is O(k~%~*) so that we need only consider the integral. We have 


ik ay < Noah) HO) 
(54) fe Dy ae git ) ds /so0r ) tan | as. 


By the methods used by STERNBERG [6], we can show that the integral in the 
right member of (54) is bounded by O(k*). This and (54) complete the proof 
of the theorem. 

The above theorems have shown that w% is asymptotic to a solution of the 
reduced wave equation. Among other properties of this solution, we have seen 
that it is a sum of the solution of the prescribed boundary value problem and 
a remainder whose boundary values are small. In Chapter III of [5] it is shown, 
under appropriate conditions, that a solution of the reduced wave equation is 
small with its boundary values in a certain sense. With this result we may 
sharpen our statements. Thus we have the following theorems: 


Theorem 4.4. Let 8 be a smooth simply closed curve and let Imk>0. If 
2Im k is larger than some constant C(%) then 


O(|R| N+) 
(Imk)? * 


uN =u + 


Theorem 4.5. Let % be a circle of radius a, let k be real, and let k> - Then 


uN’ —u+O(k Nt), 
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Proof of Theorem 4.4. From Theorem 4.1 we have for Im k>0 


(55) uN =u + RN + O(|k|-%-*)) 
and 
(56) | Re |s<O(e|-*~). 


RN as given in (45) is a sum of layers. A computation with (45) shows that 
R% satisfies the radiation condition. Thus by Theorem 3 of ChapterI of [4] 


RN a 

67) RP =+ f [HP e|P—o) -— RS Apel P— Q)| as 
B 

From (57) 


58 lub| RY (P)| <* | lub| RY (P 
4 


Pe 


) (k| P— Q|) dst 


B 
sc a PN 
+|/ fixeras|/ f aa as), 
B B 


Let C(%) be larger than the right member of (3) of Chapter III of [5]. Then the 
estimate of Lemma 2 of Chapter III of [5] is applicable. Inserting the estimate 
of the lemma into (58), we obtain 


(59) lub | RY(P)| <5 ]O(|R|-*- ) {ban a (k| P — Q)) dst 


of fiseieeo ea HE + aft) fiers fl EPa 


; we proceed as follows. Setting f=(—1)%e'**”Aoy in (41), 


N 
To estimate eS 


we obtain 
(60) Aiw’ (P) = 47k% RN (P = yan (—1) WE Ao ae (k| P — Q|)d 


From (43) and (9) and (10), (60) becomes 


M 


(61) RY (P) = — fk¥ Ya colle e¥ Aoy HO (k| P — Q|)d 


Gay e a 


From this, since by Lemma 3.4 | Aoy| is bounded, we obtain 


(62) 7 RX) soda) 400K) ff 2 apelP—opas 
D 


Using the technique of STERNBERG’s estimate [6], we can show that 


ra) 
eee ay 
f on Hy 


(63) 
[| HP |2as<0(R|-# 
S 
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and 
| 8 p71) 
(64) J} hee 
D 


Now insert (64) into (62). Insert the resulting estimate and (63) into (59). We 
obtain 


(65) lub| RY (P)| <0 (||-*—8) O(| |) + 


oars | Stel) 4+ oda om as 


es O(|R| NTH) 


dS <0 (\h|*) 


<O(|k|-'**s) ae O 


(Imk)? ~~ (Im)? 
Inserting this into (56), we obtain 
we OUR) 
(66) Us =u + “(mA 


which proves the theorem. 

The proof of Theorem 4.5 is parallel to the proof of Theorem 4.4. The estimate 
of Lemma 3 of Chapter III of [5] replaces the estimate of Lemma 2 of Chapter III 
of [5]. The proof of Theorem 4.5 will therefore be omitted. 


B. Regular Domains with Infinite Boundaries. Since the estimates of 
Chapter III of [5] are not known to be valid when % is infinite, Theorems 4.4 
and 4.5 are not necessarily true in this case. 

An examination of the proofs of the remaining theorems of part A of this 
section reveals that they are correct if the representation of Corollary 1 of Chapter 1 
of [5] for wv and w’ is applicable when & is infinite (see (41), (50), and (51)). 
Theorem 4.3 requires one additional consideration which we shall mention below. 
Corollary 2 of Chapter I of [5], which provides this representation when % is 
infinite, requires that certain boundary integrals exist. By the remark preceding 
Corollary 2 of Chapter I of [5] these Pa will converge and the proofs of 
Theorems 4.1, 4.2, and 4.3 will be valid for % infinite if 
| aud | 


| on, i < constant. 


OUy 
én (2 


N 
Ug 


(67) Jim R max { 
Now in the proof of Theorem 4.3 in the case that N<M we used (see (54)) 


(68) [ dae : "(| P — Ql) | dss o(), 
N¥1 
(67) implies (68) for % infinite. 
We formulate the result of this discussion as 
Theorem 4.6. Under the additional condition (67), Theorems 4.1, 4.2, and 4.3 
are valid when & 1s infinite. 


For N= WM we may make the following stronger statement. 
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Theorem 4.7. If the representations of (41), (50), and (51) are valid, then 
Theorems 4.1, 4.2, and 4.3 are true when B is infimte. 


C. Non-regular Domains. A domain 2 which relative to 2 is non-regular 
may be divided into a regular domain t and a non-regular domain %. D= RMN. 
In the following figure the exterior of a circle with CALs, a9=1)6,=—0, 7 en 
is decomposed in such a manner. On the left side of the circle we choose the 
upper sign in (12) while on the right side we choose the lower sign. The caustic 
of the former rays is indicated in the figure. The caustic of the latter rays 1s the 


a 


Bead caustic 
of rays 
emanating 
from lett 
side of 
circle 


Fig. 1. Example of non-regular domain 


point at *== — oo. The common boundary of % and %t are the two rays L. The 
common boundary is chosen in this manner in order to make - 7, €,7, bounded 
in ft. 

Let S,=ROB. Let T=LUBy be the boundary of R. 

Our theorems are inapplicable to %t because the caustic of the rays in % is 
at ~=-— oo. However, the domain it with boundary & falls into the case of 
part B of this section with some contingencies. The first of these concerns the 
representations used in (41), (50), and (51). The integrals over 8 are replaced 
by integrals over Y. Some assumption must be made about the convergence 


of these integrals. The restriction of (67) would suffice but is artificial since & 
is in the interior of D, the domain of definition of w. 


In Theorem 4.1, the conclusion that the boundary values of the remainder 
solution R” are less than O(k~N~*) must be replaced by the conclusion that 
the values of RY along Bp are of this order of magnitude. 


With these alterations Theorems 4.6 and 4.7 are valid for the domain Q. 
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Uber das Newtonsche Verfahren 


JOHANN SCHRODER 


Vorgelegt von L. COLLATZ 


Es werden Iterationsverfahren der Gestalt 
Un+y = T (Uy; 4) (n == Q; ie 2, med a (0.1) 


in einem (verallgemeinerten) metrischen Raum untersucht. Die x, seien dabei 
nach einer gegebenen Vorschrift aus den vorangehenden u,; bestimmt. Diese 
Form hat das (allgemeinere) Newtonsche Verfahren 


Unti = Uy, — [Gixa)] ve Gu, (n = 0,1, 2,.. .) (0.2) 


zur Lésung einer Gleichung Gu= 9 in einem Banachschen Raum. Spezialfalle 
sind das gewohnliche und das vereinfachte Newtonsche Verfahren (x, =u, bzw. 
x, =w, fest). Es interessieren jedoch auch andere x,. Das Iterationsverfahren 


Pu+1 = Mn+1 A Pn> L®y =) 4 (n = (6) ‘he pee oe .) (0.3) 


(1 geeigneter fester Zeilenvektor) zur Bestimmung des gréBten Eigenwertes 1/1 
einer Matrix A und des zugehérigen Eigenvektors m laBt sich z.B. als Newton- 


sches Verfahren der Form (0.2) mit den Elementen wu, = tS) x= ( ‘A deuten 
(§§ 4,5) : 


Gesucht sind Elemente u*, fiir die 


n 


u* —lim T(u*,x,) oder u*=T(u*,x,) fiir allen =4,2,3,... (0.4) 
gilt’. Jeder dieser Bedingungen entspricht beim Newtonschen Verfahren (0.2) 
unter geeigneten Voraussetzungen Gu* =, 


Die Untersuchung des Verfahrens (0.1) (in § 2) geht aus von einer Abschatzung 
des Abstandes 6 (T(w, x), T(v, y)) durch eine Funktion F mehrerer Variabler [vel. 


* Vortrag auf der GAMM-Tagung 1957 in Hamburg. 

* Auch beim gewohnlichen Newtonschen Verfahren (x,,=u,) erscheint es zweck- 
maBig, die Lésung so zu charakterisieren. Zundchst lage es nahe, dies Verfahren in 
der Form u,41= Tu, mit Tu=T(u, uv) zu schreiben und dann nach Elementen w* 
mit w* = T'u* zu fragen. Jedoch soll der Fall G(,*)=O nicht ausgeschlossen werden, 
fiir welchen Tu* = T(u*, u*) nicht erklart ist. Che perden deutet auch der Ein- 
deutigkeitsbeweis darauf hin, da® es sinnvoll ist, die Stellen *x,, an denen die Ab- 
leitung genommen wird, von den Naherungen U, ZU unterscheiden, denn auch beim 


gewohnlichen Newtonschen Verfahren treten im zweiten Vergleichsverfahren (0.6) 


zwei verschiedene Folgen t, und Er Ommaltite 
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(2.2)|. Diese Funktion F verbindet das Verfahren (0.1) mit einem Vergleichs- 
verfahren entsprechender Gestalt 


Ogee = ene e) (7. == .03 Ae Me ele (0.5) 


Aus der Konvergenz des Vergleichsverfahrens wird auf die Konvergenz des 
Ausgangsverfahrens geschlossen, und auch die Fehlerabschatzung fiir (0.1) ge- 
schieht mit den GréBen des Vergleichsverfahrens. Eindeutigkeitsaussagen be- 
kommt man ahnlich mit Hilfe eines zweiten Vergleichsverfahrens 


UG T (a, a) (n = 0, 1,'2, .. .) (0.6) 
mit den gleichen é,. 
Als Abstand 6 und als Elemente a, &, t in (0.5) und (0.6) kénnen Zahlen oder 


auch allgemeinere GréBen (Elemente eines halbgeordneten linearen Raumes %, 
§ 1) verwendet werden. Bei dem Verfahren (0.3) benutzt man z.B. mit Vorteil 


den Abstand 6(u, v) = »[u—v], worin [wu] fiir w= ( ") den Vektor iW bedeutet 
[2 ad 

(|| || mumerische Norm). Man erhalt dann fiir Eigenwert und -vektor getrennte 

Fehlerabschatzungen. 

Fir das Newtonsche Verfahren (0.2) kann man eine Abschatzung der ge- 
forderten Form (2.2) herleiten (§ 3), falls die Differenz der Ableitungen G/,— Gj,) 
sich [nach (3.4¢), (3.24)] majorisieren 1a4Bt durch eine Differenz von Ableitungen 
eines im Raum % der Elemente o definierten Operators G. Die zweite Ableitung 
Gi wird dabei nicht bendtigt. Das zugehérige Vergleichsverfahren (0.5) ist dann 
gerade ein Newtonsches Verfahren von (0.2) entsprechender Form fiir die Glei- 
chung Go = 6. 

Als Beispiel wird (in § 4) das Verfahren (0.3) fiir lineare Operatoren A in 
einem Banachschen Raum behandelt (/ = lineares Funktional) und (in §5) dann 
auf Matrizen spezialisiert und an numerischen Beispielen erprobt. 

Ein anderes Verfahren der Gestalt (0.1) erhalt man, wenn man in (0.2) G’ 
durch Gj mit einem Naherungsoperator G, fiir G ersetzt (Zusammenhang mit der 
Stérungsrechnung). Dies Verfahren laBt sich ahnlich wie (0.2) behandeln (§ 6). 

Iterationsverfahren der Form w,.,—= Tu, werden als Spezialfall mit erfaBt, 
das Vergleichsverfahren (0.5) lautet dabei On+1=1 0,. Die iibliche Lipschitz- 
bedingung fiir T bei numerischem Abstand fiihrt auf ein Vergleichsverfahren 
mit linearer Funktion 7. Mit allgemeineren GréBen als Abstand und nicht- 
linearem Operator T lassen sich die Eigenschaften des Operators T oft besser 
beriicksichtigen. Man bekommt damit weitere Existenzaussagen und genauere 
Fehlerabschatzungen. Insbesondere kann man das vereinfachte Newtonsche 
Verfahren (x, =w) als Verfahren dieser einfacheren Form mit Tu=T(u, w) be- 
handeln. Es lat sich durch 7, =O charakterisieren, und eine solche Bedingung 
wie auch das Verschwinden weiterer Ableitungen TY) = --- =7(%)=O kann 
durch geeignete Vergleichsoperatoren T beriicksichtigt werden (siehe [25] und § 7). 

Funktionalanalytisch wurde das Newtonsche Verfahren zuerst von KANTORO- 
vitcH behandelt. Eine ausfiihrliche Darstellung findet man in [13]. Die ent- 
sprechenden Ergebnisse fiir das Newtonsche Verfahren bei gewdhnlichen 
Gleichungen bewies (z.T.) Ostrowski [2/|. Das Literaturverzeichnis enthalt 


156 JOHANN SCHRODER: 


eine Reihe von Arbeiten (ALTMAN [1], BarTLe [2], Feny6 [7, 8, 9], Hirasawa 
[10], Kantorovitcu [12, 13, 14, 15, 16], MERTVECOVA [17], Mysowsx1 [18, 19, 
20), SaALEKow [22], STEIN [26]) tiber das Newtonsche Verfahren oder Varianten 
im Banachschen Raum oder allgemeineren Raumen. Diese Verfahren werden 
darin in einer von der hiesigen Darstellung beziiglich Herleitung und Art der 
Ergebnisse abweichenden Form untersucht. Dieser Mitteilung inhalthch am 
nachsten stehen noch die Arbeiten [/5, 16], in denen KANTOROVITCH ein 
allgemeines Majorantenprinzip fiir Iterationsverfahren u,.,;=Tu, in eimem 
linearen Raum mit allgemeinerer Norm (anderer Definition) auf das vereinfachte 
Newtonsche Verfahren (mit x, —=%,)) unter Benutzung von Gi, anwendet und 
dann durch wiederholte Anwendung des Prinzips auch Aussagen tiber das ge- 
wohnliche Newtonsche Verfahren (x,—=wu,) bekommt. L.CoLiarz [3, 4] be- 
handelte das vereinfachte Newtonsche Verfahren bei nichtlinearen Randwert- 
aufgaben und bei algebraischen und transzendenten Gleichungen funktional- 
analytisch als Iterationsverfahren mit kontrahierendem Operator. 


§ 1. Bezeichnungen und Definitionen 


1.4. In §2 seien u,v,... die Elemente eines (vollstandigen) t-metrischen 
Raumes Mim Sinne der in [24] und [25] gegebenen Definition. Jedem Elemente- 
paar u, vCR sei also ein Element 6(u, v) eines halbgeordneten linearen Raumes ‘Vt 
als (verallgemeinerter) Abstand zugeordnet. Der Abstand sei definit und sym- 
metrisch und geniige der Dreiecksungleichung. 

In den §§3 bis 7 verwenden wir den Begriff des Yt-normierten Raumes N. 
Das sei ein lineares System von Elementen wu, v,... («+v,au erklart, « reell 
oder komplex), bei dem jedem Element uw ein Element y[w|]E® als (verallge- 
meinerte) Norm zugeordnet ist. Fiir die Norm gelte 


yiuJ>o (dh. 20,+0) fir w+9@ (0,9: Nullelemente von % bzw. %), 
v[u+v] <v[ul+r[o], vilaw)] =|a|v[w]. (1.1) 


Mit 6(u, v) =v[u—v] ist R dann ein Yt-metrischer Raum. 

Die dabei fiir die Elemente 0, ,... des (Abstands- oder Normenraumes) 
geforderten Voraussetzungen sind z. B. erfiillt, falls die Operationen o+o, xe 
(mit reellen Zahlen «) und @= 0 (ftir gewisse o, die dann_,,positiv genannt 
werden) mit den itblichen Rechenregeln erklart sind, % also linear und halb- 
geordnet ist, und auBerdem eine (numerische) Norm || || definiert ist, welche 
[auBer den iiblichen, (1.1) entsprechenden Forderungen| folgenden Bedingungen 
genugt: 


a) aus 0 Xo <o folgt ||o|| <||@| 


) 


b) aus. 0, 20 (4 = 4,:253)42)iund: lim Ilo —e,|| =0 folgt oo. (4.2) 


In den §§ 3, 6 und 7 wird vorausgesetzt, daB N ein Raum dieser spezielleren 
Art ist. Durch ||«||=||»[w]|| ist dann in 9 eine numerische Norm erklart. 
Begriffe wie Konvergenz, Stetigkeit, Vollstandigkeit, ... sind bei % und 8 im 
Sinne der Normen || @|| bzw. || || zu verstehen. Ist dann u (t)€ Reine fiiro<t<1 
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stetige Funktion des Parameters #, so ist auch y[w(f)] stetig, und falls % und N 
vollstandig sind, gilt 1 
»| fu ad < J(u (()] at. 
0 


In einer halbgeordneten linearen oe % 1aBt sich z.B. eine Norm definieren, 
falls ein Element CC®, ¢>0, existiert, fiir das es zu jedem o€% eine kleinste 


Zahl «, mit 
al % M1 =e ces oeas, 


(1.2): 
Beispiel. Mt: im Intervall [0,1] stetige Funktionen @ (x) ; 


C= C(x) >0, [lel] =e =max le, 


2. Die Ableitung eines Operators sei folgendermaBen definiert. St und & 
seien lineare, (numerisch) normierte Raume. G bedeute einen Operator, welcher 
Dc in S abbildet. Es existiere eine lineare Mannigfaltigkeit CCR mit den 
Eigenschaften : 

a) fiir vu, vED ist u— VEX, 

b) ist w€D, dann auch u+h fiir hE@ Fat Iles 8) (1.3) 
(Spezialfall: £= RK, D offen). 
Einen linearen (nicht notwendig beschrankten) Operator L, der & in © ab- 


bildet, nennen wir eine auf & definierte Ableitung Gi,) des Operators G an der 
Stelle u, falls 


in |G +h) —Gu—Lh|| =0 
gilt. 

Beispiel. R,S: Im Intervall [0,1] stetige Funktion w(x), || || = max | u(x)|. 

D: Funktionen w(x) € 2 mit stetiger zweiter Ableitung w’’(x), die die Rand- 
bedingungen (0) =a, w(1) =0 erfiillen und innerhalb eines Gebietes & der x, y- 
Ebene verlaufen, ©: 0S «<1, p(x) <y<y(x). 

@: Funktionen f(x) CR mit stetiger zweiter Ableitung h(x), die die Rand- 
bedingungen / (0) = A(1) = 0 erfiillen. 

Gu=—u'"'+f(x, u), f(x, y) erklart in G mit geeigneten Differenzierbarkeits- 
eigenschaften. 


Gy) h=—h” + j, (x, 4) h 


§ 2. Allgemeiner Existenz- und eee: oe 


2.1. t sei ein vollstandiger Jt-metrischer Raum. %, D und %, D seien Teil- 
mengen von }t bzw. Yt (=) mit § CD und ScD. wsei ein as Beneceement 
€%. (Im allgemeinen ist es zweckmabig, dafiir eine Naherung fiir die gesuchte 
Lésung zu verwenden, z.B. u.) 5 enthalte das Nullelement o und % alle Ele- 
mente x€D mit d(x, w)€ %. Liegt 9 Soin % bzw. 9, dann sei auch jedes positive 
o<o aus % bzw. ®. 


2 Fiir Vektoren definiert so HOUSEHOLDER [11] eine Norm, der allgemeine Fall 
wird in [25] benutzt. 
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T (u, x) bedeute eine Abbildung der Elemente ED, «EF in KR. Wir unter- 
suchen Iterationsverfahren der Art (0.1) 


Uy oye L(y hq) Oy, 2p ee) UyED, MEV. (2.1) 
Es gelte eine Abschatzung der Form 


6(T(u, x), T(v, y)) <F [6(u, »), 6(u, wv), 6(v, w); 0(x, y), 6(x, w), 6(y, w); 
6 (u, x), 0(v, vy); 6(T(v, y), 4), 6(T(, ¥), v) | (2.2) 
fiir alle Elemente wu, v, x, y, fiir welche beide Seiten der Ungleichung erklart sind. 


Dabei bedeute F[o!, 0?,...,0!]=F[e'] eine Abbildung in % der Elemente 
{ol, 02, .-., °} = {o"} mit PEM, of Zoli =1, 2,..., 10) und* 


GG? UOMO UOT Ore: (2.3a) 
ot + 08, CS. (2.3 b) 
F sei positiv und monoton: 
OSFilocl=Fle) fir cae se. BG=-detee 10) (2.4) 
F hange im Falle @4=0 nicht von ¢° ab. 


Es existiere eine Abbildung T(o, &) der Elementepaare 0€9, ECS in N 
derart daB 
Le, 5) Ope NO Os ote pe 
oo, C39; 16.) Oe, | (2.5) 
mit allen 9, 0, €,7 gilt, fiir welche die darin vorkommenden GroBen erklart sind. 


Fiir den Konvergenzbeweis verwenden wir ein Vergleichsverfahren in ‘Vt: 
mea le Sa) (nm = 0, 1, 2, aes 20€, EES (2.6) 


und fiir den Eindeutigkeitsbeweis ein zweites Verfahren 


Trt1=1 (tm, &,) (w=0, 14, 2, ...) (2a) 
mit den gleichen &,,. 
Dabei gelte 
0 (Up, @) SQ, (2.8) 
0 (uy, %) SQ,— Op, Ah. T (Qo, Fo) 20(%4, %) + Qo 1. (2.9) 


2.2. Mit Hilfe der Abbildung F werden Beziehungen zwischen Elementen € Xt 
auf Beziehungen zwischen Elementen €% tibertragen, wird eine Verbindung 
zwischen dem zu untersuchenden Verfahren (2.1) und den Vergleichsverfahren 
(2.6), (2.7) hergestellt. Gegeniiber (2.2) sind in (2.5) die Elemente u, v, w, x, y 
durch @, o, 0, €, 7, der Operator T durch ig: der Abstand zweier Elemente € it 
durch die Differenz der entsprechenden Elemente €% ersetzt. 


8 Im allgemeinen kann man D=X wahlen, die Forderung (2.3a) entfallt dann. 
* Diese Forderung stellt bei den in den §§ 3, 6 und 7 behandelten Verfahren nur 
eine immer leicht zu erfiillende Normierungsbedingung dar. 
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Im Spezialfall des Iterationsverfahrens u,.,—= Tu, vereinfacht sich vieles. 
Tu ist fiir T(u, x) und To fiir T(g, €) zu schreiben, und es werden dann nur 
Abschatzungen (2.2) betrachtet, in denen F nur von g!, 9? und 9? abhangt. Ab- 
schatzungen dieser einfacheren Gestalt lassen sich beim vereinfachten Newton- 
schen Verfahren herleiten ([25] und §7). Auch die Ungleichung 6(Tu, Tv) < 
K6(u,v) (K linearer Operator [24] oder Lipschitzkonstante [27], [5] bei numeri- 
schem Abstand 6(u,v)) hat eine solche Form. (2.5) und (2.9) sind dabei erfiillt 


fiir To =8(uy, UM) +00—Keo+ Ke. 


2.3. Wir betrachten nur Vergleichsverfahren mit monotoner Folge &,, die 
durch die Folge @,, majorisiert wird: 
EL i SOAS. OF Gish eee == 0) (2.10) 


Soweit die 9, nach (2.6) mit Elementen &; erklart sind, welche den Ungleichungen 
(2.10) geniigen, gilt: 


Ont =O (i Olean Oay=—.0),3 (2.11) 
Denn fiir »=0,1 ist diese Ungleichung nach (2.10) bzw. (2.9) richtig, und gilt 
sie bis =f (= 1), so folgt mit (2.5) und (2.4): 
Cp+1— Op = TQ», ame T(Qp-1, Ep—1) SF [Op — Op—1; Op» Op—13 
Ep — Sp—a, Sp Sp—13 Op — Ep» Op—1 — Sp—1 Op — Op» Op — Op-i] 2O- 
E,41=Op41 Oder &,.;= 0, sind daher z.B. Elemente, mit denen die Unglei- 


chungen (2.10) fiir» —£-+1 gelten, falls sie bis n=# erfiillt waren. Allgemeiner 
gentigt die Teilfolge 


CFs Oy aC NSS ee oy ary (G0 OF AGD, So) ee 2-12) 
(speziell €, =0, 09 oder e,,) den Ungleichungen (2.10), soweit die Iteration damit 
durchfihrbar ist. 


2.4. Auch an die Folge x, in (2.1) werden [(2.10) entsprechende] zusatzliche 
Forderungen gestellt. x, sei jeweils so gewahlt, dal 


O(U), FESO, Sa (SOAS 20 eee (2.13) 
Oe) ee ee == Oat ee 0) (2.13b) 

gilt. (2.13b) ergibt: 
ONG Pay 59) bp Sais (0, = ALO ARS m= 0) (2.14) 


und insbesondere 6(%,,,@)<&,,, so daB jedes in D liegende x,, aus & ist, falls 


é 10 % liegt. 
Soweit die Folgen u,,, 0, mit Elementen x; bzw. &; erklart sind, welche den 
Ungleichungen (2.10) bzw. (2.13) gentigen, gilt 


O (gy Mp) SS Oy Cn (Mi, W= A, 0/1 62,4. ns ul = w). -(2.15) 


Fiir m, n= —1, 0, 1 ergibt sich dies mit 


6 (Up, U4) SO (Ms, Xo) + O(%, X41) So Ey + 5 —€4 = 0p — 0-1 
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und (2.9). Ist (2.15) bis m=p (21) richtig, so folgt mit (2.13) und (2.14) 


O (psy, Mp) = O(T (up, %p), T (mp1, pe 
<F[O (uy, Up—1), 0 (Up, @), O(Mp—1, &); O (Xp, Xp—s), O(%p, W), O(%p-1, W); 
O (Up, Xp), O(Up—1, Xp—1); 0 (Mp, Up), 6 (Up, Up—x)] 
ef (Gp — Ope p epi asp = Epa Sp. Spa 
Cp ae Cp—1 Sper: Op — Op» Op — Op—1] 
= Le; é,) ae Tor E,—1) = Op+1 Op 
und damit auch (2.15) fir m=p-+1. 

Unter welchen Voraussetzungen kann man x,€D nach Ermittelung der vor- 
hergehenden a), “,,..., Up€D (P= 0) so wahlen, daB die Ungleichungen (2.13) 
auch fiir n= gelten, wenn sie (im Falle >0) bis n=p—1 erfiillt sind? 

Im Falle (2.12) &,=@,, gentigen wegen (2.15) z.B. die Elemente x, —w,, 
den Ungleichungen (2.13), soweit die darin vorkommenden Gr6fen erklart sind. 
Praktisch bedeutet das: Wenn man das Verfahren w,.,—= T(u,, u,,,) untersuchen 
will, wird man @,4,;= INO 0,,) als Vergleichsverfahren wahlen. Damit werden 
also z. B. Newtonsche Verfahren (0.2) behandelt, bei denen man in die Ableitung 
immer das gleiche Element, etwa x, =, (vereinfachtes Verfahren), oder das 
zuletzt berechnete x, =u, (gewOhnliches Verfahren) einsetzt oder aber all- 
gemeiner als Elemente x, die Naherungen uw; verwendet und sie der Folge x, =, 
entsprechend wechselt. 

Bei allgemeinerer Folge é,€%} laBt sich immer ein solch gesuchtes x, ED 


angeben, wenn 9 ein Banachscher Raum, D konvex und 6(u, v) =||#—v|| ist: 
Ep — Sp—1 Op — Sp 
u,+_—>—" x,_,, falls ee 
Xp = Cree 4 ae Gre ops eat Op Sp 


Up » falls @,=&,_). 

Man beachte dabei, daB wegen (2.13) (us, %p_1) <0, —&p-_, Bilt. 

Ein Beispiel mit einer anderen Folge x, wird in § 4 behandelt. 

2.5. Sind die unter 2.14 genannten Voraussetzungen erfiillt, so gilt der folgende 

Satz 1. a) Es gebe durch die Iterationsvorschrift (2.6) verkniipfte Elemente 
0,€®, & €% (1 =0, 1, 2,...), welche den Ungleichungen (2.10) gentigen, und die 
Folge 0, (n=1, 2,3,...) besitze ein Grenzelement oNa= lint T(o*, e) ED. D ent- 
halte alle Elemente uC KR mit 6 cae 

0(u, 4) < 0* — Q. (2.16) 


Dann ergibt die Iterationsvorschrift (2.1) mit Elementen x,€®D, welche die 
Forderungen (2.13) erfiillen’, eine Folge u,CD (n=1, 2,3,...), die gegen ein 


°® braucht nur die Elemente u¢T(%, D), welche der Bedingung (2.16) geniigen, 
und deren Haufungselemente zu enthalten. Eine noch schwichere Forderung, welche 
fiir die unmittelbare Anwendung auf Differentialgleichungen geeignet ist, wird in 
[25] (Zusatz 1a) benutzt. 

ie Diese Forderung sichert, da® die Iteration (2.1) unbeschrankt durchfiihrbar ist. 

Sie ist immer erfillt, wenn man D=t setzen kann. Fiir das Iterationsverfahren 
u,+1—= Tu, ist das bei Verwendung eines linearen Vergleichsverfahrens im allgemeinen 
nur moglich und zweckmaBig, wenn T selbst linear ist, bei nichtlinearen Vergleichs- 
verfahren dagegen in weit gréBerem Umfange. 

’ Also unter der Voraussetzung, daB solche %, existieren (vgl. 2.4). 
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Element u* = lim T(u*, x,)€D konvergiert, und es gilt die Fehlerabschatzung 
n—->Co 


0(u*, u,)<0*—o, (n=0,1,2,...). (2:47) 


b) Ist durch (2.7) fiir irgendein Element t,€® mit 09XT, eine Folge t,€D 


(n= 1, 2,3, ...) definiert, die ebenfalls gegen o* konvergiert, und gibt es ein Element 
uC, fiir das8 
aoe O(U, Up) ST — Oo (2.18) 


WaT (un, %,)\ far W= 0,1; 2) 35%. (2.19) 


und 
gilt, so folgt u=u*. 


Beweis. a) Wegen 0* — 0, = 0 geniigt w =u, der Ungleichung (2.16), liegt also 
in D. Sind alle u; aus D mit 7<p (P=1), so folgt aus (2.15) mit m=1, daB 
auch w=u,,, und damit jedes u; der Ungleichung (2.16) geniigt, also €®D ist. 
Da die Folge 0, konvergiert und # vollstandig ist, besitzt die Folge u, wegen 
(2.15) ein Grenzelement u*€%. Dieses erfiillt infolge der Stetigkeit des Ab- 
standes 6 ebenfalls die Ungleichung (2.16). Die Fehlerabschatzung (2.17) ergibt 
sich durch Grenziibergang m-—> oo in (2.15). 


Mit (2.13) und (2.47) schatzt man fiir » = 0,1, 2,... ab®: 
O(T(u*, x,), u*) <6 (T(u*, x,),T (Uy, %)) + b(t sy, ¥*) 
<F[d(u*, u,), d(u*, w), O(u,, w); 0, O(X,, w), O(x,, w); 
O(u*, x,), O (Uys Xn) ®, O (Uys, Un) | +0 (4,41, #*) 
SP lQ* = Ons O*, Ons Or Fn» Eni OF — Ens On — Ens *s Onda — On] + O* — On ta 
Sel Ome) On. te) 1 0" Ont Ot pem) toh eat On) t 
Der letzte Ausdruck konvergiert gegen 0, also ist w* = lim T(u*, x,). 


n—>Co 


b) Es gebe ein Element u€9, fiir das (2.18) und (2.19) gelten, und die Folge 
T,€®D konvergiere gegen g*. Dann ist also 


0 (u, Uy) = Tr — On (2.20) 
fiir »=0 richtig. Gilt diese Ungleichung bis » =f (20), so schatzt man ab: 


6(u, w) <d(u, Uy) + 6 (Up, W) ST, — Op + Op, 
6 (U, Xp) KO (u, Up) + O (Uy, %p) STs — Op + Op — Sp 
und damit ® 
6 (6, Up 43) = 6(T(u, xp), T(up, %p)) 
<F[d(u, up), 6(u, w), O(uy, w); 0, A(X, w), O(%>, w); 
O(u, Xp), O(Up, Xp) *, O(Mp +1, Up) | 
Ons Tpri Ops On Sp Spree SpoOpe op ate Open Cel 
<T (ty, is T (ep, &5) = Tp+1— Opti: 
8 Die Menge der w€% mit J(u, Up») <T) — Om beiirgendeinem m (m= 1, 2, 3, ..., ©; 
Ups = U*, Oo = 0*) ist in der Menge (2.18) enthalten: 
6 (U, Uy) SO (U, Up) +4 (Um, Uo) ST — Om + Om — Co = To — Qo- 
9 Wegen 04=o0 hangt F von @° nicht ab. 
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(2.20) gilt dann also auch fiir n»=p-+1 und damit fiir alle n=0, 1,05. .eebe 
beide Folgen 0, und t, gegen @* streben, ist nach (2.20) «= limu,=—w*. 


Zusatz 1. Guibt es ein Element te, mit dem 
Oot und teT (rt, &) fiir jedes EC 5 (2.21) 


ist, so laBt sich die Iteration (2.6) mit jeder der Bedingung (2.10) gentigenden 
Folge &,€% unbeschriinkt durchfiihren, und fiir die Folge 0, gilt daber 


Ont ag Pere ran be (2.22) 


Dieser Zusatz laBt sich oft zur bequemen Fehlerabschatzung verwenden. 
Ferner kann man aus der Monotonie und Beschranktheit der Folge @, vielfach 
auf deren Konvergenz schlieBen. In manchen Fallen gilt das erste Konvergenz- 
kriterium in MN, in anderen muB man dabei benutzen, daB die Folge 0, €T (B, 3) 
ist?’ 

Beweis des Zusatzes: Die Elemente g, seien mit (2.10) gentigenden ECR 
(j=0,1,2,...,P—1) bis n=p (]O) erklart und <tr. [Fir p=0 ist dies der 
Fall, (2.10) entfallt dabei.] Dann ist auch 0p€D und mit jedem E,€%, fiir das 
E, ,X<f,< 0, gilt (also z.B. &,=&,_,), nach (2.6) ein g,,, definiert. Fiir dies 
Element gilt ®: 


.— Ope ie $,) hi T (dp, 5») 
ZE[t— Op, T Op; 0 Fp, ps T — Sp» Op — Eps * Opti — Op] =O. 


Damit ist die Behauptung bewiesen. 


§ 3. Anwendung auf das Newtonsche Verfahren 


3.4. Das Newtonsche Verfahren (0.2) soll nun mit Hilfe der Ergebnisse von 
§2 untersucht werden. Ht, S und % seien Banachraume. Dabei sei Kt Yt-normiert 
in dem in § 1 beschriebenen Sinne und die numerische Norm in ®t durch || «|| = 
||» [w]|| definiert. 


Zu lésen sei eine Gleichung 
Gu=® (= Nullelement von %) (3.1) 


mit einem Operator G, der eine konvexe Teilmenge DCR in S abbildet. Es 
existiere eine abgeschlossene™ lineare Mannigfaltigkeit @CR, die mit D durch 
die Eigenschaften (1.3) verkniipft ist, und G besitze fiir jedes w€D eine auf & 
definierte Ableitung G;,). 


Es gebe ein Element w€® derart, daB der Operator [Gi] 4 existiert und 
auf ganz © definiert und beschrankt ist. Der Operator 


Hu = [Gw]1Gu (3.2) 


MY So wird z.B. in [6] vorgegangen. 
*\ Die Abgeschlossenheit von 2 wird bei (3.7) bendtigt. Man kann auf sie ver- 


zichten, wenn man verlangt, daB &¢(A(, — Hj) & gilt. Dabei bedeutet (Hx) — Hie) 
die stetige Fortsetzung des Operators (Hi, — Ay) auf die Menge &, welche aus & 
durch Zufiigen der Haufungselemente entsteht. 
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bildet dann D in & ab und besitzt eine auf & definierte Ableitung 
Hw = [Gwe] * Cw, (3.3) 
deren Wertebereich ebenfalls in @ liegt. 

Fir diesen Operator H soll nun eine Majorante existieren im Sinne der fol- 
genden Forderungen. Es gebe einen Operator H, der eine konvexe, offene?? Teil- 
menge DCX in Y¥ abbildet, die mit @=o0 auch jedes positive Element o< 0 
enthalt (mindestens ein g > 0 soll in D liegen). H besitze fiir alle 0 €® eine auf 
N definierte, in jeder Teilmenge o<o<w (w€9) gleichmaBig beschrankte Ab- 
leitung H(,), die stetig von 9 abhangt. Es gelte ferner 

a) Ho=o, Hj,=O (Nulloperator), 

b) 0S (Aie+0) Hig) tS (Ae) — Hey) t 

Hite oaoe el. ono CD, i410: 
insbesondere also nek, \t=o fiir o<e ED), += 0, 
c) » [(H; y— Hy) Ns (i v(u—v]+y»r[v—w]) —Hiy- =i) y[h] 38 
fiir u,v€D, hEL, v[u—v]+r[v—w]€9. (3.4) 


% bedeute eine (nicht leere) Menge von Elementen &>o0 aus 9, fiir welche 
der Operator 


ye = [E -Hy)? (3.5) 


auf 9% erklart und beschrankt ist und die mit @2o auch jedes positive o< 0 
enthalt (E = Einheitsoperator). 


3.2. Aus (3.4c) folgt 
(Hig — Hw) ASH oxy v4] fir «ED, HEL, y[x—w]E§. (3.6) 
Fiir «€ D, v[x—w]€ % ist daher der Operator 


(oe) 


= |E —[G (Gig — Gy) (A =[E— (By He) = 2D (BiB) _G-7) 


j=1 
auf ganz & definiert und Yt-beschrankt, d.h. beschrankt im Sinne der Norm y: 


»[R (2) hl 1s Del ie ns > ote—ep) 9A 
is H, es w))] *y[h] y (3.8) 
Damit besitzt G(x) fiir jedes solche x die auf ganz © erklarte Inverse 


[Gi] t= R(x) [Geo]. (3.9) 


12 Die Forderung, daB Dd offen sein soll, und auch die Voraussetzungen tiber D 


lassen sich mildern. 
13 Geniigt z.B. im Falle y[w]=||u|| der Operator G statt H einer solchen Un- 


gleichung, so kann man auch fiir H eine entsprechende herleiten, indem man abschatzt: 
|| ly = Ae) A] S [GIN (Gly — Ge) AIT 

Jedoch ist es zweckmaBig, (3.4c) direkt von H zu verlangen, wie man am Beispiel 

von § 4 erkennt. 


Nie 
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Mit % sei nun irgendeine Teilmenge von D bezeichnet, fiir deren Elemente x 
[Gi] existiert und auf © definiert ist und die alle x©D mit y[x—w]E@ 5 
enthalt. T(u, x) bedeute dann die durch 


T(u, x) =u — [Gy] 1Gu (3.10) 


definierte Abbildung der Elemente w€D, «€% in Rt. Das damit gebildete Itera- 
tionsverfahren (0.1) ist mit dem Newtonschen Verfahren (0.2) fiir die Gleichung 
(3.4) identisch. 

Die Gleichung Gu* =9@ ist jeder der Bedingungen (0.4) aquivalent, wenn 
y[x,—w] eine durch ein Element w€® beschrankte Folge €§ ist (in Satz 1 
und 2 ist »[x,—w]<&,<o*=). Wir brauchen nur zu zeigen, daB Gu*=@ 
aus u* = lim T(u*, x,) folgt, d. h. aus lim hy, = Toit 


= [Gi] 1G U* eR (ay) [Gi Oe 
Nach (3.7) ist 
[Gini Gut = Rs, = (Hw) ¥ FH ,)) | hy, 
und damit wegen (3.6) und (3.4b) 
v[[Gw] 4G u*] <[E + Hite, —wy] ¥ Pn] S[E + Diy] oa], 
also [G(w)] + Gu* =@ und folglich Gu* = 9. 


3.3. Wir beweisen nun, daB fiir die so definierte Abbildung 7T(u, x) mit dem 
Abstand 6(u, v) =yv[u—v] eine Abschatzung (2.2) gilt, welche allen in 2.1. ge- 
nannten Voraussetzungen geniigt. Da wir die Existenz einer zweiten Ableitung 
G’’ nicht vorausgesetzt haben und méglichst allgemeine Vergleichsprobleme zu- 


lassen wollen, werden im Beweis einige sonst vermeidbare Grenziiberginge er- 
forderlich. 


Mit (3.10) und (3.9) erhalt man fiir u,v, x, yED; »[x— vw], v[y—w] ES 
T(u, x) — Tv, y) = [Gig] {— Gu —Gi, (TO, ») —u)} 
= [Gin {— Gu — Gig (TO, 9) —) + G04 Gy (LO, 9)—0)} 
= R(x) if, (3.11) 
f=—Au+dHv Hy (z—u) + Hy (z —v), a=! (Uy). 


f laBt sich darstellen als f= lim f(s), 


s>-+-0 


f(s) =— Hut Hv +—| H (we +.s(2 u)) +Hx +H (y+ s(z v)) Hy| 


1 
es If 1 , / 
Ta v) dt 4 ={t El eaetaa)) (x—y—s(u- v)) + Hq (x—y)] dt 
0 
1 


aa 1 , , 
~~ ft |] (u—v) dit— ; 5 | Herse-m— He] (x—y—s(u—uv)) dt, CG AQ} 


0 
@=tba(1—t)v, X= tes (eae 
Hy hangt infolge (3.4c) stetig von u ab. 
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Voriibergehend wollen wir mit 0’ die folgenden speziellen Werte bezeichnen: 
g=r[u—v], o&=r[u—w], o =v[v—w]; of =r[x—y], S=r[x—w], 
=[y—w]; Pv =»[u—x], & =r[v—y]; oP =v[z—4], o=v[z—V]. (3-13) 
Damit gelten dann fiir O0<t<1 die Abschatzungen 
»[@—#] =o[t(w—2) + 1-1) Wy] <te’ +1’, 
UA Oe? Ok) ty — wie arte oo (ar at! of ye ol 8) 0%, 
y[z—%] = v[t(z—u) + (4 —2) (ze —v)] St0® + (1-2) 0”. 


Mit diesen Ungleichungen und (3.4) erhalt man von (3.12) ausgehend (fiir s > 0) 
rif S06. 0), (3.14) 
9(50) = f [Eee rorey+a— etter) — Hiren] oat + 
0 


ow 


1 
1 ry! , 
= ~f [A (s¢or+—n o) trott+e’) — Lerer+er] (Q* + SQ!) at 
0 


soweit dieser Ausdruck g(s, 0‘) erklart ist. 


Wir fassen die 9’ nun wieder als unabhangige Variable auf. g(s, o') ist bei 
gentigend kleinem s > 0 definiert fiir 


Geen Oe 0 0 no Cone Otc e. (3-15) 
Wegen der Eigenschaften des Operators H gilt in diesem Definitionsbereich 
0<(s,0')< G(s, o') fiir o<o’<o' und geniigend kleines s. (3.16) 
@(s, o') 1aBt sich darstellen als 
| 9(s, o') = Pile!) + Pals, @') + Gals, @') 
mit 


71 (0') = f Aiororerta—n (e°+e°)) @ dt, 
0 
hy 
Ss 


f Bou (¢o°-+(1—#) 0!) +¢0!+ 0°) (s (0° as o"°) +0*) — Heeate o*| at 


=— [AG e+ of + 0) —H(so" + 0°) —H(ot +o) +Ho'), 
1 
P3 (s, o') cae if [7s (tor+ 1 °°) +t04+ 0°) (o" ce 0° ate o”°) — Hie o%) 20" dt. 
0 


Die beiden letzten Ausdriicke streben fiir s-+-+0 gegen die Grenzelemente 
92(o') = H (e.¢%) 0° — Hi) om 
bzw. 


1 
G0) =) Hearn 6 — © 4 e@.) at 
0 
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Insgesamt hat man also 
Jim 96s, @) = pile) + v2le) + P0(@) = P(e). 
p(o') ist erklart, falls (3.15) gilt. Aus (3.16) folgt in diesem Bereich 
o<g(o')<(o') fir oso'<o’. (3.17) 


Im Falle 9=o0 hangt @ nicht von 9 (und auch nicht von 9°) ab; die 0° (und @!%) 
enthaltenden Bestandteile heben sich weg. 
Durch Grenziibergang s->-++0 erhalt man aus (3.14) mit den Werten (3.13) 


vif}<(o'). (3.18) 


Der fiir o<& es erklarte Operator [E =i @)t ist positiv, wie aus der Dar- 
stellung (3.5) folgt. Ferner nimmt er wegen (3.4b) mit € monoton zu. 
Insgesamt erhalt man so mit (3.8) und (3.18) von (3.11) ausgehend 


v[T(u, x) —Tv, )) = R(®) A <[E Ain) 9 SE Hin) (0), 


soweit die durch (3.13) gegebenen Werte 0‘ den Bedingungen (3.15) gentigen und 
auch 0°€% ist. Das ist eine Abschatzung der Form (2.2) mit der Funktion 


F [g'] = [E — Hin] (0). 


Wegen der bereits bewiesenen Eigenschaften des Operators [E —H (] 2 und der 
Funktion (o') braucht nur noch gezeigt zu werden, daB (2.5) mit einer ge- 
eigneten Abbildung 7 gilt, welche auch (2.9) erfiillt. 

Setzt man in F die Werte 


g=o-0; of=t—n; ge, 8=n; oe —o—k, CO iG eG ne ra wee 
ein, so lassen sich die auftretenden Integrale berechnen, und man erhalt 
F[o'] =[E —H)" (Ho —Ho + Hig 6 — 0) —Hiy€—9)). (3.19) 


Mit dem fiir 9€D, €€% bei beliebigem Element 0, =>”[u, — Up] erklarten Operator 


T(e,4=e+(E—Ai,)7[(E —H gy) 0 


Go H Q, Or H@] (3.20) 
sind die Forderungen (2.5) und (2.9) erfiillt, wovon man sich tiberzeugt, indem 
man ¢=T (a, 7) in (3.19) einsetzt und T (Qo, &) bildet. 

3.4. Es sei G der durch 


~ 


Ge = (E —H,) o,+ e—Ho—e+He (3.21) 


in D definierte Operator, insbesondere ist 


Ge=o,—o0+Ho fir f= e)=0. 


Auch fiir G statt A gilt die Forderung (3.4c). Mit diesem Operator nehmen die 
Vergleichsverfahren (2.6) und (2.7) die Gestalt (3.22) bzw. (3.26) an. — Ahnlich 
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wie unter 3.2 zeigt man, daB die Bedingung o* = lim T(o*, é,) der Gleichung 


G o* =o aquivalent ist, wenn &, eine durch ein Element w €9® beschrankte Folge 
€® ist (in Satz 1 und 2 ist &, <o* =a). — Ist r>0 ein Element as) mit Gt<o, 
so ergibt sich fiir jedes ECS: 


tet+[(E—H]4Grt=T (cr, 8). 


Die damit unter den in 3.4 genannten Voraussetzungen aus Satz 1 und 
Zusatz 1 fiir das Newtonsche Verfahren folgenden Ergebnisse lauten: 


Satz 2. a) Die Iteration 
On+1 = On — [Gey] * G0, (n = 0, A 2; aa 0<0,€D (3.22) 

set mit einer der Forderung 
nee an TO Oe Fe) (3.23) 


gentigenden Folge &, € § unbeschréinkt durchfiihrbar, die Folge 0, konvergiere gegen 
eine Losung ox ED der Gleichung G o* =0. D enthalte alle Elemente uC R mit 4 


y[u — u,]< o* — 9,. 
Dann laBt sich nach der Vorschrift 
Un = Uy — [Gilt Gu, (w= 0,1,2,...) fiir ugED mit e(uy,w)< ey (3.24) 
eine Folge u,€D berechnen, falls man x, €D jeweils so bestimmen kann, daB 
Dit ty eg oye! Kye hy bys gey (WO, 15 25.3 Kay =) 9-25) 


gilt, Diese Folge u, konvergiert gegen eine Losung u*€D der Gleichung Gu=@Q, 
und es gilt die Fehlerabschatzung 


y [u* — u,,|S 0* — (== OAR Oe 


b) Ist fiir ein T)ED mit 0y< To. durch 


Ge TENG (Gt, © (Uh Oy tye) (3.26) 
eine Folge t,€® definiert, die ebenfalls gegen 0* konvergiert, — diese Voraus- 
setzung ist fiir t= 0* erfiillt —, so gibt es aufer hichstens u=u* keine (in D 
liegende) Lésung der Gleichung Gu=Q, welche der Ungleichung 

v[u — U9] ST) — Oo (3.27) 


geniigt 16, 
Zusatz 2. Gibt es ein Element t CD mit 
= Gt<o, 
so ist die Iteration (3.22) mit jeder (3.23) gentigenden Folge E,€% unbeschrinkt 
durchfiihrbar, und es gilt oe mabye 


Dabei liegt jede (3.23) geniigende Folge &,, in §, falls tC § ist. 


149 braucht nur die Menge der Elemente u=v— [Gi] Gu (vED, x€%), welche 
vy [uw — u,) <e*— Q, erfiillen, und deren Haufungselemente zu enthalten. Vel. FuBnote ®. 

15 Dazu siehe 2.4. ae 

16 Auch u = u* braucht jedoch (3.27) nicht zu erfiillen. Siehe das Beispiel von § 4. 
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3.5. Die Abb. 1 bis 3 veranschaulichen verschiedene Falle fiir die Vergleichs- 
verfahren bei numerischem Abstand!’, — Die Funktion Go sei auf der ganzen 
Menge = der reellen Zahlen definiert. « bedeute die gréBte Zahl (eventuell 
auch «=- oo), fiir die im Intervall O<e<« 


Abb. 1 


Abb, 2 


Abb, 3 


Abb, 1—3. Vergleichsverfahren bei numerischem Abstand 


or =a =9" 1, (n=912...) 


(Zh seat, Ci iaG 


gilt. Dieses Intervall [0, «) 
konnen wir als Teilmenge § 
verwenden. Wir nehmen weiter 
o, >0 an. Damit ist fir 


0<& = ES 
Go.> 0 und Gon) <0 


Besitzt dann die Gleichung 
Go =0 eine Lésung e* >0 — 
wie in den drei skizzierten 
Fallen — so konvergiert jede 
nach (3.22) gebildete Folge 
gegen diese. Als &, kann dabei 
jede positive, monotone Folge 
verwendet werden, bei der &, 
links“ von o, legt (in den 
Bildern wste3— e-). 

Ist ty irgendeine Zahl mit 
ere < oder aber. 7) ==6" is0 
konvergiert die (mit den glei- 
chen &,) nach (3.26) gebildete 
Folge t, ebenfalls gegen o*. 
Besitzt also z.B. die Funktion 
Go eine zweite Nullstelle ort 
(Abb. 1), so kann man jedes Ty, 
mit 0*< t)<o** verwenden 
und der Eindeutigkeitsaussage 
von Satz 2 die Form geben: fiir 
jede eventuell in D existierende 
Lésung u** = u* der Gleichung 
Gu=9Q gilt 


[| w** — ug|| 2 Q** — Qo. 


§ 4. Berechnung des groBten Eigenwertes eines linearen Operators 


4.1. Gegeben sei eine Eigenwertaufgabe 


(4.1) 


fiir einen linearen Operator A, der einen Banachraum %’ in sich abbildet. Das 


™ Die vertikale und horizontale Koordinatenachse haben in den Abb. 1 bis 3 
jeweils verschiedenen Mafstab, es ist G’ 


2) 
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gesuchte Eigenelement werde mit einem auf t’ definierten linearen Funktional / 


durch 
lp=1 (4.2) 


,normiert’’. Die Ergebnisse des §3 sollen auf das Iterationsverfahren 


en tA Pn> 1 Qn = 1 (n = Onl Kae .) (4.3) 


angewendet werden. Im allgemeinen konvergiert die Folge 1/u,, gegen den (dem 
Betrage nach) groBten Eigenwert von A und die Folge », gegen ein zugehdriges 
Eigenelement. 


4.2. Wir schreiben im folgenden meist wu’ fiir das Element q, w’’ fiir die Zahl 


1/A. St sei die Menge der Kompositionen u = a GS =’. G bedeute den auf 
der Menge D der w€ # mit J wv’ =1 durch e 


Gu=wu—u"'Aw (4.4) 


definierten Operator. Die Gleichung Gu=@ (w€D) ist dann (fiir A=0) der 
Aufgabe (4.1), (4.2) aquivalent. 
Jedem u€ Rt wird der (zweidimensionale) Vektor 


roa =| 
ja" 
als verallgemeinerte Norm zugeordnet. Als Norm der Elemente 9 = = EN 
@ 
|, |e" |, |"). 
G besitzt fiir jedes w€D die auf der abgeschlossenen Menge & der w€ Ht mit 
lu’ = 0 definierte Ableitung 


Gh =f! — hi Aw — wu" Al, 


, 


). Damit wird || || = ||» [2] || = max (|| 


wahlen wir || || = Max ( 


Sei % die Menge der Elemente x€D mit x’’=0,7Ax’=+-0, so existiert fiir +E 


patty ay 
[Gin] = es ae lt 
et, 
und es gilt damit 
wu’ rea A x’ — Byy(u" u’) 
[Ginl oo = 1 heen) 
Php ra = V(x) (U u ) 
Dabei ist 
y, ; 1 ze h 
Bap =Ap—S4 Ax, MP = Fy AY fiir men. (4.5) 
Iaube Ze a wird 
AW Au’ 
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Das Newtonsche Verfahren (3.24) mit dem Operator G aus (4.4) und den Elementen 


= (%), =| (w= 04,21), =O (4.6) 
uU 0 


ist daher dem zu untersuchenden Verfahren (4.3) aquivalent. 
4.3. Es sei nun w= i) = 4 ein festes Element €3, B= By) und r=7%q). 
Es gebe Konstanten a, b, mit denen 
|Bells@|le|| und |rg|sollel| fir lp~=o (4.7) 
gilt. Der Operator H = [G;,| 4G besitzt die Ableitung 
; h’ Sy 8; h”’ ae oa hy 
Hw) h Te ( ? ff ) : 
r(h' u' —u" h’) 
Man schatzt ab: 


rt ayy ap ee Oh ea 


|r (h''(w' yest v’) ae (u”’ — #2’) h’) 


Ny 


’ , / ” tr , a , ry! 
= (|r | || = | ie | Yo h 1) (3) = (1, (.—si+»{0—o}) —4 b po_w)) 7A 


mit dem in = erklarten einfachen Operator 


dessen Ableitung 
ry! ee FB eee a 
Ho =(0'"'9’ +@ “1(5) 


ist. H besitzt alle in 3.1 geforderten Eigenschaften. 
Man berechnet 


Ty’ j , RIN sot ie eT a . 
[Hilo = (60' + 40e"))* (0 0' + 0 (5) (j =4, 2,3,...). 
Als % kénnen wir daher die Menge der Elemente € = ( é mit 
0 


' f0, O50 24 
wahlen. Fiir €€% ist dann 


[E—Hy)'o=04+ > [Alico =o4 : > o"' & al 
ver 1— bE b 


Damit sind fiir das hier behandelte Beispiel alle in 3.14 vorkommenden 
Begriffe erklart. Die dort genannten Voraussetzungen sind erfiillt. 


4.4, Mit eine ae 
i Ug — = 
c= f= o—w | ee Ga A 
0 | my | 


a I Dest» , wu a 
Ge=o,+ %—e+(0'e ~eai)(5). 


wird 
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Unter der Voraussetzung ; 
und 


= On eiaus 170 ==.0 


1—b 99 — Yb(o1 + 00) — aor 
\ Go — YP(o1 + @0) Yani | => falls ab =0 (4.9) 


besitzt die Gleichung Go=o zwei (eventuell zusammenfallende) positive Lésungen 


=a tert (5), o** =0,+ O+y** el (4.10) 


wobei y* und y**(=> 0) bzw. die kleinste und gr6éBte Wurzel der Gleichung 
aby? — (1 —6 (0, + Qo) —401) y +o, 0; =0 (4.11) 
bedeuten (y** = + oo im Falle ab=0). Wegen oj +0 ist 


bot’ <1. (4.12) 
Fiir die Elemente 


sind die Forderungen (3.23) erfiillt, und ebenso gilt damit (3.25), soweit die vor- 
kommenden Gr6Ben erklart sind. Man beweist dies mit vollstandiger Induktion. 
Das Verfahren (3.22) lautet damit 


1 rn [| @ 
Ont+1 — 97 Lo 1— bo), (On 00) O71 ( ) (4.13) 


Sind die Voraussetzungen (4.8), (4.9) erfiillt, so ist diese Iteration nach Zusatz 2 
(mit t = 0*) unbeschrankt durchfiihrbar, denn wegen (4.12) sind alle &,(< @, < @*) 
aus %. Da die Folge e, Monoton und beschrankt ist, konvergiert sie. Ihr Grenz- 
element ist @*, was man hier am einfachsten direkt aus (4.13) folgert. 

5.5. Eindeutigkeitsaussagen fiir die Aufgabe (4.1), (4.2) erhalt man bequemer 
mit dem vereinfachten Newtonschen Verfahren fiir den Operator (4.4). Dieses 


/ 


lautet mit “(1 (w= O40, 2.02.3) 


ui eed — ua), 
d. h. 
Vere Mag A Gy 1 et On = Pols Lo, = 1 “=O 4-2 ..)5 
stimmt also fiir »=0 mit (4.3) tiberein. Die zugehdrigen Vergleichsverfahren 
(3.22), (3.26) haben [unter der Voraussetzung (4.8)] die Form 


aioe = 05) oof 
omar cot +795] mit y,— (=H (= 0,1,2,...), 


On = 0, bzw. =T, gesetzt. Fiir die Zahlen y, erfolgt hieraus 


Ynti=H(Yn) (w= 0,4, 2, ..-) (4.14) 
mit 1 , a” , ”r 
hy) = aie. (6,0, + (bo, +40;)y+ aby’). 
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Aus der Konvergenz der Folge y, ergibt sich die der Folge w, und umgekehrt. 
Die Gleichung y=/(y) ist [unter der Voraussetzung (4.8)] mit (4.11) identisch. 
Sind also die Voraussetzungen (4.8), (4.9) erfiillt, so besitzt diese Gleichung die 
beiden Lésungen y* und y** (=> 0), und die Folge y, (4.14) konvergiert gegen y*, 
wenn 0<7)<y** ist. Der Grenzwert der Folge w, ist dann g*. Fir y= Qo 
wird y»=0. Die Folge (w, =) @, des ersten Vergleichsverfahrens konvergiert 
also gegen 9*. Die Folge (w, =) t, des zweiten Verfahrens hat denselben Grenz- 
wert o*, wenn 


On, — comeles Ly (4.15) 


Sty 
Nach der Eindeutigkeitsaussage von Satz 2 gilt daher fiir eine eventuell vor- 
handene weitere Lésung u** += u*: 


aem* — ug | |u**”"| = (1 =F, 00) eet 


denn andernfalls wiirde mit den Zahlen (to — @o) = ||w**’— uo||, T =|u**"'| die 
Ungleichung (4.15) erfiillt sein, die Folge 1,, also gegen o* konvergieren, und 
auBerdem (3.27) fiir w= u**=+ u* gelten, was nicht moglich ist. 

4.6. Wir fassen das aus diesen Uberlegungen mit Satz 2 fiir die Aufgabe (4.1), 
(4.2) und das Verfahren (4.3) folgende Ergebnis zusammen. 


Satz 3. p und gy seien Elemente mt lAw+0, lAg +0 und ly=lg)=1 
(z. B. auch p=). Es gebe Konstanten a,b derart, dag 


lA 
49 — 74244 salle 


; Tay 4elsellel fir lp=0 (4.16) 


st, und mut diesen Konstanten gelie 


1— b||y — poli > 0 
und 


es y : si PO set alls #10 10) 
J7=Fle—v— VoUTe—val + low) ~ Vall {= 9 oe ne an 


Dann konvergieren die Folgen p, und my, gegen ein Element g* bzw. eine 
Zahl w* 4-0, welche die Gleichungen 


1 
erfiillen; es gilt die Fehlerabschdtzung 


“|le*— || Say*, — |u*—p,|<dy* (4.19) 
und allgemein 


lle*— ell Se —e,,  |u*—w|Set”—0f (n =0,1,2,...), 
worn y* die kleinste Lésung der Gleichung 


ee a 
aoe" (11 — Poll + Ilo — vl) — 4| tl] y + |e. — oll lel = 0" (4.20) 
bedeutet und o*, 0, durch (4.10) bzw. (4.13) definiert sind. 
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Gibt es 1m Falle ab==0 eine weitere Lésung p**, w** der Gleichungen (4.18), 
so gilt fiir diese 


llo** — poll |m**| = (1 — b|| Po — ll) »* (4.21) 


mut der zweiten Wurzel y** der Gleichung (4.20). Im Falle ab=0 existiert keine 
weitere Losung. 

Am einfachsten und genauesten werden diese Ergebnisse fiir y=q,). Damit 
ist ||@) —y||=0 und (4.16) geht tiber in 


|A~—(CAgG) gill <allell, [ml |24e| <ollpl| fir Fp=0. (4.22) 


Jedoch ist es fir die numerische Rechnung zweckmaBig, auch den Fall p+ q, 
zuzulassen. Wenn man wahrend der Iteration mehrmals eine Fehlerabschatzung 
durchfiihren will, braucht man die Konstanten a und b dann nur einmal zu 
berechnen. 

§ 5. Spezialisierung auf Matrizen 


5.1. Der Spezialfall einer Matrizeneigenwertaufgabe (4.1), (4.2) fiir N-dimen- 
sionale Vektoren = (qj) mit einer N x N-Matrix A =(a;;) und einem Zeilen- 
vektor /= (/;)) soll noch etwas ausfiihrlicher behandelt werden. In diesem Falle 
ist 


mit der Matrix 
=) Aa a ees 
B=A—f(lA) und 8B @ Ay. 


Fir y=q, wird B= 9,. Die Abschatzung dieser Matrix B ist der wesentlichste 
Teil der Anwendung des Satzes 3. Sie erfordert am meisten Rechenarbeit und 
beeinfluBt die Fehlerabschatzung am starksten. Insbesondere muB a|4,| <1 er- 
reicht werden. Der Umfang der Rechenarbeit und die GréBe der Konstanten a 
hangen stark von der Wahl des Vektors / und der Norm |/g|| ab. Zur Bildung 
des Produktes f(/A) sind jedoch bei bekanntem f immer héchstens 2 N? Multi. 
plikationen erforderlich. 
Als besonders einfaches Beispiel wahlen wir fiir / den Einheitsvektor 


ba emit ve yia=i0,; (A ak LN) (ey) 
und einer geeigneten festen Nummer k (1 < k < N). Die ,,Normierungsbedingung” 


(4.2) wird dann 
Py = 1, 


und die Elemente der Matrix B sind leicht zu berechnen: 


b= a,,— Boy Mj 
Die k-te Zeile und Spalte von B brauchen nicht ermittelt zu werden. Die k-te 
Zeile enthalt nur Nullen und die k-te Spalte wird in der Abschatzung (4.16) 
nicht benétigt. Im Beispiel 2 (in 5.3) sind zur Berechnung des Produktes 6 (/ A) 
nur N—1 (=7) Multiplikationen mit der Zahl 2 erforderlich, da die a, ;(7 ==) 
nur Werte 0, 1 oder 2 annehmen. 
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Wahlt man als Norm | 
le] = mex (5.2) 


mit gegebenen Gewichten d;>0, so wird fiir den Vektor / aus (5.1) 


1<isNn = 
sas j+h 
1 he 1 
Paar elie Gee fiir v=o), (5.3) 
ih 
wobei C = (c,) =D4BD und D= (a; 6,,), also 
Ga 
Cig = a, 88: 


ist. 
Im Falle einer hermetischen Matrix A erhalt man mit p=, 1= Go: 


ee 
B=A—— : 
Ih Ya 
Zur Berechnung von B sind dann N? Multiplikationen erforderlich. Ware exact 


Ag=— @,, so wiirde B die gleichen Eigenvektoren und Eigenwerte besitzen 
[en 


wie A bis auf den Eigenwert 1/u,, an dessen Stelle B den Eigenwert 0 hatte. 


Im allgemeinen konvergiert die Folge 1/u,, gegen den dem Betrage nach gréBten 
Eigenwert 4“, und zwar um so schneller, je kleiner das Quotient |A®/A@| mit 
dem betragsmaBbig zweitgréBten Eigenwert A® ist. Kennt man die Lage der 
Eigenwerte ungefahr, so laBt sich dieser Quotient oft mit Hilfe einer Spektral- 
verschiebung verkleinern. Bei der Berechnung der Konstanten a kann sich er- 
geben (siehe Beispiel 2), daB mit einer Matrix A —aE statt A ein kleinerer Wert 
a|,| zu erhalten ware, offenbar ein Hinweis dafiir, daB eine solche Verschiebung 
des Spektrums um « auch fiir die Iteration giinstiger gewesen ware. 

Man wird dann nachtraglich noch einen Iterationsschritt mit der Matrix 
A—=A—aE durchfiihren. Beginnt man dabei mit Po=Yo, SO erhalt man 


Pe TAG, TAQ eee (5.4) 
Rene SBR TO Ge ea Le —o 
Py = 4 AG = Cae (A Gy —% Mp) = Bain (5.5) 


Die zugehérige Matrix B hatte z. B. fiir den Vektor / in (5.1) und p=qp die 
Elemente 


65; = (45; —Piay%e;) —%06;; fir 7,7-k. (5.6) 


Die Werte @,(;, werden zwar etwas besser sein als die P1(;), aber sich von diesen 
im allgemeinen nur in den letzten Stellen etwas unterscheiden, welche Verande- 
rung durch geeignete Abrundung beriicksichtigt werden kann. B braucht also 


nicht neu berechnet zu werden (siehe Beispiel 2), abgesehen von der Subtraktion 
der Zahl « in der Hauptdiagonalen. 
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5.2. Beispiel 1. Am Beispiel der stark unsymmetrischen Matrix 


jhe at: 
7, eA VO i Pe 
Cuter 2 


wird erprobt, welchen EinfluB die Wahl des Vektors / hat, und zwar werden die 
Einheitsvektoren /=e*(k=1, 2,3) und y= q, benutzt. Die Norm sei (5.2) mit 
d;=1 (¢=1, 2,3). Tabelle 1 und 2 enthalten die Ergebnisse. 


Tabelle 1 (Aujgabe 1). (Die Vektoren yy unterscheiden sich nur um einen Faktor) 


v Po M1 Pi- Po B 9: (1) 

‘Loa tl 1 6) * * * * 

O | 0,694 403} 0,694 568 0,000165 te — 0,389136 | — 0,083 704 | 0,422840 

0 | 0,420683) 0,420660) — 0,000023 “ 0,158 680 0,738020 | 0,896 700 

O | 1,440085] 1,439742) — 0,000 343 | —0,439742 0,120 516 | 0,560258 

Aad 1 6) Hy * * * 
0,605 820} 0,605 642; — 0,000178 | — 0,605 642 ss 0,788 716 | 1,394358 

O | 2,377082) 2,377 215 0,000 133 1 SOP OS ‘3 1,377215 

O | 1,650653] 1,651138 0,000485 1 —0;'654 138 5 1,651138 


Tabelle 2. (In der Spalte fiiy Aufgabe 2 beziehen sich die Werte auf die verschobene 
Matrix A=A—2E und die Norm mit den Gewichten d;= A;) 


Aufgabe 1 


Aufgabe 2 


k=2 

My 0,273 908 0,273 843 0,273 923 0,10006 
Ile Pol| 0,000 165 0,000 343 0,000485 0,000 75 

0,896 700 1,394358 1,651 138 6,803 
3 1,369 540 0,821 529 0,273 923 1,13401 
1—Yb || %—%|| — Va] 4 0,489 374 0,365 288 0,315947 0,145 80 
Fehler- {|| y*— i= 0,000 054 0,000 213 0,000 401 0,001 60 
abschtzg.| |u* —44| < 0,000 082 0,000 125 0,000 066 0,000 26 

ee — is 0,000 01 5 0,000 033 0,000045 — 

Exact { fey 0,000017 0,000048 0,000032 a 


5.3. Beispiel 2. An diesem Beispiel wird unter anderem der Einflu8 der Norm 
|| || erklart. Die numerischen Ergebnisse findet man in den Tabellen 2 und 3. 
Zur Lésung der Eigenwertaufgabe 


aU 


aU 


One 


dy? 


+ AU=0 


im Gebiet %, 


wird das gewohnliche Differenzenverfahren verwendet (Abb. 4). 


rungswerte g,,) in den Gitterpunkten P;(7=1, 2, ... 


U = 0 auf dessen Rande J’ 


Fiir die Nahe- 
, 8) des in Abb. 5 skizzierten 


Gitters erhalt man leicht unmittelbar [23] ein Gleichungssystem A gy =Aqg mit 


der in Tabelle 3 angegebenen Matrix A. 
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Es sei J der Einheitsvektor /=¢!. Der Vektor g» wurde von einer mit groBerer 
Maschenweite berechneten Naherung ausgehend nach 3 Iterationsschritten er- 
halten [23]. Mit p =p ergibt sich die Matrix B in einfacher Weise. Die Zeilen 
summen @;(1) =>)’ |b;,| dieser Matrix zeigen, daB fir die Norm (5.2) mit d;=1 


i] c se : 
die Konvergenzbedingung a|j,| <1 wegen || © 7y nicht erfullt ist. 


Abb. 4. Grundgebiet der Randwertaufgabe Abb. 5. Gitter fiir das Differenzenverfahren 
Wir verschieben das Spektrum. Fiir A =A—2E erhilt man die neuen Werte 
fix, G, aus (5.4), (5.5). G, unterscheidet sich von g, so wenig, daf B von B 
nur in der vierten und fiinften Dezimale abweicht, abgesehen davon, daB die 
Hauptdiagonalelemente von B um die Zahl 2 zu verringern sind [vgl. (5.6)]. 
Ferner verwenden wir, da die Summen der Zeilen 2 und 3 besonders groB 
sind, die Norm (5.2) mit d;= A, 


Pe pe TF 0,65. Ade consis 

Tabelle 3 

A l’ we =e Pe ere P1()— Yo (3) 
Ai 
l 
4 tod 002/20 22 1 1 1 | 4 ) 0 

DEAN OROR OFA MAO) O 0,924 80 0,92470 | 0,92468 — 0,00012 — 0,00007 
DLORAP OTA ONO 2 O 0,924 83 0,924 71 0,924 68 —0,00015 — 0,00009 
OOM R252 OROmO O 0,31643 0,31622 0,31617 — 0,000 26 — 0,000 26 
DOs yr 30 he le Oy Il O 0,655 00 0,65468 0,654 61 — 0,000 39 —0,00039 
2D PW CY Gr ae 1 @ O 0,651 59 0,65096 0,65084 —0,00075 — 0,00075 
2, 2) (0) Oy Oy jl 4h Gh O 0,696 14 0,696 62 0,696 71 0,000 57 0,000 57 
Ab (OD)! PY LO DO 2 ah O 1,06901 1,069 61 1,069 73 0,000 72 0,000 72 


* % # ¥ * & * * 
OV Oe® 7 WS: ro: = 
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Die damit nach (5.3) ESS Zahlen 8,(A,) besitzen ein Maximum 4, das 


zusammen mit 6 = ind, a,,A 


i,| die Konvergenzbedingung (4.17) 


at) 
erfiillt. Die Free Dniane. cutee in Tabelle 2 und 3. 


§ 6. Stérungsrechnung 


Die praktische Durchfiithrung des Newtonschen Verfahrens (0.2) kann manch- 
mal schwierig sein, weil darin [G(,)] + auftritt. Bei dem nun beschriebenen Ver- 
fahren wird dieser Operator durch [Gow]? ersetzt, wobei Gy ein G_ ,,benach- 
barter“‘, mdglichst ,,einfacher“‘ Operator ist. 


6.1. R, S und ¥ seien wie in § 3. Zur Lésung einer Gleichung 


Gu=8 mit G=G,+6G, (6.1) 
wird das Iterationsverfahren 
Un+y = Uy, 4 [Go (Nie Gu, (n a 0, ty 2; oD 5) (6.2) 


angesetzt. Gp) bedeute dabei einen Operator, der alle in 3.1 von G geforderten 
Eigenschaften besitzt. Die zugehdrigen Operatoren H, H und die zugehorige 
Abschatzungsfunktion F werden durch einen Index 0 gekennzeichnet, die iibrigen 
Bezeichnungen 9, %, %, ... — welche sich nun also alle auf Gy beziehen — weiter 
ohne Index verwendet. 


Der Operator G, bilde D in © ab, wir setzen H, = [G,)| 1G,. Es gebe einen 
® in 9 abbildenden Operator H, derart, daB gilt: 

a) H,o =o, 
b) 0S 4, (9 +0) —H,o SH, (0' + 0’) — Ayo’ 

fir o<oSo’ CD, 0<050' Ed, 
c) v[H,w — H,v] <H, (v[w — 0] + »[v — w])— FA, (» [v — v)) 

fir u,v€D, v[u—v|]+r[v—wle®. (6.3) 
6.2. T(u, x) sei fir wED, «ES durch 

T(u, 2) =u — [Gq]? Gu 


erklart. Dann sind die Verfahren (0.1) und (6.2) aquivalent. AuBerdem ist jede 
der Bedingungen (0.4) mit Gu*=@ gleichbedeutend, falls »[x, —w]<@ E® ist 
(Beweis wie in 3.2). 
Es ist 
T(u, x) —T(v, y) = [Gow] 4{— Gou + Gov — Gow (% — 4) + Gog (z — v)}+ 
+- [Gow] t{— G4+G,r}, z=T,y). 
Die Norm » des ersten Summanden laBt sich nach 3.3 durch Jj [o*] mit den Werten 
(3.13) abschatzen. Die Norm des zweiten Summanden ist 
v[R(x) (— Hye + Hy») ]<(E — Ao) (Ai (e+ 0°) — Hy 0°) = Ale’). 


Also hat man 


y[T(u, x) —Tv, y) | <Flo)] + Alo] =F le’). 
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F[o'] hat alle in 2.1 geforderten Eigenschaften, (2.5) gilt fiir den (2.9) er- 
fiillenden Operator 


F(o, )=0 + [E—Ho ey]? ((E—Foey] + 09—H o,—0 | H 9), 6 (#1, UM) SO. 


os G=G6,4+46, (6.4) 
a Coie Hy ] % + Go — Hees 0 +Hye, 
Ge = —H,o,+ Ko 
erhalten die Vergleichsverfahren (2.6), (2.7) die (6.2) entsprechende Form 
Ont = On — Ga 2G Op DLWe Gera ic, el 2G Te i =0,4, 227-08 


Aus 0<tE® und Gr<o folgt wie in 3.4 fiir jedes BER: t2eT(t, &). Damit 
gilt: 

Satz und Zusatz 4. Sind die in 6.1 genannten Voraussetzungen erfillt, so 
gelten wortlich die Aussagen des Satzes 2 und des Zusatzes 2, wenn darin G und 
G die durch (6.1) bzw. (6.4) definierten Operator bedeuten und man in (3.22), (3.26) 
und (3.24) G’ und G’ durch Go bzw. Go ersetzt. 


§ 7. Polynome als Majoranten beim Verfahren u,.,;= Tu, 


Fiir das einfachere Verfahren u,,,—= Tu, werden Vergleichsverfahren herge- 
leitet, bei denen To ein Polynom in @ ist. 

# und XN seien Banachraume, R sei Yt-normiert und ||w|| = ||y [w]|| (siehe § 1). 
T bedeute einen Operator, welcher eine offene konvexe Teilmenge DCR in KR 
abbildet. Die (auf % definierten) Ableitungen Ti) (G=4, 2,3, ...) werden’ de- 
finiert als Operatoren, welche die Menge ft’ der 7-tupel h,, fa, ...,h; mit h,EeR 
(unabhangig von der Reihenfolge der h,) linear beziiglich jeder Komponente h, 
in }t abbilden und den Beziehungen 


jim | 258 — Te? — TB h|| =o 
geniigen. 

Mit L; bezeichnen wir Operatoren, welche die Menge der i-tupel 0, 0?, ..., 0% 
mit ofEN (unabhangig von der Reihenfolge) in 9 abbilden, linear beziiglich 
jeder pombe positiv (L; 0! o®... oe’ So fiir o1, o?,..., e’ 0) und beschrankt 


sind (IZ: CORO tee etl le?l|..- |le’]|). Eimen Se Operator nennen 
wir eine St- ae von 7, falls 


y [73 hy hy... ASL, vy [hy] v [he]... v[h,] 
fir alle Ay, hy,...,h,ER gilt. Ti) hei®t dann durch L; N-beschrankt. 


‘Satz 5. eee der Operator T fiir jedes wCD (auf KR definierte) Ableitungen 
Dee te) —1), die fiir ein festes u=wCD durch L; N-beschrankt sind 38, 


8 Diese Voraussetzung entfallt also fiir m—=1. 
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und gentigt die (m—1)-te Ableitung einer Lipschitzbedingung 


vl(Te ie “)— Live SO) Dy Teas My, 1) SL» (U—v] » [Ay] v [he]... ¥ Pm —1] 
fir u,veD, her 
(m= 1), so gilt 
y[Tu — Tr] <F[y[u —v], »[u— vw], v[v —w]] 
mit 


m 4 1 on = ah ae 
wc ol > ie ey Ce) = Fle—v,¢,6l=lo—16, 
‘= k=1 
To=>'-L;(o')+9 (HEN, beliebig). (7.4) 


Fir 0=0,+ 0 a + L;(@o)‘ und » [u, — uo] <o, erfiillt T die Forderung (2.9). 

i=1 

Fiir den Zahlenraum %t wurde dieser Satz in [25] § 3 bewiesen. Dieser Be- 
weis gilt auch in dem hier behandelten allgemeineren Fall, wenn man statt || || 
die Norm y[u] verwendet und einige geringe Anderungen vornimmt. Zum Bei- 
spiel ist darauf zu achten, daB der Ausdruck L;o!0? ... o* seinen Sinn verliert, 
wenn man L, weglaBt, was bei Zahlen 0’ geschehen konnte. Der Beweis sei 
daher hier nicht noch einmal wiederholt. 

Das an der Stelle w angesetzte vereinfachte Newtonsche Verfahren ist von 
der Form u,,,=Tu, mit Tj) =O. Unter den Voraussetzungen des Satzes 5 
(z. B. fiir m= 2) wiirde man also in (7.1) eimen Vergleichsoperator T erhalten, 
bei dem das in @ lineare Glied fehlt. 

Ist G(u) eine gewohnliche Funktion im Raum der reellen oder komplexen 
Zahlen mit geeigneten Differenzierbarkeitseigenschaften, so gilt fiir das Itera- 
tionsverfahren 


Unt1 = = T(u,,) 
mit " )) 
Ti) = 6 — (era + SE Cw) +S ree CWP (7.2) 


zur Lésung der Gleichung G(u) =0: 
DOG ed aN) eX), 


so da8 dafiir in (7.1) auch das quadratische Glied fehlen wiirde. 
Das Verfahren (7.2) stimmt fiir w =, im ersten Schritt mit dem verbesserten 
Newtonschen Verfahren 
— , — Gn) 1 [G(Un) PG" (Un) 
Mat = Un Glug) 2 [Wind P 73) 


iiberein, so daB die Fehlerabschatzung fiir (7.2) auch fiir (7.3) verwendet werden 
kann. 
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